本文版权归ljh2000和博客园共有,欢迎转载,但须保留此声明,并给出原文链接,谢谢合作。
本文作者:ljh2000
作者博客:http://www.cnblogs.com/ljh2000-jump/
转载请注明出处,侵权必究,保留最终解释权!
//It is made by ljh2000
#include <iostream>
#include <cstdlib>
#include <cstring>
#include <cstdio>
#include <cmath>
#include <algorithm>
using namespace std;
typedef long long LL;
const int MAXN = 200011;
int n,len,a[MAXN],p[MAXN],ans;
inline int getint(){
int w=0,q=0; char c=getchar(); while((c<‘0‘||c>‘9‘) && c!=‘-‘) c=getchar();
if(c==‘-‘) q=1,c=getchar(); while (c>=‘0‘&&c<=‘9‘) w=w*10+c-‘0‘,c=getchar(); return q?-w:w;
}
inline void manacher(){
//memset(p,0,sizeof(p));
for(int i=1;i<=n;i++) p[i]=0;
int maxR=0,id=0;
for(int i=1;i<=n;i++) {
if(i<maxR) p[i]=min(p[2*id-i],maxR-i);
else p[i]=1;
for(;i+p[i]<=n && a[i-p[i]]==a[i+p[i]];p[i]++) ;
if(i+p[i]>maxR) { maxR=i+p[i]; id=i; }
}
}
inline void work(){
int T=getint(); int Case=0;
while(T--) {
len=getint(); a[0]=-2; a[1]=-1; n=1; for(int i=1;i<=len;i++) { a[++n]=getint(); a[++n]=-1; }
manacher(); ans=1;
//+2的原因是对称中心只能取在区分符号上
for(int i=3;i<n/*!!!*/;i+=2) //枚举the second part
for(int j=ans;j<=p[i];j+=2) //枚举一个部分的长度,利用ans剪枝
if(p[i+j-1]>=j) //如果the third part的回文长度大于等于j,即j可以作为答案(part2和part3对称)
ans=j;
Case++; ans/=2; ans*=3;//只算了两个part的长度
printf("Case #%d: %d\n",Case,ans);
}
}
int main()
{
work();
return 0;
}
原文:http://www.cnblogs.com/ljh2000-jump/p/6266113.html