[POJ1226 Substrings]

[关键字]：后缀数组 字符串

[题目大意]：给定n 个字符串，求出现或反转后出现在每个字符串中的最长子串。

//=========================================================================

[分析]：和前几道题一样的做法，先见字符串和它反转后的串连在一起并用不同的字分割（小于原串的任意字符且不相同）最后加上一个最小的字符结束，再求后缀数组，二分答案对height数组分组判定是否有至少一组后缀包含所有字串。

[代码]：

View Code

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

const int MAXN=20010;
const int MAXL=120;

int test,x,n,df;
int top[MAXN],temp[MAXN],sa[MAXN],rank[MAXN],height[MAXN];
int s[MAXN],belong[MAXN];
char s1[MAXL],s2[MAXL];
bool v[MAXL];

void Make()
{
     int i,j,len,m;
     m=n<256?256:n;
     memset(top,0,sizeof(top));
     for (i=0;i<n;++i) ++top[rank[i]=s[i]&0xff];
     for (i=1;i<m;++i) top[i]+=top[i-1];
     for (i=0;i<n;++i) sa[--top[rank[i]]]=i;
     for (len=1;len<n;len<<=1)
     {
         for (i=0;i<n;++i)
         {
             j=sa[i]-len;
             if (j<0) j+=n;
             temp[top[rank[j]]++]=j;
         }
         sa[temp[top[0]=0]]=j=0;
         for (i=1;i<n;++i)
         {
             if (rank[temp[i]]!=rank[temp[i-1]] ||
                rank[temp[i]+len]!=rank[temp[i-1]+len])
                top[++j]=i;
             sa[temp[i]]=j;
         }
         memcpy(rank,sa,sizeof(sa));
         memcpy(sa,temp,sizeof(sa));
         if (j>=n-1) break;
     }
}
     
void Lcp()
{
     int i,j,k;
     for (j=rank[height[i=k=0]=0];i<n-1;++i,++k)
         while (k>=0 && s[i]!=s[sa[j-1]+k])
               height[j]=k--,j=rank[sa[j]+1];
} 

bool Judge(int limit)
{
     int i,d=0;
     memset(v,0,sizeof(v));
     v[0]=1;
     if (!v[belong[sa[0]]]) ++d;
     v[belong[sa[0]]]=1;
     for (i=1;i<n;++i)
     {
         if (height[i]<limit)
         {
                             d=0;
                             memset(v,0,sizeof(v));
                             v[0]=1;
         }
         if (!v[belong[sa[i]]]) ++d;
         if (d>=x) return 1;
         v[belong[sa[i]]]=1;
     }
     return 0;
}

void Solve()
{
     int mid,l=1,r=120,ans=0;
     while (l<=r)
     {
           mid=(l+r)/2;
           if (Judge(mid)) ans=mid,l=mid+1; else r=mid-1;
     }
     printf("%d\n",ans);
}

int main()
{
    freopen("in.txt","r",stdin);
    freopen("out.txt","w",stdout);
    scanf("%d",&test);
    while (test--)
    {
          df=n=0;
          scanf("%d",&x);
          for (int i=1;i<=x;++i)
          {
              scanf("%s",s1);
              //memset(s2,0,sizeof(s2));
              int len=strlen(s1);
              for (int j=0;j<len;++j)
                  s[n++]=(int)s1[j]+200,s2[len-j-1]=s1[j],belong[n-1]=i;
              s[n++]=++df;
              for (int j=0;j<len;++j)
                  s[n++]=(int)s2[j]+200,belong[n-1]=i;
              s[n++]=++df;
              //printf("%s %s\n",s1,s2);
          }
          s[n-1]=0;
          if (x==1) {printf("%d\n",strlen(s1));continue;}
          Make();
          Lcp();
          //for (int i=0;i<n;++i) printf("%d ",s[i]);
          //printf("\n");
          Solve();
    }
    return 0;
}