hdu3068

http://www.elijahqi.win/2017/07/16/hdu3068/
最长回文

Time Limit: 4000/2000 MS (Java/Others) Memory Limit: 32768/32768 K (Java/Others)
Total Submission(s): 21988 Accepted Submission(s): 7968

Problem Description
给出一个只由小写英文字符a,b,c…y,z组成的字符串S,求S中最长回文串的长度.
回文就是正反读都是一样的字符串,如aba, abba等

Input
输入有多组case,不超过120组,每组输入为一行小写英文字符a,b,c…y,z组成的字符串S
两组case之间由空行隔开(该空行不用处理)
字符串长度len <= 110000

Output
每一行一个整数x,对应一组case,表示该组case的字符串中所包含的最长回文长度.

Sample Input
aaaa

abab

Sample Output
4
3

Source
2009 Multi-University Training Contest 16 - Host by NIT

Recommend
lcy | We have carefully selected several similar problems for you: 1358 3336 3065 3746 3067

复习一遍manacher 贴了下别人的讲解

当 mx - i > P[j] 的时候，以S[j]为中心的回文子串包含在以S[id]为中心的回文子串中，由于 i 和 j 对称，以S[i]为中心的回文子串必然包含在以S[id]为中心的回文子串中，所以必有 P[i] = P[j]，见下图。

当 P[j] > mx - i 的时候，以S[j]为中心的回文子串不完全包含于以S[id]为中心的回文子串中，但是基于对称性可知，下图中两个绿框所包围的部分是相同的，也就是说以S[i]为中心的回文子串，其向右至少会扩张到mx的位置，也就是说 P[i] >= mx - i。至于mx之后的部分是否对称，就只能一个一个匹配了。

对于 mx <= i 的情况，无法对 P[i]做更多的假设，只能P[i] = 1，然后再去匹配了

下面给出原文，进一步解释算法为线性的原因

#include<cstdio>
#include<cstring>
#define N 120000
char str1[N],a[N<<1];
int p[N<<1],id,mx;
inline int min(int a,int b){
    return a<b?a:b;
}
int main(){
    freopen("hdu3068.in","r",stdin);
    freopen("hdu3068.out","w",stdout);
    int ans;
    while (~scanf("%s",str1)){
        ans=0;
        int n=strlen(str1);a[0]='0';
        for (int i=1;i<=n;++i){
            a[2*i]=str1[i-1];a[2*i-1]='#';
        }
        a[2*n+1]='#';
        //printf("%s",a);
        n=2*n+1;
        mx=id=0;
        for (int i=1;i<=n;++i){
            if (i<mx){
                p[i]=min(p[2*id-i],mx-i);
            }else p[i]=1;
            while (a[i+p[i]]==a[i-p[i]]&&(i-p[i]>0)&&(i+p[i]<=n)) p[i]++;
            if (i+p[i]>mx) mx=i+p[i],id=i;
            if (ans<p[i]) ans=p[i];
        }
        printf("%d\n",ans-1);
    }
    return 0;
}