BZOJ2084: [Poi2010]Antisymmetry

$n leq 500000$的01串，1跟0配，问最长回文子串。

$0=1$，$1 eq 1$，$0 eq 0$，然后二分哈希或manacher或回文树。

//#include<iostream>
#include<cstring>
#include<cstdio>
//#include<time.h>
//#include<complex>
//#include<set>
//#include<queue>
#include<algorithm>
#include<stdlib.h>
using namespace std;

#define LL long long
int qread()
{
    char c; int s=0,f=1; while ((c=getchar())<'0' || c>'9') (c=='-') && (f=-1);
    do s=s*10+c-'0'; while ((c=getchar())>='0' && c<='9'); return s*f;
}

//Pay attention to '-' , LL and double of qread!!!!

int n;
#define maxn 500011
char s[maxn];
struct PT
{
    struct Node{int ch[2],pre,len,cnt;}a[maxn];
    int last,size;
    PT() {last=size=1; a[1].pre=0; a[1].len=0;}
    void insert(int c,int p)
    {
//        cout<<p<<' '<<last<<"!!!!!
";
        int y=last,len;
        while (y)
        {
//            cout<<y<<' '<<a[y].len<<' '<<a[y].pre<<endl;
            len=a[y].len;
            if (p-len-1>=1 && s[p-len-1]!=s[p]) break;
            y=a[y].pre;
        }
        if (!y) {last=1; return;}
//        cout<<y<<' '<<a[y].len<<' '<<a[y].pre<<endl;
        if (a[y].ch[c]) {last=a[y].ch[c]; a[last].cnt++; return;}
        int x=last=++size;
        a[x].len=a[y].len+2; a[y].ch[c]=x; a[x].cnt=1;
        while (y)
        {
            y=a[y].pre; len=a[y].len;
            if (p-len-1>=1 && s[p-len-1]!=s[p])
            {
                a[x].pre=a[y].ch[c];
                break;
            }
        }
        if (!a[x].pre) a[x].pre=1;
    }
//    void test()
//    {
//        for (int i=1;i<=size;i++) cout<<a[i].len<<' '<<a[i].pre<<' '<<a[i].cnt<<endl;
//    }
    void hei()
    {
        for (int i=size;i;i--) a[a[i].pre].cnt+=a[i].cnt;
    }
}pt;

int main()
{
    scanf("%d",&n); scanf("%s",s+1);
    for (int i=1;i<=n;i++) pt.insert(s[i]-'0',i);
//    pt.test();
    pt.hei();
    int ans=0; for (int i=2;i<=pt.size;i++) ans+=pt.a[i].cnt;
    printf("%d
",ans);
    return 0;
}