BZOJ 2565 最长双回文串（回文自动机）
题意
给一个长度为N的字符串S。对于一个字符串AB，如果A和B都是回文串，那么称AB是一个双回文串。求问S最长双回文子串的长度？
N <= 100000
题解
正反双向构造回文自动机，得到某一个点为结尾和开始的最长回文串长度，记为f[i],g[i];
答案就是f[i]+g[i+1]的最大值。
1 #include<iostream> 2 #include<cstring> 3 #include<cstdio> 4 #include<cmath> 5 #include<algorithm> 6 using namespace std; 7 const int N=100500; 8 int len[N],cnt[N],to[N][30],s[N],fail[N],last,tot,n; 9 int ans,L,p1[N],p2[N]; 10 char s1[N],s2[N]; 11 void init(){ 12 tot=1;n=0; 13 len[0]=0;len[1]=-1; 14 fail[0]=1;s[0]=-1;last=0; 15 memset(to,0,sizeof(to)); 16 } 17 void add(int c,int p[]){ 18 s[++n]=c; 19 int cur,now,tmp; 20 for(cur=last;s[n-len[cur]-1]!=c;cur=fail[cur]); 21 if(!to[cur][c]){ 22 tot++; 23 len[tot]=len[cur]+2;now=tot; 24 for(tmp=fail[cur];s[n-len[tmp]-1]!=c;tmp=fail[tmp]); 25 fail[now]=to[tmp][c]; 26 to[cur][c]=now; 27 } 28 last=to[cur][c]; 29 cnt[last]++; 30 p[n]=len[last]; 31 } 32 int main(){ 33 scanf("%s",s1+1); 34 int L=strlen(s1+1); 35 for(int i=1;i<=L;i++)s2[i]=s1[L-i+1]; 36 init(); 37 for(int i=1;i<=L;i++)add(s1[i]-'a',p1); 38 init(); 39 for(int i=1;i<=L;i++)add(s2[i]-'a',p2); 40 for(int i=1;i<=L-1;i++){ 41 ans=max(ans,p1[i]+p2[L-i]); 42 } 43 printf("%d",ans); 44 return 0; 45 }
