【Manacher算法】poj3974 Palindrome

Manacher算法教程：http://blog.csdn.net/ggggiqnypgjg/article/details/6645824

模板题,Code 附带注释:

 1 #include<cstdio>
 2 #include<algorithm>
 3 #include<cstring>
 4 using namespace std;
 5 char b[10000001],a[10000001];
 6 char tmp[10000001];
 7 int n,f[10000001],zu;//f[i]表示插入一堆#后，以i为中心的最长回文子串の半径长度，
 8                 //故其减1后就是原串的最长回文子串的答案
 9                 //也是原串中以i为开端的最长回文串の长度。
10 int id,maxid,ans;
11 int main()
12 {
13     while(1)
14       {
15           scanf("%s",tmp);
16           if(tmp[0]==‘E‘)
17             break;
18           zu++;
19           n=strlen(tmp);
20           memset(b,0,sizeof(b));
21           memset(a,0,sizeof(a));
22           memset(f,0,sizeof(f));
23           ans=id=maxid=0;
24           for(int i=1;i<=n;i++)
25           b[i]=tmp[i-1];
26         a[1]=‘#‘;
27         for(int i=1;i<=n;i++)
28           {
29             a[i<<1]=b[i];
30              a[i<<1|1]=‘#‘;
31            }
32         a[0]=‘-‘;
33         a[(n+1)<<1]=‘+‘;
34         n=n<<1|1;
35         f[1]=1;
36         id=1;//用id这个变量记下取得这个最优maxid时的id值
37              //即右端扩展到maxid+1时，该回文串中心的位置
38         maxid=2;//maxid是曾经扫描到的回文串中，匹配到的最远的位置+1
39         for(int i=2;i<=n;i++)
40           {
41              if(maxid>i)//算法核心：防止重复匹配
42                  f[i]=min(f[(id<<1)-i]//以 关于id的对称点 为中心的的最长回文串长
43                                    //因为，分别在id两侧的两半回文串是完全一样的
44                        ,maxid-i    // 但是，以id的对称点为中心的最长回文串有可能超出
45                                    //以id为中心的最长回文串的范围，所以，限制其无法超出
46                                    //此范围
47                                    ); //(id<<1)-i 为 i 关于 id 的对称点
48               else
49               f[i]=1;//否则f[i]=1
50               for(;a[i+f[i]]==a[i-f[i]];f[i]++);
51               ans=max(ans,f[i]-1);
52             if(f[i]+i>maxid)//注意是f[i]+i，不是f[i]+i-1，因为maxid是匹配到的最远位置+1
53               {
54                 maxid=f[i]+i;
55                 id=i;//若以i为中心时，回文串可以扩展到更远的地方，更新id
56               }
57           }
58         printf("Case %d: %d\n",zu,ans);
59       }
60
61     return 0;
62 }