uva 11107Life Forms

题目：给你n(<100)个长度不超过1000的字符串，问至少是n/2个字符串的子串的字符串的最大长度，有多组按字典序输出所有的字符串。

分析：首先可以二分答案，然后判断长度为l的字符串是不是至少是n/2个字符串的子串

如何判断？

将n个字符串拼接在一起，中间用不相同的特殊字符隔开，然后再判断长度为l的时候，遍历一遍高度数组，要求连续>=l，

统计每个原字符串是不是出现过，我这边用了一个set去判断当前已经出现过多少种。当高度数组<l的时候，就判断一下当前

set的大小是不是超过n/2，然后从新开始统计。

#include <iostream>
#include <cstdio>
#include <cstring>
#include <algorithm>
#include <queue>
#include <string>
#include <vector>
#include <cmath>
#include <set>
#define maxn 200*1000+10
#define maxm 2010
#define maxt 110
#define INF 0x3f3f3f3f
#define mod 1009
#define MAX_STATE 2500
using namespace std;
int n,k;
int N;
int Rank[maxn];
int temp[maxn];
int LCP[maxn];
int SA[maxn];
string T;
//比较(rank[k,i],rank[k,i+k])与(rank[k,j],rank[k,j+k])
bool cmp_sa(int i,int j){
    if(Rank[i]!=Rank[j])return Rank[i]<Rank[j];
    else{
        int ri = i+k<=n?Rank[i+k]:-1;
        int rj = j+k<=n?Rank[j+k]:-1;
        return ri<rj;
    }
}
//计算后缀数组O(nlognlogn)
void construct_sa(string S,int *sa){
    n = S.length();
    //长度为1的字符串,Rank直接取字符的编码
    for(int i=0;i<=n;++i){
        sa[i]=i;
        Rank[i]=i<n?S[i]:-1;
    }
    //利用对长度为k的字符串排序计算长度位2k的顺序
    for(k=1;k<=n;k*=2){//注意这里不是 int k
        sort(sa,sa+n+1,cmp_sa);
        //计算新的rank，暂存到temp中
        temp[sa[0]]=0;
        //调整rank,相同的字符串的rank时一样的
        for(int i=1;i<=n;++i){
            temp[sa[i]] = temp[sa[i-1]]+ (cmp_sa(sa[i-1],sa[i])?1:0);
        }
        //存回rank
        for(int i=0;i<=n;++i){
            Rank[i]=temp[i];
        }
    }
}
//计算LCP，O（n）
void construct_lcp(string S,int *sa,int *lcp){
    n = S.length();
    for(int i=0;i<=n;++i)Rank[sa[i]]=i;
    int h=0;
    lcp[0]=0;
    for(int i=0;i<n;++i){
        //计算字符串从i位置开始的后缀及其在后缀数组前一个的后缀的LCP
        int j = sa[Rank[i]-1];
        //将h先减去首字母的1的长度，保持前缀相同前提下不断增加
        if(h>0)h--;
        for(;j+h<n&&i+h<n;++h){
            if(S[j+h]!=S[i+h])break;
        }
        lcp[Rank[i]-1]=h;
    }
}

int ID[maxn];
bool judge(int len,int flag){
    set<int>used;
    used.insert(ID[SA[0]]);
    for(int i=1;i<=n;++i){
        while(i<n&&LCP[i-1]>=len){
            used.insert(ID[SA[i]]);
            i++;
        }
        if(used.size()*2>N){
            if(flag)return true;
            else{
                for(int j=0;j<len;++j){
                    cout<<T[SA[i-1]+j];
                }
                cout<<endl;
            }
        }
        used.clear();
        used.insert(ID[SA[i]]);
    }
    return false;
}
void solve(int L,int R){
    if(!judge(1,1)){
        cout<<"?"<<endl;
        return;
    }
    while(L<R){
        int mid = (L+R)/2;
        if(judge(mid,1))L = mid+1;
        else R = mid;
    }
    judge(L-1,0);
}
int main(){
    //freopen("in.txt","r",stdin);
    //freopen("out.txt","w",stdout);
    ios::sync_with_stdio(false);
    int CASE=0;
    while(cin>>N){
        if(N==0)break;
        if(CASE!=0)cout<<endl;
        CASE++;
        string S;
        T="";
        if(N==1){
            cin>>S;
            cout<<S<<endl;
            continue;
        }
        int mmax =0;
        int cnt=0;
        for(int i=0;i<N;++i){
            cin>>S;
            int sl = S.length();
            mmax = max(mmax,sl);
            for(int j=0;j<sl;++j){
                ID[cnt]=i;
                cnt++;
            }
            T+=S;
            T+='z'+i+1;
            cnt++;
        }
        construct_sa(T,SA);
        construct_lcp(T,SA,LCP);
        solve(1,mmax+1);
    }
}