UVA 11107 生命的形式（後綴數組+LCP）

把所有輸入的字符串拼起來，二分答案，每次判斷是否有一個長度爲p的串在超過一半的串中連續出現，判斷方法是掃描height數組，因爲height數組中，相同串長度都聚集在一起。

 

#include<bits/stdc++.h>
using namespace std;
const int maxn=100005;
char s[maxn];
int sa[maxn],t[maxn],t2[maxn],c[maxn];
int belong[maxn];
int n;
vector<int> pos[1005];
void build(char *s,int n,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=0;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=0;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;
	}
	
}

int rank[maxn],height[maxn];

void getheight(int n){
	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;
	}
}



bool check(int mid,int N){
	bool vis[110];
	memset(vis,0,sizeof(vis));
	int cnt=0;
	bool ret = false, in = false;

	for(int i=0;i<=N;i++){
		if(height[i]>=mid){
			if(!vis[belong[sa[i]]]){
				cnt++;
				vis[belong[sa[i]]]=1;
			}
			if(cnt>(n>>1)){
				ret=true;
				if(!in){
					in=true;
					pos[mid].push_back(sa[i]);
				}
			}
		}
		else
			memset(vis,0,sizeof(vis)),vis[belong[sa[i]]]=1,in=false,cnt=1;
			
	}
	return ret;
}

char in[maxn];
int main()
{
    int N;
    bool between = false;
    while(scanf("%d", &n) && n)
    {
        N = 0;
        if(between) printf("\n");
        else between = true;
        for(int i = 0; i <= 1000; i++) pos[i].clear();
        for(int i = 0; i < n; i++)
        {
            scanf("%s", in);
            int tmp = strlen(in);
            for(int j = 0; j < tmp; j++)
                belong[N] = i, s[N++] = in[j] - 'a' + 1;
            s[N++] = 27 + i;
        }
        N--;
        s[N] = 0;
        build(s, N + 1, 130);
        getheight( N);
        int L = 0, R = 1000, ans = 0, mid;
        while(L <= R)//二分查找最大的長度mid
        {
            mid = (L + R) >> 1;
            if(check(mid, N))
            {
                ans = mid;
                L = mid + 1;
            }
            else
                R = mid - 1;
        }
        if(ans == 0)
            printf("?\n");
        else
           
            for(unsigned int i = 0; i < pos[ans].size(); i++)
            {
                for(int j = 0; j < ans; j++)
                    printf("%c", s[pos[ans][i] + j] - 1 + 'a');
                printf("\n");
            }
    }
    return 0;
}