【2019 Multi-University Training Contest 2】【HDU 6599 I Love Palindrome String】【回文自动机+马拉车】

题目链接：
http://acm.hdu.edu.cn/showproblem.php?pid=6586
题意：
这个题目要求求回文串的一半也是回文串的这种子串的数量（根据不同长度分类）。
题解：
如果只用manacher() 也是可以判断的但是会T （如果给出的字符有很多一样比如aaaa…aaaa这种）
这里先用回文树跑出本质不同的回文串（这样我们在cheak的时候会省去重复的判断）
然后再将原串manacher跑一遍（也可以字符串哈希）
然后用manacher跑出来的值判断一下回文串的一半子串是不是回文串。
代码：

#include<bits/stdc++.h>
using namespace std;
typedef  long long ll;
typedef  unsigned  long long ull;
const int N =3e5+50;
ll ret[N];
char s[N];
char new_s[N*2];
int pi[N*2];
bool Check(int l,int r)
{
    int len=r-l+1;
    if(pi[l+r+2]-1>=len)return true;
    return false;
}
bool check(int l,int r){
    int len=r-l+1;
    if(len<=2)
        return true;
    else
        return Check(l,(l+r)/2);
}
struct Palindromic_Tree{
    int nxt[N][30],fail[N],cnt[N];
    int num[N],len[N],s[N],id[N];
    int last,n,p;
    int newnode(int l){
        memset(nxt[p],0,sizeof(nxt[p]));
        cnt[p]=num[p]=0;
        len[p]=l;
        return p++;
    }
    
    void init(){
        p=0;
        newnode(0),newnode(-1);
        last=n=0;
        s[0]=-1;
        fail[0]=1;
    }
    
    int get_fail(int x){
        while(s[n-len[x]-1]!=s[n]) x=fail[x];
        return x;
    }
    
    void add(int c){
        c-='a';
        s[++n]=c;
        int cur=get_fail(last);
        if(!nxt[cur][c]){
            int now=newnode(len[cur]+2);
            fail[now]=nxt[get_fail(fail[cur])][c];
            nxt[cur][c]=now;
            num[now]=num[fail[now]]+1;
        }
        last=nxt[cur][c];
        cnt[last]++,id[last]=n;
    }
    void Count()
    {
        for(int i=p-1;i>=0;i--)cnt[fail[i]]+=cnt[i];
        for(int i=2;i<p;i++){
            if(check(id[i]-len[i],id[i]-1)){
                ret[len[i]]+=cnt[i];
            }
        }
    }
}pam;

int Init()
{
    int len=strlen(s);
    new_s[0]='$';
    new_s[1]='#';
    int j=2;
    for(int i=0;i<len;i++)
    {
        new_s[j++]=s[i];
        new_s[j++]='#';
    }
    new_s[j]='\0';
    return j;
}
void Manacher()
{
    int len=Init();
    int id;
    int mx=0;
    for(int i=1;i<len;i++)
    {
        if(i<mx)
            pi[i]=min(pi[2*id-i],mx-i);
        else
            pi[i]=1;
        while(new_s[i-pi[i]]==new_s[i+pi[i]])
            pi[i]++;
        if(mx<i+pi[i]) {
            id=i;
            mx=+i+pi[i];
        }
    }
}


int main()
{
    while(~scanf("%s",s)){
        memset(ret,0,sizeof(ret));
        pam.init();
        Init();
        Manacher();
        int len=strlen(s);
        for(int i=0;i<len;i++)
        {
            pam.add(s[i]);
        }
        pam.Count();
        printf("%lld",ret[1]);
        for(int i=2;i<=len;i++)
        {
            printf(" %lld",ret[i]);
        }
        printf("\n");
    }
    return 0;
}