HDOJ3613解题报告【拓展KMP】

题目地址：

　　http://acm.hdu.edu.cn/showproblem.php?pid=3613

题目概述：

　　给出一个字符串s以及每个小写字母的价值，现要求将s分为两个部分，对于每部分，如果是回文串则这个子串的价值为所有字母的价值之和，否则为零，请合理分开这个串使总价值最大。

大致思路：

　　EX_KMP。其实也是没有看懂啦。

　　不过看了下这位大神的BLOG：http://blog.csdn.net/dyx404514/article/details/41831947

　　迷迷糊糊的过了……

代码：

#include <iostream>
#include <cstdio>
#include <cstdlib>
#include <cmath>
#include <vector>
#include <ctime>
#include <map>
#include <stack>
#include <queue>
#include <cstring>
#include <algorithm>
using namespace std;

#define sacnf scanf
#define scnaf scanf
#define maxn  500010
#define maxm 26
#define inf 1061109567
#define Eps 0.00001
const double PI=acos(-1.0);
#define mod 7
#define MAXNUM 10000
void Swap(int &a,int &b) {int t=a;a=b;b=t;}
int Abs(int x) {return (x<0)?-x:x;}
typedef long long ll;
typedef unsigned int uint;

char s1[maxn],s2[maxn];
int extend1[maxn],extend2[maxn],Next[maxn];
int val[maxm],sum[maxn];

void Get_Next(char *S,int len)
{
    int po=0;
    Next[0]=len;
    while(po+1<len&&S[po]==S[po+1]) po++;
    Next[1]=po;po=1;
    for(int i=2;i<len;i++)
    {
        if(Next[i-po]+i<Next[po]+po) Next[i]=Next[i-po];
        else
        {
            int j=max(0,Next[po]+po-i);
            while(i+j<len&&S[j]==S[i+j]) j++;
            Next[i]=j;po=i;
        }
    }
}

void EKMP(char *S,char *T,int *extend,int len)
{
    Get_Next(s1,len);
    int po=0;
    while(po<len&&S[po]==T[po]) po++;
    extend[0]=po;po=0;
    for(int i=1;i<len;i++)
    {
        if(Next[i-po]+i<extend[po]+po) extend[i]=Next[i-po];
        else
        {
            int j=max(0,extend[po]+po-i);
            while(i+j<len&&S[i+j]==T[j]) j++;
            extend[i]=j;po=i;
        }
    }
}

int main()
{
    //freopen("data.in","r",stdin);
    //freopen("data.out","w",stdout);
    //clock_t st=clock();
    int T;scanf("%d",&T);
    while(T--)
    {
        for(int i=0;i<26;i++) scanf("%d",&val[i]);
        scanf("%s",s1);int len=strlen(s1);sum[0]=0;
        for(int i=0;i<len;i++)
        {
            sum[i+1]=sum[i]+val[s1[i]-'a'];
            s2[len-i-1]=s1[i];
        }

        EKMP(s2,s1,extend1,len);
        EKMP(s1,s2,extend2,len);

        int ans=0;
        for(int i=1;i<len;i++)
        {
            int temp=0;
            if(extend1[len-i]==i) temp+=sum[i];
            if(extend2[i]==len-i) temp+=sum[len]-sum[i];
            ans=max(ans,temp);
        }
        printf("%d
",ans);
    }
    //clock_t ed=clock();
    //printf("

Time Used : %.5lf Ms.
",(double)(ed-st)/CLOCKS_PER_SEC);
    return 0;
}