【hdu3948-不同回文串的个数】后缀数组

题意：求不同回文串的个数 n<=10^5

题解：

先按照manacher的构造方法改造一遍串，然后跑一遍manacher。

如ababa--> $#a#b#a#b#a#@

然后跑一遍后缀数组。

对于一个后缀sa[i]~cl（cl为字符串的总长），我们本来是要加上以sa[i]为中心的回文串的个数p[sa[i]]。

但是这可能有重复！

我们可以维护一个tmp，也就是上图中蓝色的框。tmp表示以字符sa[i-1]为中心已经被统计过的回文串的个数。

到了当前的sa[i]，tmp=min(tmp,h[i]);

每次如果p[x]<=tmp，就continue；

否则，ans+=（p[x]-tmp）/2;（/2是因为有#）

按照最上面的构造方法，每个最长回文串p[i]必定以#开头和结尾，所以可以直接除以2，可以画个图看看。

一开始以为z的ascll码<=100,get_sa()那里参数小了。。TLE了好几遍才发现。。

 #include<cstdio>
 #include<cstdlib>
 #include<cstring>
 #include<iostream>
 using namespace std;

 const int N=;
 int cl,sl,p[N],rk[N],Rs[N],sa[N],wr[N],y[N],h[N];
 char c[N],s[N];

 int minn(int x,int y){return x<y ? x:y;}
 int maxx(int x,int y){return x>y ? x:y;}

 void get_sa(int m)
 {
     for(int i=;i<=cl;i++) rk[i]=c[i];
     for(int i=;i<=m;i++) Rs[i]=;
     for(int i=;i<=cl;i++) Rs[rk[i]]++;
     for(int i=;i<=m;i++) Rs[i]+=Rs[i-];
     for(int i=cl;i>=;i--) sa[Rs[rk[i]]--]=i;

     int ln=,p=;
     while(p<cl)
     {
         int k=;
         for(int i=cl-ln+;i<=cl;i++) y[++k]=i;
         for(int i=;i<=cl;i++) if(sa[i]>ln) y[++k]=sa[i]-ln;

         for(int i=;i<=cl;i++) wr[i]=rk[y[i]];
         for(int i=;i<=m;i++) Rs[i]=;
         for(int i=;i<=cl;i++) Rs[wr[i]]++;
         for(int i=;i<=m;i++) Rs[i]+=Rs[i-];
         for(int i=cl;i>=;i--) sa[Rs[wr[i]]--]=y[i];

         for(int i=;i<=cl;i++) wr[i]=rk[i];
         for(int i=cl+;i<=cl+ln;i++) wr[i]=;
         p=;rk[sa[]]=;
         for(int i=;i<=cl;i++)
         {
             if(wr[sa[i]]!=wr[sa[i-]] || wr[sa[i]+ln]!=wr[sa[i-]+ln]) p++;
             rk[sa[i]]=p;
         }
         ln*=,m=p;
     }
     sa[]=;rk[]=;
 }

 void get_h()
 {
     int k=,j;
     for(int i=;i<=cl;i++) if(rk[i]!=)
     {
         j=sa[rk[i]-];
         if(k) k--;
         while(c[i+k]==c[j+k] && i+k<=cl && j+k<=cl) k++;
         h[rk[i]]=k;
     }
     h[]=;
 }

 void manacher()
 {
     int id=,mx=;
     p[]=;
     for(int i=;i<=cl;i++)
     {
         if(i+p[*id-i]- < mx) p[i]=p[*id-i];
         else
         {
             p[i]=maxx(,mx-i+);
             while(c[i+p[i]]==c[i-p[i]] && i+p[i]<=cl && i-p[i]>=) p[i]++;
             if(i+p[i]->mx) mx=i+p[i]-,id=i;
         }
     }
 }

 void solve()
 {
     int x,tmp=,ans=;
     for(int i=;i<=cl;i++)
     {
         x=sa[i];
         tmp=minn(tmp,h[i]);
         if(p[x]<=tmp) continue;
         ans+=(p[x]-tmp)/;
         tmp=p[x];
     }
     printf("%d\n",ans);
 }

 int main()
 {
     freopen("a.in","r",stdin);
     // freopen("me.out","w",stdout);
     int T;
     scanf("%d",&T);
     for(int TT=;TT<=T;TT++)
     {
         scanf("%s",s+);
         sl=strlen(s+);
         cl=;
         c[++cl]='$';
         int i;
         for(i=;i<=sl;i++)
             c[++cl]='#',c[++cl]=s[i];
         c[++cl]='#';c[++cl]='@';
         // for(int i=1;i<=cl;i++) printf("%c",c[i]);printf("\n");
         // for(int i=1;i<=cl;i++) printf("%d ",p[i]);printf("\n");
         manacher();
         get_sa();
         get_h();
         printf("Case #%d: ",TT);
         solve();
     }
     return ;
 }