【HDOJ6599】I Love Palindrome String（PAM，manacher）

题意：给出一个由小写字母组成的长为n的字符串S，定义他的子串【L,R】为周驿东串当且仅当【L,R】为回文串且【L，(L+R)/2】为回文串

求i=【1，n】 所有长度为i的周驿东串的个数

n<=3e5

思路：PAM把所有回文串找出来，记录一下在原串S中的位置和长度，最后check每个结点是不是周驿东串

check部分可以正反两次哈希，N=3e5应该要双哈

我是用manacher预处理出每个位置为中心的回文串的最大长度，check就直接看一下(l+mid)/2这个位置的长度有没有r-mid+1

std里还判了右半部分，迷惑行为

#include<bits/stdc++.h>
using namespace std;
typedef long long ll;
typedef unsigned int uint;
typedef unsigned long long ull;
typedef pair<int,int> PII;
typedef pair<ll,ll> Pll;
typedef vector<int> VI;
typedef vector<PII> VII;
typedef pair<ll,ll>P;
#define N  700010
#define M  200010
#define fi first
#define se second
#define MP make_pair
#define pi acos(-1)
#define mem(a,b) memset(a,b,sizeof(a))
#define rep(i,a,b) for(int i=(int)a;i<=(int)b;i++)
#define per(i,a,b) for(int i=(int)a;i>=(int)b;i--)
#define lowbit(x) x&(-x)
#define Rand (rand()*(1<<16)+rand())
#define id(x) ((x)<=B?(x):m-n/(x)+1)
#define ls p<<1
#define rs p<<1|1

const ll MOD=1e9+7,inv2=(MOD+1)/2;
      double eps=1e-4;
      ll INF=1ll<<60;
      ll inf=5e13;
      int dx[4]={-1,1,0,0};
      int dy[4]={0,0,-1,1};

char s[N];
int ans[N],aid[N],p[N],a[N],ID;

int read()
{
   int v=0,f=1;
   char c=getchar();
   while(c<48||57<c) {if(c=='-') f=-1; c=getchar();}
   while(48<=c&&c<=57) v=(v<<3)+v+v+c-48,c=getchar();
   return v*f;
}

void manacher(char *ch,int len)
{
    rep(i,1,len*2+10) a[i]=-1;
    int n=2; a[1]=26; a[2]=27;
    rep(i,1,len)
    {
        a[++n]=ch[i]-'a'; aid[i]=n;
        a[++n]=27;
    }
    a[++n]=28;
    rep(i,0,n) p[i]=0;
    int mx=0,id=0;
    rep(i,2,n-1)
    {
        if(mx>i) p[i]=min(p[2*id-i],mx-i);
         else p[i]=1;
        while(a[i+p[i]]==a[i-p[i]]) p[i]++;
        if(p[i]+i>mx)
        {
            mx=p[i]+i;
            id=i;
        }
    }
    rep(i,2,n-1) p[i]--;
}

int isok(int l,int r)
{
    int mid=(l+r)>>1;
    int t=(p[(aid[l]+aid[mid])>>1]>=mid-l+1);
    return t;
}

struct pam
{
    int q,p,cnt[N],num[N],f[N],len[N],id[N],ch[N][26];

    pam(){}

    void init()
    {
        rep(i,0,ID)
        {
            len[i]=f[i]=num[i]=id[i]=cnt[i]=0;
            rep(j,0,25) ch[i][j]=0;
        }
        ID=1;
        f[0]=f[1]=1; len[1]=-1;
    }

    void add(int x,int n)
    {
        while(s[n-len[p]-1]!=s[n]) p=f[p];
        if(!ch[p][x])
        {
            int q=++ID,k=f[p];
            len[q]=len[p]+2;
            while(s[n-len[k]-1]!=s[n]) k=f[k];
            f[q]=ch[k][x];
            ch[p][x]=q;
            num[q]=num[f[q]]+1;
        }
        p=ch[p][x];
        cnt[p]++;
        id[p]=n;
    }

    void solve()
    {
        per(i,ID,0) cnt[f[i]]+=cnt[i];
        rep(i,0,ID)
         if(len[i]>0)
          ans[len[i]]+=isok(id[i]-len[i]+1,id[i])*cnt[i];
    }
}pam;



int main()
{
    while(scanf("%s",s+1)!=EOF)
    {
        int n=strlen(s+1);
        pam.init();
        rep(i,1,n) ans[i]=0;
        rep(i,1,n) pam.add(s[i]-'a',i);
        manacher(s,n);
        pam.solve();
        rep(i,1,n-1) printf("%d ",ans[i]);
        printf("%d
",ans[n]);
    }

    return 0;
}