poj3294 Life Forms

题面在这里

题意：

给你n个串，找出最长的出现在超过n/2个串中的公共子串。

做法：

首先二分一个长度x，然后问题转化为如何检验是否存在超过n/2个串的公共子串长度为x。
对于公共子串这个问题，马上想到后缀数组，因为height[]数组就可以巧妙地解决这个问题。
首先把所有串用一些不同的没有出现过的字符连接起来。跑这个大串的sa。
把所有后缀按照height分组，把height>=x的连续一段取出来，那么这一段任意两个串肯定存在一个长度>=x的公共子串。
于是我们只要看这一段里面是否存在超过>n/2个出现在不同串里的后缀就可以了。

易错点：

不要背错板子，排序写成了顺序循环（背错板子这种sb错误应该也只有我会犯了。。）

代码：

/*************************************************************
    Problem: poj 3294 Life Forms
    User: bestFy
    Language: C++
    Result: Accepted
    Time: 938MS
    Memory: 3436K
    Submit_Time: 2018-01-18 10:58:36
*************************************************************/

#include<cstdio>
#include<cstring>
#include<iostream>
#include<algorithm>
#include<cmath>
#include<cstdlib>
#include<cctype>
using namespace std;
typedef long long ll;

const int N = 200010;
char str[1010];
int n, total, all;
int s[N], pos[N], sa[N], rk[N], tp[N], tong[N], h[N], vis[105];

inline void ssort() {
    for(int i = 0; i <= all; i ++) tong[i] = 0;
    for(int i = 1; i <= n; i ++) tong[rk[tp[i]]] ++;
    for(int i = 1; i <= all; i ++) tong[i] += tong[i-1];
    for(int i = n; i >= 1; i --) sa[tong[rk[tp[i]]] --] = tp[i];//倒序！！
}
inline void get_sa() {
    for(int i = 1; i <= n; i ++) rk[i] = s[i], tp[i] = i;
    all = 130; ssort(); int w = 1; all = 1;
    while(all < n) {
        int t = 0;
        for(int i = n-w+1; i <= n; i ++) tp[++ t] = i;
        for(int i = 1; i <= n; i ++) if(sa[i] > w) tp[++ t] = sa[i]-w;
        ssort(); for(int i = 1; i <= n; i ++) tp[i] = rk[i];
        rk[sa[1]] = all = 1;
        for(int i = 2; i <= n; i ++)
            rk[sa[i]] = (tp[sa[i]] == tp[sa[i-1]] && tp[sa[i]+w] == tp[sa[i-1]+w])?all:++ all;
        w <<= 1;
    } int k = 0;
    for(int i = 1; i <= n; i ++) {
        if(k) k --; int j = sa[rk[i]-1];
        for(; i+k <= n && j+k <= n && s[i+k] == s[j+k]; k ++);
        h[rk[i]] = k;
    }
}
inline void print(int x, int y) {
    for(int i = x; i <= x+y-1; i ++) putchar(s[i]-100+'a');
    puts("");
}
inline bool check2(int l, int r) {//检查这一段区间是否含有超过一半的串
    memset(vis, 0, sizeof vis); int sum = 0;
    for(int i = l; i <= r; i ++)
        if(!vis[pos[sa[i]]]) sum ++, vis[pos[sa[i]]] = 1;
    return 2*sum > total;
}
inline bool check(int x, bool flag = 0) {
    int st = 0, ed = 0;
    bool ok = 0;
    for(int i = 1; i <= n; i ++)
        if(h[i] >= x) ed ++;
        else {
            if(check2(st, ed)) {
                ok = 1;
                if(flag) print(sa[st], x);
            }
            st = ed = i;
        }
    if(st < ed) {
        if(check2(st, ed)) {
            ok = 1;
            if(flag) print(sa[st], x);
        }
    }
    return ok;
}
int main() {
    while(~scanf("%d", &total) && total) {
        n = 0; int l = 1, r = 0, mid;
        for(int i = 1; i <= total; i ++) {
            scanf("%s", str+1); int l = strlen(str+1); r = max(r, l);
            for(int j = 1; j <= l; j ++) s[++ n] = str[j]-'a'+100, pos[n] = i;
            s[++ n] = i, pos[n] = i;
        }
        get_sa(); int ans = -1;
        while(l <= r) {
            mid = l+r>>1;
            if(check(mid)) ans = mid, l = mid+1; else r = mid-1;
        }
        if(ans == -1) puts("?"); else check(ans, 1);
        puts("");
    }
    return 0;
}