【BZOJ】2342: [Shoi2011]双倍回文-manacher&set

版权声明：欢迎转载(请附带原链接)ヾ(๑╹◡╹)ﾉ" https://blog.csdn.net/corsica6/article/details/84930970

传送门:bzoj2342

---

题解

回文自动机很好做，但这里讲一下manacher的做法。

设偶回文中心 $i$的回文半径为 $R_i$。

枚举回文中心 $i$，找到最小的 $j$满足: $j\geq i-\dfrac{R_i}{2}$且 $R_j+j\geq i$

将回文中心按 $R_j+j$升序排序并加入set， $1-n$枚举 $i$时 $erase$所有 $R_j+j<i$，二分找出现存的最小的 $j$即可。

---

代码

#include<bits/stdc++.h>
using namespace std;
const int N=1e6+100;

int n,m,f[N],rk[N],nt[N],ans;
char s[N],t[N];
set<int>st;

inline void manacher()
{
    int i,j,ct=0,mx=0;
    for(i=1;i<=n;++i){
        for(j=mx>i?min(mx-i,f[(ct<<1)-i]):1;t[i-j]==t[i+j];++j);
        f[i]=j;if(i+j>mx) mx=i+j,ct=i;
    }m=0;
    for(i=3;i<=n;i+=2) nt[++m]=f[i]>>1;
    m--;
}

inline bool cmp(const int&x,const int&y){return x+nt[x]<y+nt[y];}

int main(){
    int i,j,x;
    scanf("%d%s",&m,s+1);
    t[0]='!';
    for(i=1;i<=m;++i) t[++n]='#',t[++n]=s[i];
    t[++n]='#';t[++n]='?';
    manacher();for(i=1;i<=m;++i) rk[i]=i,st.insert(i);
    sort(rk+1,rk+m+1,cmp);
    set<int>::iterator it;
    for(i=j=1;i<=m;++i){
        for(;j<=m && rk[j]+nt[rk[j]]<i;++j) st.erase(rk[j]);
        x=*st.lower_bound(i-(nt[i]>>1));ans=max(ans,(i-x)<<2);
    }
    printf("%d\n",ans);
    return 0;
}