UVA 11107 Life Forms（后缀数组+二分）

思路：

后缀数组+二分：

我们先把n 个字符串用不同的字符连接起来，注意最后要加上一个字符。

然后求这个字符串的后缀数组。

然后二分答案，判断长度len 是否可行。 判断方法是同样的套路，给sa数组进行分块，每一块都大于等于长度len，然后判断每一块是否有超过一半的串在这一块中出现， 可以开一个flag数组， 记录每一块的字符串， 最后统计flag数组即可。

这样我们二分得到了最长长度len。

然后打印解。

同样的方式 枚举每一块，如果这一块合法，  因为这一块每一个字符串 前len 个字符都是一样的， 直接取前len个字符就是答案。

最后排序即可。

#include <cstdio>
#include <cstring>
#include <algorithm>
#include <vector>
#include <string>
#define Siz(x) (int)x.size()
using namespace std;

const int maxn = 1000000 + 10;

int t1[maxn], t2[maxn], c[maxn];

bool cmp(int* r, int a,int b,int l){
return r[a] == r[b] && r[a+l] == r[b+l];
}

void da(int str[], int sa[], int Rank[], int lcp[], int n, int m){
++n;
int i, j, p, *x = t1, *y = t2;
for (i = 0; i < m; ++i) c[i] = 0;
// puts("hha");
for (i = 0; i < n; ++i) c[x[i] = str[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 (j = 1; j <= n; j <<= 1){
p = 0;
for (i = n-j; i < n; ++i) y[p++] = i;
for (i = 0; i < n; ++i) if (sa[i] >= j) y[p++] = sa[i] - j;
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] ] = cmp(y, sa[i-1], sa[i], j) ? p-1 : p++;

}

if (p >= n)break;
m= p;

}

int k = 0;
n--;
for (i = 0; i <= n; ++i) Rank[sa[i] ] = i;
for (i = 0; i < n; ++i){
if (k)--k;
j = sa[Rank[i]-1 ];
while(str[i+k] == str[j+k])++k;
lcp[Rank[i] ] = k;
}
}

int lcp[maxn], a[maxn], sa[maxn], Rank[maxn];
int n,k;

char s[1000 + 10];

vector<int>bit;
int cnt, id;
int flag[100 + 10];
int nn;
bool judge(int m){
int cp = n / 2;
bit.clear();
bit.push_back(sa[1]);
for (int i = 2; i <= cnt; ++i){
if (lcp[i] >= m) bit.push_back(sa[i]);
else {
memset(flag,0,sizeof flag);
for (int j = 0; j < Siz(bit); ++j){
int v = bit[j] / (nn+1);
flag[v] = 1;
}
int sum = 0;
for (int j = 0; j < n; ++j) sum += flag[j];
if (sum > cp) return true;
bit.clear();
bit.push_back(sa[i]);
}
}

memset(flag,0,sizeof flag);
for (int j = 0; j < Siz(bit); ++j){
int v = bit[j] / (nn+1);
flag[v] = 1;
}
int sum = 0;
for (int j = 0; j < n; ++j) sum += flag[j];
if (sum > cp) return true;
return false;
}

int ks = 0;
vector<string>Vec;

void solve(int len){
if (len == 0){
puts("?");
return;
}

int cp = n / 2;
bit.clear();
bit.push_back(sa[1]);
for (int i = 2; i <= cnt; ++i){
if (lcp[i] >= len) bit.push_back(sa[i]);
else {
string tmp;
for (int j = bit[0], k= 0; k < len; ++k, ++j) {
tmp += a[j] - 1 + 'a';

}
int kk = Siz(bit);
memset(flag,0,sizeof flag);
for (int j = 0; j < Siz(bit); ++j){
int v = bit[j] / (nn+1);
flag[v] = 1;
}
int sum = 0;
for (int j = 0; j < n; ++j) sum += flag[j];
if (sum > cp){
string ans;
for (int i = bit[0], j = 0; j < len; ++j,++i){
char ch = a[i]-1+'a';
ans += ch;
}
Vec.push_back(ans);
}

bit.clear();
bit.push_back(sa[i]);
}

}
memset(flag,0,sizeof flag);
for (int j = 0; j < Siz(bit); ++j){
int v = bit[j] / (nn+1);
flag[v] = 1;
}
int sum = 0;
for (int j = 0; j < n; ++j) sum += flag[j];
if (sum > cp){
string ans;
for (int i = bit[0], j = 0; j < len; ++j,++i){
char ch = a[i]-1+'a';
ans += ch;
}
Vec.push_back(ans);
}

}

int main(){

while(~scanf("%d",&n) && n){
Vec.clear();
id=26, cnt = 0;
int len;
for (int i = 0; i < n; ++i){
scanf("%s",s);
if (i == 0) nn = strlen(s);
for (int j = 0; s[j] ;++j){
a[cnt++] = s[j] - 'a' + 1;
}
a[cnt++] = ++id;
}

a[cnt] = 0;
da(a, sa, Rank, lcp, cnt, 200);
int l = 1, r = cnt, m;
while(l <= r){
m = l + r >> 1;
if (judge(m)) l = m + 1;
else r = m - 1;
}
if (ks++) putchar('\n');
solve(r);
sort(Vec.begin(), Vec.end());
for (int i = 0; i < Siz(Vec); ++i){
printf("%s\n", Vec[i].c_str());
}
}

return 0;
}