Description
Input
输入分为两行,第一行为一个整数,表示字符串的长度,第二行有个连续的小写的英文字符,表示字符串的内容。
Output
输出文件只有一行,即:输入数据中字符串的最长双倍回文子串的长度,如果双倍回文子串不存在,则输出0。
首先由题意可知,只用考虑偶数长的回文串。这样就不用插入#,直接manacher即可。
其中p[i]表示以i和i+1的中点为对称轴,最大回文串长度的一半。
接下来考虑求解,设答案表示为i..j,k..x,x+1..y,y+1..q,这四段分别为w,wr,w,wr。只考虑x和y,则合法的条件为y-p[y]<=x【后两段是回文串】且y<=x+p[x]/2【整个是回文串】。
具体做法是先将y按y-p[y]排序,然后扫描x,把满足y-p[y]<=x的y插入set中,在查找最大的小于等于x+p[x]/2的y更新答案。
整体复杂度O(nlogn)。
#include<cstdio>
#include<cstring>
#include<algorithm>
#include<set>
using namespace std;
int l,p[500010],a[500010];
char s[500010];
set<int> se;
set<int>::iterator it;
bool cmp(const int &a,const int &b)
{return a-p[a]<b-p[b];
}
void mana()
{int i,j,k,id,mx=0;s[0]='@';s[l+1]='$';for (i=1;i<=l;i++){if (i<mx) p[i]=min(p[2*id-i],mx-i);else p[i]=0;while (s[i-p[i]]==s[i+p[i]+1]) p[i]++;if (i+p[i]>mx){mx=i+p[i];id=i;}}
}
void solve()
{int i,j,k,x,y,z,now=0,ans=0;for (i=1;i<=l;i++)a[i]=i;sort(a+1,a+l+1,cmp);se.insert(-1);se.insert(-2);for (i=1;i<=l;i++){while (now<l&&a[now+1]-p[a[now+1]]<=i){now++;se.insert(a[now]);}x=*(--se.upper_bound(i+p[i]/2));if (x<0) continue;ans=max(ans,4*(x-i));}printf("%d\n",ans);
}
int main()
{scanf("%d%s",&l,s+1);mana();solve();
}