hdu多校2 I I Love Palindrome String 回文自动机+hash

回文自动机学习链接：https://blog.youkuaiyun.com/lwfcgz/article/details/48739051  

                                    https://www.cnblogs.com/Xu-daxia/p/10544128.html#autoid-0-3-0

对于此题，回文自动机记录本质不同回文串出现次数，hash判断这个回文串的前一部分是否是回文串（因为这个串本身就是回文串，如果这个串前半部分的hash值等于后半部分的hash值，说明这个串前半部分跟后半部分相等，就说明这个串前半部分是个回文串）。

#include<bits/stdc++.h>
using namespace std;
const int maxn=3e5+10;
char s[maxn];
int hs[maxn],pp[maxn];
long long ans[maxn];
const int tt=7;
int geths(int l,int r)
{
    return hs[r]-hs[l-1]*pp[r-l+1];
}
int check(int l,int r)
{
    int len=r-l+1;
    int mid=r+l>>1;
    if(len&1)
        return geths(l,mid)==geths(mid,r);
    else
        return geths(l,mid)==geths(mid+1,r);
}
struct PAM
{
    int len[maxn],fail[maxn],nxt[maxn][26],id[maxn],num[maxn],tot,last,L[maxn];
    long long siz[maxn];
    void init()
    {
        len[0]=0;
        fail[0]=1;
        len[1]=-1;
        fail[1]=0;
        tot=1;
        last=0;
        memset(nxt[1],0,sizeof(nxt[1]));
        memset(nxt[0],0,sizeof(nxt[0]));
        memset(siz,0,sizeof(siz));
    }
    int new_node(int x)
    {
        int now=++tot;
        memset(nxt[tot],0,sizeof(nxt[tot]));
        len[now]=x;
        return now;
    }
    void ins(int c,int n)
    {
        int u=last;
        while(s[n-len[u]-1]!=s[n])u=fail[u];
        if(nxt[u][c]==0)
        {
            int now=new_node(len[u]+2);
            int v=fail[u];
            while(s[n-len[v]-1]!=s[n])v=fail[v];
            fail[now]=nxt[v][c];
            nxt[u][c]=now;
            num[now]=num[fail[now]]+1;
        }
        last=nxt[u][c];
        siz[last]++;
        L[n]=len[last];
        id[last]=n;
    }
    void solve()
    {
        for(int i=tot; i>=0; i--)
        {
            siz[fail[i]]+=siz[i];//每一个后缀出现的次数
        }
        for(int i=2; i<=tot; i++)
        {
            //cout<<i<<" "<<id[i]<<" "<<len[i]<<endl;
            if(check(id[i]-len[i]+1,id[i]))
            {
                ans[len[i]]+=siz[i];
            }
        }
        int length=strlen(s+1);
        for(int i=1; i<length; i++)
        {
            printf("%lld ",ans[i]);
        }
        printf("%lld\n",ans[length]);
    }
} pam;
int main()
{
    pp[0]=1;
    for(int i=1; i<maxn-5; i++)
        pp[i]=pp[i-1]*tt;
    while(~scanf("%s",s+1))
    {
        int len= strlen(s+1);
        memset(ans,0,sizeof(ans));
        hs[0]=0;
        for(int i=1; i<=len; i++)
            hs[i]=hs[i-1]*tt+s[i];
        pam.init();
        for(int i=1; i<=len; i++)
            pam.ins(s[i]-'a',i);
        pam.solve();
    }
    return 0;
}