【Manacher算法】poj3974 Palindrome

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

模板题,Code 附带注释:

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