UVA - 11107：Life Forms

后缀数组height+二分

#include<cstdio>
#include<cstdlib>
#include<algorithm>
#include<cstring>
#include<string>
#include<vector>
#include<iostream>
#define MAXN 100*10000+10
#define pii pair<int,int>
using namespace std;
struct Node{
    pii v;
    int p;
    friend bool operator < (const Node &p1,const Node &p2){
        return (p1.v<p2.v);
    }
}tsa[MAXN];

char s[MAXN];
int T,n;
int b[],p[MAXN],ans[MAXN];
int rk[MAXN],sa[MAXN];
void build_sa(){
    for(int i=;i<=n;i++){
        rk[i]=s[i]-;
    }
    for(int k=;k<=n;k<<=){
        for(int i=;i<=n-k;i++){
            tsa[i].v=make_pair(rk[i],rk[i+k]);
            tsa[i].p=i;
        }
        for(int i=n-k+;i<=n;i++){
            tsa[i].v=make_pair(rk[i],);
            tsa[i].p=i;
        }
        sort(tsa+,tsa+n+);
        int cnt=;
        for(int i=;i<=n;i++){
            if(!cnt||tsa[i].v!=tsa[i-].v){
                cnt++;
            }
            rk[tsa[i].p]=cnt;
        }
        if(cnt>=n)break;
    }
    for(int i=;i<=n;i++){
        sa[rk[i]]=i;
    }
}
int height[MAXN];
void build_height(){
    int k=;
    for(int i=;i<=n;i++){
        if(k)k--;
        int j=sa[rk[i]-];
        while(s[i+k]==s[j+k]){
            k++;
        }
        height[rk[i]]=k;
    }
}
void init(){
    n=;
    char tmp='z'+;
    for(int i=;i<=T;i++){
        scanf("%s",s+n+);
        int t=n;
        n=strlen(s+);
        s[++n]=(tmp++);
        for(int j=t+;j<=n;j++){
            p[j]=i;
        }
    }
    build_sa();
    build_height();
}
int check(int x){
    memset(b,,sizeof(b));
    int cnt=;
    int ret=;
    for(int i=;i<=n;i++){
        if(height[i]>=x){
            if(!cnt){
                if(p[sa[i]]==p[sa[i-]]){
                    cnt=;
                }
                else{
                    cnt=;
                }
                b[p[sa[i]]]=b[p[sa[i-]]]=;
            }
            else if(!b[p[sa[i]]]){
                b[p[sa[i]]]=;
                cnt++;
            }
            if(cnt>T/){
                ans[sa[i]]=;
                ret=;
            }
        }
        else{
            memset(b,,sizeof(b));
            cnt=;
        }
    }
    return ret;
}
vector<string> vs;
string ept;
void print(int x){
    vs.clear();
    string ts(s,n+);
    for(int i=;i<=n;i++){
        if(ans[i]){
            string t;
            t=ept;
            for(int j=i;j<=i+x-;j++){
                t+=s[j];
            }
            vs.push_back(t);
        }
    }
    sort(vs.begin(),vs.end());
    vs.erase(unique(vs.begin(),vs.end()),vs.end());
    for(int i=;i<vs.size();i++){
        cout<<vs[i]<<endl;
    }
}
void solve(){
//    for(int i=1;i<=n;i++){
//        for(int j=sa[i];j<=n;j++){
//            printf("%c",s[j]);
//        }
//        printf("\n");
//    }
//    for(int i=6;i<=n;i++){
//        printf("%d\n",check(i));
//    }
    int L=,R=n;
    while(R-L>){
        int mid=(L+R)/;
        if(check(mid)){
            L=mid;
        }
        else{
            R=mid;
        }
    }
    memset(ans,,sizeof(ans));
    if(R&&check(R)){
        print(R);
    }
    else{
        if(L&&check(L)){
            print(L);
        }
        else{
            printf("?\n");
        }
    }
}
int main()
{
//    freopen("data.in","r",stdin);
//    freopen("my.out","w",stdout);
    int flag=;
    while(){
        scanf("%d",&T);
        if(!T)break;
        if(flag){
            printf("\n");
        }
        else{
            flag=;
        }
        if(T==){
            scanf("%s",s+);
            printf("%s\n",s+);
            continue;
        }
        init();
        solve();
    }
    return ;
}