后缀数组（多个穿的最长公共子串）+KMP+poj3080

跟Life Forms很像

#include<iostream>
#include<cstdio>
#include<cstring>
#include<vector>
#include<cmath>
#include<queue>
#include<stack>
#include<map>
#include<string>
#include<set>
#include<algorithm>
using namespace std;
const int maxn=1010*100+10;
int sa[maxn],height[maxn],rank[maxn],t1[maxn],t2[maxn],c[maxn],str[maxn];
int vis[110],idx[maxn];
char word[1010];
int n,t,maxlen;
void add(int x,int id)
{
    idx[n]=id;
    str[n++]=x;
}
void build_sa(int m)
{
    int *x=t1,*y=t2;
    for(int i=0;i<m;i++)c[i]=0;
    for(int i=0;i<n;i++)c[x[i]=str[i]]++;
    for(int i=1;i<m;i++)c[i]+=c[i-1];
    for(int i=n-1;i>=0;i--)sa[--c[x[i]]]=i;
    for(int k=1;k<=n;k<<=1)
    {
        int p=0;
        for(int i=n-k;i<n;i++)y[p++]=i;
        for(int i=0;i<n;i++)if(sa[i]>=k)y[p++]=sa[i]-k;
        for(int i=0;i<m;i++)c[i]=0;
        for(int i=0;i<n;i++)c[x[y[i]]]++;
        for(int i=1;i<m;i++)c[i]+=c[i-1];
        for(int i=n-1;i>=0;i--)sa[--c[x[y[i]]]]=y[i];
        swap(x,y);p=1;
        x[sa[0]]=0;
        for(int i=1;i<n;i++)
            x[sa[i]]=(y[sa[i-1]]==y[sa[i]]&&y[sa[i-1]+k]==y[sa[i]+k])?p-1:p++;
        if(p>=n)break;
        m=p;
    }
}
void getheight()
{
    int k=0;
    for(int i=1;i<=n;i++)rank[sa[i]]=i;
    for(int i=0;i<n;i++)
    {
        if(k)k--;
        int j=sa[rank[i]-1];
        while(str[i+k]==str[j+k])k++;
        height[rank[i]]=k;
    }
}
bool judge(int l,int r)
{
    if(r-l<t)return false;
    memset(vis,0,sizeof(vis));
    int cnt=0;
    for(int i=l;i<r;i++)
    {
        int tmp=sa[i];
        if(idx[tmp]!=t&&!vis[idx[tmp]]){cnt++,vis[idx[tmp]]=1;}
    }
    return cnt>=t;
}
void print(int l,int r)
{
    for(int i=l;i<r;i++)printf("%c",str[i]-1+'A');
    printf("\n");
}
bool can(int len,bool flag)
{
    int l=0;
    for(int i=1;i<=n;i++)
    {
        if(i==n||height[i]<len)
        {
            if(judge(l,i))
            {
                if(flag){print(sa[l],sa[l]+len);return true;}
                else return true;
            }
            l=i;
        }
    }
    return false;
}
void solve()
{
    int l=0,r=60,mid;
    while(l<r)
    {
        mid=l+(r-l+1)/2;
        if(can(mid,false))l=mid;
        else r=mid-1;
    }
    if(l>=3)can(l,true);
    else printf("no significant commonalities\n");
}
int main()
{
    freopen("in.txt","r",stdin);
    int T;
    scanf("%d",&T);
    while(T--)
    {
        scanf("%d",&t);
        n=maxlen=0;
        memset(sa,0,sizeof(sa));
        for(int i=0;i<t;i++)
        {
            scanf("%s",word);
            int len=strlen(word);
            for(int j=0;j<len;j++)
                add(word[j]-'A'+1,i);
            add(100+i,t);
        }
        add(0,t);
        build_sa(100+t);
        n--;
        getheight();
        solve();
    }
    return 0;
}

也可以用KMP做，枚举第一个串的每一个字串，然后多后面的串KMP，看存不存在这个字串

#include<cstdio>
#include<cstring>
char str[20][70];
char tmp[70];
int p[70];
void getp(int m,char *b1){
    char b[70];
    strcpy(b+1,b1);
    p[1]=0;
    int i,j=0;
    for(i=2;i<=m;i++){
        while(j>0&&b[j+1]!=b[i]) j=p[j];
        if(b[j+1]==b[i]) j+=1;
        p[i]=j;
    }
}
bool kmp(char *a1,char *b1,int n,int m)
{
    char a[70],b[70];
    strcpy(a+1,a1);strcpy(b+1,b1);
    int i,j=0,cnt=0;
    for(i=1;i<=n;i++){
        while(j>0&&b[j+1]!=a[i]) j=p[j];
        if(b[j+1]==a[i]) j+=1;
        if(j==m){
             return true;
        }
    }
    return false;
}
bool check(char *s,int tot)
{
     int i,j;
     for(i=2;i<=tot;i++)
     {
         int n=strlen(str[i]+1),m=strlen(s);
         if(!kmp(str[i]+1,s,n,m))
             return false;
     }
     return true;
}
char ans[70];
int main()
{
    int t,n,i,j;
    scanf("%d",&t);
    while(t--)
    {
        scanf("%d",&n);
        for(i=1;i<=n;i++)   scanf("%s",str[i]+1);
        int len=strlen(str[1]+1);
        int L=0;
        for(i=1;i<=len;i++)
        {
            for(j=1;j<=len-i+1;j++)
            {
                strncpy(tmp,str[1]+j,i);
                memset(p,0,sizeof(p));
                getp(i,tmp);
                if(check(tmp,n))
                {            
                    if(i>=L)
                    {
                    if(strcmp(tmp,ans)<0&&i==L) strcpy(ans,tmp);
                    if(i>L)   strcpy(ans,tmp);
                       L=strlen(ans);
                    }
                }
                memset(tmp,0,sizeof(tmp));
            }
        }
        if(L>=3) puts(ans);
        else puts("no significant commonalities");
    }
}