后缀数组+LCP+多个串找公共子序列+uva11107

详见训练指南，需要注意的细节太多了

#include<iostream>
#include<cstdio>
#include<string>
#include<cstring>
#include<algorithm>
using namespace std;
const int MAXN=100010;
int n;
struct Suffix
{
    int s[MAXN],sa[MAXN];
    int height[MAXN],rank[MAXN];
    int c[MAXN],t[MAXN],t2[MAXN];
    int n;
    void clear(){n=0;memset(sa,0,sizeof(sa));}
    void build_sa(int m)
    {
        int i,*x=t,*y=t2;
        for(i=0; i<m; i++)c[i]=0;
        for(i=0; i<n; i++)c[x[i]=s[i]]++;
        for(i=1; i<m; i++)c[i]+=c[i-1];
        for(i=n-1; i>=0; i--)sa[--c[x[i]]]=i;
        for(int k=1; k<=n; k <<= 1)
        {
            int p=0;
            for(i=n-k; i<n; i++) y[p++]=i;
            for(i=0; i<n; i++) if(sa[i]>=k) y[p++]=sa[i]-k;
            for(i=0; i<m; i++) c[i]=0;
            for(i=0; i<n; i++) c[x[y[i]]]++;
            for(i=1; i<m; i++) c[i]+=c[i-1];
            for(i=n-1; i>=0; i--) sa[--c[x[y[i]]]]=y[i];
            swap(x,y);
            p=1;
            x[sa[0]]=0;
            for(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 i, j, k = 0;
        for(i = 0; i < n; i++) rank[sa[i]] = i;

        for(i = 0; i < n; i++)
        {
            if(k) k--;
            int j = sa[rank[i]-1];
            while(s[i+k] == s[j+k])
            {k++;}
            height[rank[i]] = k;
        }
    }
};
Suffix sa;
int idx[MAXN];
int flag[110];
bool good(int L,int R)
{
    memset(flag,0,sizeof(flag));
    if(R-L<=n/2)return false;
    int cnt=0;
    for(int i=L;i<R;i++)
    {
        if(idx[sa.sa[i]]!=n&&!flag[idx[sa.sa[i]]])
            flag[idx[sa.sa[i]]]=1,cnt++;
    }
    return cnt>n/2;
}
void print_sub(int L, int R) {
  for(int i = L; i < R; i++)
    printf("%c", sa.s[i] - 1 + 'a');
  printf("\n");
}

bool process(int len,bool print)
{
    int L=0;
    for(int R=1;R<=sa.n;R++)
    {
        if(R==sa.n||sa.height[R]<len)
        {
            if(good(L,R))
            {
                if(print)print_sub(sa.sa[L],sa.sa[L]+len);
                else return true;
            }
            L=R;
        }
    }
    return false;
}
void solve(int maxlen)
{
    if(!process(1,false))
        cout<<"?"<<endl;
    else
    {
        int L=1,R=maxlen,mid;
        while(L<R)
        {
            mid=L+(R-L+1)/2;
            if(process(mid,false))L=mid;
            else R=mid-1;
        }

        process(L,true);
    }
}
void add(int ch,int i)
{
    idx[sa.n]=i;
    sa.s[sa.n++]=ch;
}
int main()
{
    #ifndef ONLINE_JUDGE
        freopen("in.txt","r",stdin);
    #endif
    bool first=true;
    while(scanf("%d",&n)!=EOF,n)
    {
        if(first)first=false;
        else cout<<endl;
        char s[1010];
        int maxlen=0;
        sa.clear();
        for(int i=0;i<n;i++)
        {
            scanf("%s",s);
            int len=strlen(s);
            maxlen=max(maxlen,len);
            for(int j=0;j<len;j++)
                add(s[j]-'a'+1,i);
            add(100+i,n);
        }
        add(0,n);
        //cout<<sa.n<<endl;
        if(n==1)cout<<s<<endl;
        else
        {
            sa.build_sa(100+n);
            sa.getheight();
            solve(maxlen);
        }
    }
    return 0;
}

又写了一遍。。

#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/2)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/2;
}
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);}
                else return true;
            }
            l=i;
        }
    }
    return false;
}
void solve()
{
    int l=1,r=maxlen,mid;
    if(!can(1,false)){printf("?\n");return;}
    while(l<r)
    {
        mid=l+(r-l+1)/2;
        if(can(mid,false))l=mid;
        else r=mid-1;
    }
    can(l,true);
}
int main()
{
    //freopen("in.txt","r",stdin);
    bool first=true;
    while(scanf("%d",&t)!=EOF,t)
    {
        if(first)first=false;
        else printf("\n");
        n=maxlen=0;
        memset(sa,0,sizeof(sa));
        for(int i=0;i<t;i++)
        {
            scanf("%s",word);
            int len=strlen(word);
            maxlen=max(len,maxlen);
            for(int j=0;j<len;j++)
                add(word[j]-'a'+1,i);
            add(100+i,t);
        }
        if(t==1){printf("%s\n",word);continue;}
        add(0,t);
        build_sa(100+t);
        n--;
        getheight();
        solve();
    }
    return 0;
}