GDKOI2016 Day2 T3 项链

T3 项链

给出一个环状字符串，求删除连续一段后剩下的字符为对称的最大长度、

首先你要知道怎样的一个字符串算是对称的。一定是由两个回文串拼成的。我们可以倍长原串，然后枚举开头，这样删除的一定是枚举的串的前缀。
用manacher预处理出所有的回文中心的扩展半径。两个合法的回文中心i和j一定满足j-i<=n&j-p[j]<=i+p[i]。并且它的串的长度就是j-i。于是我们用线段树维护一个区间最左边的i+p[i]比j-p[j]大的位置即可。

#include<cstdio>
#include<cstring>
#include<algorithm>
#define fo(i,a,b) for(int i=a;i<=b;i++)
#define N 200005
#define inf 0x7fffffff
using namespace std;
int p[2*N],t[N*6],n,tot,id,mx,ans,k;
char s[N],ch[2*N];
void change(int v,int l,int r,int x,int y) {
    if (l==r) {t[v]=y;return;}
    int m=(l+r)/2;
    if (x<=m) change(v*2,l,m,x,y);else change(v*2+1,m+1,r,x,y);
    t[v]=max(t[v*2],t[v*2+1]);
}
void get(int v,int l,int r,int x) {
    if (l==r) {if (t[v]>=x) k=l;return;}
    int m=(l+r)/2;
    if (t[v*2]>=x) get(v*2,l,m,x);
    else if (t[v*2+1]>=x) get(v*2+1,m+1,r,x);
}
void find(int v,int l,int r,int x,int y,int z) {
    if (k||x>y) return;
    if (l==x&&r==y) {
        get(v,l,r,z);return;
    }
    int m=(l+r)/2;
    if (y<=m) find(v*2,l,m,x,y,z);
    else if (x>m) find(v*2+1,m+1,r,x,y,z);
    else find(v*2,l,m,x,m,z),find(v*2+1,m+1,r,m+1,y,z);
}
void ext(int x) {
    while (p[x]+x<tot&&x-p[x]>1&&ch[p[x]+x+1]==ch[x-p[x]-1]) p[x]++;
}
int main() {
    scanf("%s",s+1);n=strlen(s+1);
    fo(i,n+1,n*2) s[i]=s[i-n];ch[++tot]='#';
    fo(i,1,n*2) ch[++tot]=s[i],ch[++tot]='#';
    fo(i,1,tot*3) t[i]=inf;
    fo(i,1,tot) {
        if (i>=mx) ext(i);else {
            int j=id*2-i;
            if (p[j]<mx-i) p[i]=p[j];else 
            p[i]=mx-i,ext(i);
        }
        if (p[i]+i>mx) mx=p[i]+i,id=i;
        k=0;find(1,1,tot,max(1,i-n),i-1,i-p[i]);
        if (k) ans=max(ans,i-k);
        change(1,1,tot,i,i+p[i]);
    }
    printf("%d",ans);
}