bzoj2342 双倍回文【manacher+set】

解题思路：

开始以为只用r[i]是4的倍数的位置去检查其两侧是否也是回文就行了，但后面发现是错的，因为所求回文有可能在一个大回文内部，如aaabbaabbaaa，所以中间的r[i]会比更大;

后面看了下题解，还是太天真了。
枚举对称轴x，考虑用len(x,y)*4更新答案，则必须满足y-r[y]<=x，其中y<=x+r[x]/2。
那么我们按照y-r[y]排序后，依次把y插入set，每次询问比x+p[x]/2小的最大值就可以了。
由于我是加了’#’号的，写起来有些不同，不过大体思路是一样的。

#include<iostream>#include<cstdio>#include<cstring>#include<string>#include<algorithm>#include<cmath>#include<ctime>#include<queue>#include<vector>#include<set>using namespace std;const int N=500005;int n,m,ans;int r[N<<1];char t[N],s[N<<1];struct node{    int x,id;    inline friend bool operator <(const node &a,const node &b)    {        return a.x<b.x;    }}a[N<<1];set<int>S;void manacher(){    int mx=0,po=0;    for(int i=1;i<=n;i++)    {        if(mx>i)r[i]=min(r[po*2-i],mx-i);        else r[i]=1;        while(s[i+r[i]]==s[i-r[i]])r[i]++;        if(i+r[i]>mx)po=i,mx=i+r[i];    }}int main(){    //freopen("lx.in","r",stdin);    //freopen("lx.out","w",stdout);    scanf("%d%s",&m,t+1);    s[0]='!';    for(int i=1;i<=m;i++)        s[++n]='#',s[++n]=t[i];    s[++n]='#',s[++n]='?';    manacher();    m=0;    for(int i=1;i<=n;i+=2)a[++m].x=i-r[i],a[m].id=i;    sort(a+1,a+m+1);    int j=1;    for(int i=1;i<=n;i+=2)    {        while(j<=m&&a[j].x<i)S.insert(a[j++].id);        ans=max(ans,(*--S.lower_bound(i+(r[i]+1)/2)-i)*2);    }    cout<<ans;    return 0;}