【poj3294-不小于k个字符串中最长公共子串】后缀数组

1.注意每两个串之间的连接符要不一样。

2.分组的时候要注意最后一组啊！又漏了！

3.开数组要考虑连接符的数量。100010是不够的至少要101000。

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

 const int N=;
 int n,cl,sl,ans,tt,c[N],tl[N],tr[N],al[N],ar[N],rk[N],Rs[N],sa[N],wr[N],y[N],h[N],st[N],ed[N];
 char s[];
 bool vis[];

 int minn(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[j+k]==c[i+k] && j+k<=cl && i+k<=cl) k++;
         h[rk[i]]=k;
     }
     h[]=;
 }
 int idx(int x)
 {
     for(int i=;i<=n;i++)
         if(st[i]<=x && x<=ed[i]) return i;
     return ;
 }

 bool check(int k)
 {
     memset(vis,,sizeof(vis));
     int now=cl;
     tt=;
     for(int i=;i<=cl;i++)
     {
         if(h[i]<k)
         {
             int cnt=;
             for(int j=;j<=n;j++) if(vis[j]) cnt++;
             if(cnt>(n/)) tl[++tt]=sa[i-],tr[tt]=tl[tt]+now-;
             memset(vis,,sizeof(vis));
             now=cl-sa[i]+;vis[idx(sa[i])]=;
         }
         else
         {
             now=minn(now,h[i]);
             vis[idx(sa[i])]=;
         }
     }
     int cnt=;
     for(int j=;j<=n;j++) if(vis[j]) cnt++;
     if(cnt>(n/)) tl[++tt]=sa[cl-],tr[tt]=tl[tt]+now-;
     if(tt) return ;
     return ;
 }

 int main()
 {
     freopen("a.in","r",stdin);
     int T=;
     while()
     {
         T++;
         scanf("%d",&n);
         if(n==) return ;
         cl=;ans=;
         for(int i=;i<=n;i++)
         {
             scanf("%s",s+);
             sl=strlen(s+);
             if(i>) c[++cl]=i;
             st[i]=cl+;
             for(int j=;j<=sl;j++) c[++cl]=s[j];
             ed[i]=cl;
         }
         if(n==) {printf("%c\n",c[]);continue;}
         get_sa();
         get_h();
         // for(int i=1;i<=cl;i++) printf("%c",c[i]);printf("\n");
         // for(int i=1;i<=cl;i++) printf("%d ",sa[i]);printf("\n");
         // for(int i=1;i<=cl;i++) printf("%d ",rk[i]);printf("\n");
         // for(int i=1;i<=cl;i++) printf("%d ",h[i]);printf("\n");
         int l=,r=cl,mid;
         while(l<r)
         {
             mid=(l+r+)/;
             if(check(mid))
             {
                 l=mid;
                 ans=tt;
                 for(int i=;i<=tt;i++) al[i]=tl[i],ar[i]=tr[i];
             }
             else r=mid-;
         }
         if(T>) printf("\n");
         if(ans)
         {
             for(int i=;i<=ans;i++)
             {
                 for(int j=al[i];j<=ar[i];j++) printf("%c",c[j]);
                 printf("\n");
             }
         }
         else printf("?\n");
     }
     return ;
 }