2342: [Shoi2011]双倍回文
Description
Input
输入分为两行,第一行为一个整数,表示字符串的长度,第二行有个连续的小写的英文字符,表示字符串的内容。
Output
输出文件只有一行,即:输入数据中字符串的最长双倍回文子串的长度,如果双倍回文子串不存在,则输出0。
Sample Input
16
ggabaabaabaaball
ggabaabaabaaball
Sample Output
12
HINT
N<=500000
第一次写manacher。。。
其实manacher还是很好理解的,就是已知的最长的回文串起点为id,然后对于第i位的字符和i与id的对称点j(j显然已经被计算过了)如果j的回文串包含在已知的最长回文串中,那么显然i的回文串也是包含在已知的当中,如果没有,就暴力向后匹配。。
然后对于这道题。。先用manacher,然后再枚举终点。。
这里有一个优化就是用已知的ans剪枝。。然后就可以过了。。。
1 #include<iostream> 2 #include<cstdlib> 3 #include<cmath> 4 #include<cstring> 5 #include<cstdio> 6 #include<algorithm> 7 #include<string> 8 #include<map> 9 #include<queue> 10 #include<vector> 11 #include<set> 12 #define inf 1000000000 13 #define maxn 500000+5 14 #define maxm 10000+5 15 #define eps 1e-10 16 #define ll long long 17 #define for0(i,n) for(int i=0;i<=(n);i++) 18 #define for1(i,n) for(int i=1;i<=(n);i++) 19 #define for2(i,x,y) for(int i=(x);i<=(y);i++) 20 #define for3(i,x,y) for(int i=(x);i>=(y);i--) 21 #define for4(i,x) for(int i=head[x],y=e[i].go;i;i=e[i].next,y=e[i].go) 22 using namespace std; 23 int n,p[maxn]; 24 char s[maxn],st[maxn]; 25 int read(){ 26 int x=0,f=1;char ch=getchar(); 27 while(ch<'0'||ch>'9'){if(ch=='-')f=-1;ch=getchar();} 28 while(ch>='0'&&ch<='9'){x=10*x+ch-'0';ch=getchar();} 29 return x*f; 30 } 31 int main(){ 32 //freopen("input.txt","r",stdin); 33 //freopen("output.txt","w",stdout); 34 int n=read();scanf("%s",st+1); 35 for1(i,n)s[i<<1]=st[i];n<<=1;n++; 36 for(int i=1;i<=n;i+=2)s[i]='#'; 37 int id=0,mx=0; 38 for(int i=1;i<=n;i+=2){ 39 if(mx>i)p[i]=min(p[2*id-i],mx-i); 40 while(i-p[i]>0&&i+p[i]<=n&&s[i-p[i]]==s[i+p[i]])p[i]++; 41 if(i+p[i]>mx)mx=i+p[i],id=i; 42 } 43 int ans=0; 44 for(int i=1;i<=n;i+=2){ 45 int j=(p[i]-1)>>1; 46 if(j&1)j--; 47 for(;j>0&&2*j>ans;j-=2) 48 if(p[i-j]-1>=j&&p[i+j]-1>=j)ans=max(ans,2*j); 49 } 50 printf("%d",ans); 51 return 0; 52 }