【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+,inv2=(MOD+)/;
       double eps=1e-;
       ll INF=1ll<<;
       ll inf=5e13;
       int dx[]={-,,,};
       int dy[]={,,-,};

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

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

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

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

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

     pam(){}

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

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

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

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

     return ;
 }