Manacher-模版题poj3974 hdu3068

来源：互联网发布：aes算法原理编辑：程序博客网时间：2024/05/22 16:59

当i>mx时，就要重头开始匹配啦

其余看代码的注释好好理解啊

//Manacher模版题poj3974#include<cstdio>#include<cstring>#include<cstdlib>#include<iostream>#include<algorithm>using namespace std;#define maxn 2001000int g[maxn];char ss[maxn],s[maxn];//ss为原串，s是ss相邻字符间加一个奇怪字符的串//g[i]表示以字符s[i]为中心的最长回文长度int mymax(int x,int y){return (x>y)?x:y;}int mymin(int x,int y){return (x<y)?x:y;}int main(){//freopen("a.in","r",stdin);//freopen("a.out","w",stdout);int len,l,ans,mx,p,Case=0;while (1){scanf("%s",ss+1);if (strcmp(ss+1,"END")==0) break;l=strlen(ss+1);len=0;s[len]='|';for(int i=1;i<=l;i++){s[++len]=ss[i];s[++len]='|';}mx=0;p=0;ans=0;g[0]=1;//p为当前找到的最长回文串的（中心）位置//mx表示目前到达的最远位置，即最长回文串的最右端for(int i=1;i<=len;i++){if (mx>i) g[i]=mymin(g[2*p-i],mx-i);//当mx>i时,2*p-i是i关于p的对称点,它的g值是已经算好了的,//如果i+g[2*p-i]还在mx的范围里，那么g[i]=g[2*p-i]//否则就是说,目前可以确定的匹配好的长度为mx-i,从mx-i开始匹配else g[i]=1;//不然就从头开始匹配while (i+g[i]<=len && s[i-g[i]]==s[i+g[i]]) g[i]++;//匹配ingif (i+g[i]>mx) {mx=i+g[i];p=i;}//更新mx,pans=mymax(ans,g[i]);}printf("Case %d: %d\n",++Case,ans-1);//g[i]中的最大值-1即为原串的最长回文子串长度}return 0;}

0 0