【Manacher算法】poj3974 Palindrome

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

模板题,Code 附带注释:

 #include<cstdio>
 #include<algorithm>
 #include<cstring>
 using namespace std;
 char b[],a[];
 char tmp[];
 int n,f[],zu;//f[i]表示插入一堆#后，以i为中心的最长回文子串の半径长度，
                 //故其减1后就是原串的最长回文子串的答案
                 //也是原串中以i为开端的最长回文串の长度。
 int id,maxid,ans;
 int main()
 {
     while()
       {
           scanf("%s",tmp);
           if(tmp[]=='E')
             break;
           zu++;
           n=strlen(tmp);
           memset(b,,sizeof(b));
           memset(a,,sizeof(a));
           memset(f,,sizeof(f));
           ans=id=maxid=;
           for(int i=;i<=n;i++)
           b[i]=tmp[i-];
         a[]='#';
         for(int i=;i<=n;i++)
           {
             a[i<<]=b[i];
              a[i<<|]='#';
            }
         a[]='-';
         a[(n+)<<]='+';
         n=n<<|;
         f[]=;
         id=;//用id这个变量记下取得这个最优maxid时的id值
              //即右端扩展到maxid+1时，该回文串中心的位置
         maxid=;//maxid是曾经扫描到的回文串中，匹配到的最远的位置+1
         for(int i=;i<=n;i++)
           {
              if(maxid>i)//算法核心：防止重复匹配
                  f[i]=min(f[(id<<)-i]//以 关于id的对称点 为中心的的最长回文串长
                                    //因为，分别在id两侧的两半回文串是完全一样的
                        ,maxid-i    // 但是，以id的对称点为中心的最长回文串有可能超出
                                    //以id为中心的最长回文串的范围，所以，限制其无法超出
                                    //此范围
                                    ); //(id<<1)-i 为 i 关于 id 的对称点
               else
               f[i]=;//否则f[i]=1
               for(;a[i+f[i]]==a[i-f[i]];f[i]++);
               ans=max(ans,f[i]-);
             if(f[i]+i>maxid)//注意是f[i]+i，不是f[i]+i-1，因为maxid是匹配到的最远位置+1
               {
                 maxid=f[i]+i;
                 id=i;//若以i为中心时，回文串可以扩展到更远的地方，更新id
               }
           }
         printf("Case %d: %d\n",zu,ans);
       }

     return ;
 }