hiho一下~week_1:最长回文子串

题目意思：
题目链接：http://hihocoder.com/contest/hiho1/problem/1
给出一个字符串（n<=10^6），求出这个字符串的最长回文字串。输入有多组数据，只需要输出对应的结果即可。
样例输入
3
abababa
aaaabaa
acacdas
样例输出
7
5
3
解析：当我们看到n<=10^6时，证明一般的模拟算法肯定是过不去的，对于时间复杂度O(n^3),O(n^2)的算法肯定是要TL的，因此需要提到Manacher算法，在O(n)时间求出最长回文字串。看一下小hi的提示：
提示三：
小Ho答道：“我想想，如果以第5个字符为中心的最长回文子串的长度是5的话，这就告诉了我[3, 7]这一段是一个回文子串，所以呢？”

小Hi继续提示道：“假设这时候你想要计算以第6个字符为中心的最长回文子串的长度，你有没有什么已知的信息了？”

小Ho边想边说道：“唔，首先第6个字符和第4个字符是一样的，第7个字符和第3个字符是一样的，而第5个字符本身就肯定和第5个字符一样，那么如果[3, 5]这一段是回文子串的话，那么[5, 7]这一段肯定也是回文子串。
Manacher算法：引用文章

AC代码：

/**
 *Author: xiaoran
 *Solution:
 *
 */
#include<iostream>
#include<cstdio>
#include<string>
#include<cstring>
#include<cmath>
#include<algorithm>
#define LL long long
using namespace std;
const int MAXN=1000005;
char s[MAXN];

bool huiWen1(char *s,int low,int high){//枚举起始位置
    while(low<high){
        if(s[low]!=s[high]) return false;//不是回文
        else low++,high--;
    }
    return true;//是回文
}
int huiWen2(char *s,int mid){//枚举中心位置mid
    int len=strlen(s);
    //奇数
    int max1=1;
    int k=max(mid,len-mid-1);
    for(int i=1;i<=k;i++){
        if(s[mid-i]==s[mid+i]) max1+=2;
        else break;
    }

    //偶数
    int max2=0;
    k=max(mid,len-mid-1);
    for(int i=0;i<k;i++){
        if(s[mid-i]==s[mid+i+1]) max2+=2;
        else break;
    }

    return max(max1,max2);
}

char Ma[MAXN*2];
int Mp[MAXN*2];
//Mp[i]表示第i个字符的回文半径,即rad[i]尽可能大
//且满足Ma[i-Mp[i],i-1]=Ma[i+1,i+Mp[i]]

//Manacher算法。时间O（n），空间O（n）最长回文串
int Manacher(char *s,int len){
    int k=0;
    //初始化字符串,将所有字符串转化为奇数长度
    Ma[k++]='$';
    Ma[k++]='#';
    for(int i=0;i<len;i++){
        Ma[k++]=s[i];
        Ma[k++]='#';
    }
    Ma[k]='\0';//字符串结束

    //其中id表示最大回文子串中心的位置，mx则为id+P[id]，也就是最大回文子串的边界
    int mx=0,id=0;
    for(int i=0;i<k;i++){
        Mp[i]=mx>i?min(Mp[2*id-i],mx-i):1;
        while(Ma[i+Mp[i]]==Ma[i-Mp[i]]) Mp[i]++;
        if(i+Mp[i]>mx){
            mx=i+Mp[i];
            id=i;
        }
    }
}

int main()
{

    //freopen("E:/input.txt","r",stdin);
    //freopen("E:/output.txt","w",stdout);
    int t;
    scanf("%d",&t);
    while(t--){
        scanf("%s",s);
        int len=strlen(s);
        int lenm=1,maxlen=len;
        int ok=0;
        /**超时O(n^3)
        for(int i=0;i<len;i++){
            for(int j=i+1;j<len;j++){
                if(huiWen1(s,i,j)&&lenm<j-i+1){
                    lenm=j-i+1;
                }
            }
        }
        cout<<lenm<<endl;
        **/

        /**
        //从大到小枚举长度
        while(!ok){//还是超时，O(n^3)系数可能会小一些
            for(int i=0;i+maxlen<=len;i++){
                if(huiWen1(s,i,i+maxlen-1)){
                    ok=1;//找发最长回文
                    break;
                }
            }
            --maxlen;//遍历下一个长度
        }
        cout<<maxlen+1<<endl;
        **/

        /**
        //枚举中心位置（考虑奇偶）O(n^2)超时
        for(int i=1;i<len-1;i++){
            lenm=max(lenm,huiWen2(s,i));
        }
        cout<<lenm<<endl;
        **/
        Manacher(s,len);
        int ans=0;
        for(int i=0;i<2*len+2;i++){
            ans=max(ans,Mp[i]-1);
        }
        printf("%d\n",ans);
    }
    return 0;
}

注：参考文章和kuangbin的代码。（小猫的重学，请大神见谅！）