hdu 3068最长回文

使用o(n)的manacher算法，详见代码

#include<cstdio>  
#include<iostream>  
#include<algorithm>
#include<math.h> 
#include<string.h>  
#include<vector> 
#include<queue>
#include<iterator>
#include<vector>
#include<set>
#define dinf 0x3f3f3f3f
typedef long long ll;
//const int Max=(1<<16)+10;
using namespace std;
#define SIZE 1000005


char tmp[SIZE],s[SIZE];//tmp存变换后的串，s存原串 
int p[SIZE];//存以mid为中心的回文串的长度 


void init()
{
    int slen=strlen(s);
    tmp[0]='$';        //头标记防止越界 
    for(int i=1,j=0;i<=2*slen+1;i+=2)
    {
        tmp[i]='#';
        tmp[i+1]=s[j++];
    }
    tmp[2*slen+2]='#';//尾标记 
}

int manacher()
{
    int mid=0,mx=0,ans=0;//mx为当前计算回文串最右边字符位置最大值 
    int tlen=strlen(s)*2+1;
    for(int i=1;i<=tlen;i++)
    {
        if(mx>i)
            p[i]=min(mx-i,p[2*mid-i]);
        else        //i>mx从头开始匹配 
            p[i]=1;
        while(tmp[i-p[i]]==tmp[i+p[i]])
            p[i]++;
        if(p[i]+i>mx)//计算得到的右端点值大于mx则进行更新 
        {
            mx=p[i]+i;
            mid=i;
        }
        ans=max(ans,p[i]);
    }
    return ans-1;//p[i]-1即最长回文串的长度 
}

int main()
{
    while(scanf("%s",&s)!=EOF)
    {
        init();
        printf("%d
",manacher());
    }
    return 0;
}