【回文自动机】ural2040

题意：依次添加字符询问增加的新的回文串的个数。

以前以为回文自动机和后缀自动机很像，都是增量维护，其实这东西和ac自动机更像一点。

回文自动机的本质也和ac自动机差不多，把每个回文子串的后半部分加到trie树中去（奇数就把中间也取进去），用fail指针表示当前节点代表的回文子串的最长回文后缀。

考虑现在已经处理好了串s的回文自动机，如何来处理s+‘c’，这里可以认识到，加了一个字符之后至多增加一个回文串，而且这个回文串就是s+'c'的最长回文后缀，假设为s'，如果有别的回文子串加入，一定也是这个s‘的一个回文后缀，那么用s'的对称轴对称一下就会发现这个回文子串一定在之前出现过了。

那么就用last表示当前s的最长回文子串在回文自动机中的节点，这里可以记一下last这个串的长度，然后就可以判断这个串添上’c'是不是还是回文子串，不是的话就沿fail指针找到满足的，找到满足的节点x之后，就可以通过next判断是否已经存在这个回文子串，如果存在就不用管，如果不存在就可以新添一个节点进去，同时维护一下新的点的fail指针，维护方法同ac自动机类似，不停找fail，直到满足条件。

参考培根的实现，对于奇数回文子串和偶数回文子串可以用一些技巧来维护，用两个节点root-1和root0来分别表示奇数的和偶数的起始，将root-1的儿子的fail指向root0，root0自己指向root-1，因为单个字符作为回文串的，加一个字符之后需要check是否是双字符回文串，而最终至少添一个字符之后其本身都是一个回文串因此最终要去root-1。

#include <cstdio>
#include <cstdlib>
#include <cstring>
#include <iostream>
#include <algorithm>
const int maxn=5005000;
using namespace std;
char ch[maxn],ans[maxn];
int ss,rood,reven,last;
int next[maxn][2],rt[maxn],l[maxn];
void origin()
{
    for (int i=0;i<=ss;i++) {
        for (int j=0;j<=1;j++) 
            next[i][j]=0;
        l[ss]=-2,rt[ss]=0;
    }
    ss=0;
    rood=++ss,l[rood]=-1,rt[rood]=rood;
    reven=++ss,l[reven]=0,rt[reven]=rood;
    last=rood;
}
bool palin(int x,int id,int chr)
{
    return ch[id-1-l[x]]-'a'==chr;
}
void add(int &last,int chr,int id)
{
    for (;!palin(last,id,chr);last=rt[last]) ;
    if (next[last][chr]) last=next[last][chr];
    else {
        int x=last;
        ++ss;
        next[x][chr]=ss,l[ss]=l[x]+2;
        if (x==rood) rt[ss]=reven;
        else {
            x=rt[x];
            for (;x!=rood && (!next[x][chr] || !palin(x,id,chr));x=rt[x]) ;
            if (next[x][chr] && palin(x,id,chr)) rt[ss]=next[x][chr];
            else rt[ss]=rood;
        }
        last=ss;
    }
}
int main()
{
    gets(ch+1);
    origin();
    int len=strlen(ch+1);
    for (int i=1;i<=len;i++) {
        int tmp=ss;
        add(last,ch[i]-'a',i);
        if (tmp!=ss) ans[i]='1';
        else ans[i]='0';
    }
    ans[len+1]='\0';
    puts(ans+1);
/*    for (int i=1;i<=ss;i++) {
        cout<<i<<endl;
        for (int j=0;j<=1;j++)
            cout<<(char)(j+'a')<<' '<<next[i][j]<<' ';
        cout<<endl;
    }
    for (int i=1;i<=ss;i++) cout<<rt[i]<<' ';cout<<endl;*/
    return 0;
}