spoj 220 Relevant Phrases of Annihilation (后缀数组应用)

题意：

给出n个串，求着n个串最长公共子串，并且子串在各个串中至少重复两次，并且不能重叠的重复。

题解：

这里又用到了后缀的分组，我们首先拼接n个串然后跑一边后缀数组，二分长度，接着根据height分组，然后统计每组的后缀所属的串的个数，并且要判断是否存在至少重复两次的子串且不重叠，统计所属串的个数可以用set存，然后获取size即可，对于后一个问题，我们将每次分得的组统计所属每个串的SA[i]，得到每个串SA[i]的最大值最小值的差，如果满足差>=k（k是我们枚举的长度）那么这个方案成立。

注意每次组统计之前都要初始化最大值最小值，结果忘了，wa了n久。

#include<iostream>
#include<math.h>
#include<stdio.h>
#include<algorithm>
#include<string.h>
#include<vector>
#include<queue>
#include<map>
#include<set>
using namespace std;
#define B(x) (1<<(x))
typedef long long ll;
const int oo=0x3f3f3f3f;
const ll OO=1LL<<61;
const int MOD=10007;
const int maxn=200000;
int rank[maxn],SA[maxn],height[maxn];
int t1[maxn],t2[maxn],t3[maxn],t4[maxn];
char str[maxn];
int ss[maxn],which[maxn],n;
set<int>g;
int mx[12],mi[12];

void Swap(int*& x,int*& y){

    int *temp=x;
    x=y;
    y=temp;
}

int cmp(int t[],int a,int b,int l){

    return t[a]==t[b]&&t[a+l]==t[b+l];
}

void build_SA(int s[],int len,int up){

    int *k1=t1,*k2=t2,*r=t3,*cnt=t4;
    for(int i=0;i<up;i++)cnt[i]=0;
    for(int i=0;i<len;i++)cnt[k1[i]=s[i]]++;
    for(int i=1;i<up;i++)cnt[i]+=cnt[i-1];
    for(int i=len-1;i>=0;i--)SA[--cnt[k1[i]]]=i;

    for(int d=1,p=1;p<len;d<<=1,up=p){

        p=0;
        for(int i=len-d;i<len;i++)k2[p++]=i;
        for(int i=0;i<len;i++)if(SA[i]>=d)k2[p++]=SA[i]-d;
        for(int i=0;i<len;i++)r[i]=k1[k2[i]];

        for(int i=0;i<up;i++)cnt[i]=0;
        for(int i=0;i<len;i++)cnt[r[i]]++;
        for(int i=1;i<up;i++)cnt[i]+=cnt[i-1];
        for(int i=len-1;i>=0;i--)SA[--cnt[r[i]]]=k2[i];

        Swap(k1,k2);
        k1[SA[0]]=0;
        p=1;
        for(int i=1;i<len;i++){
            k1[SA[i]]= cmp(k2,SA[i-1],SA[i],d) ? p-1 : p++;
        }
    }
}

void get_height(int s[],int len){

    for(int i=1;i<=len;i++)rank[SA[i]]=i;
    for(int i=0,p=0;i<len;i++){

        int j=SA[rank[i]-1];
        while(s[i+p]==s[j+p])p++;
        height[rank[i]]=p;
        if(p)p--;
    }
}

void cmin(int& a,int b){
    if(a==-1||b<a)a=b;
}

void cmax(int& a,int b){
    if(a==-1||b>a)a=b;
}

bool ok(int tag,int len){

    int i=0;
    while(1){

        while(i<=len&&height[i]<tag)i++;
        if(i>len)break;
        memset(mx,-1,sizeof mx);
        memset(mi,-1,sizeof mi);
        g.clear();
        g.insert(which[SA[i-1]]);
        cmax(mx[which[SA[i-1]]],SA[i-1]);
        cmin(mi[which[SA[i-1]]],SA[i-1]);
        while(i<=len&&height[i]>=tag){

            g.insert(which[SA[i]]);
            cmax(mx[which[SA[i]]],SA[i]);
            cmin(mi[which[SA[i]]],SA[i]);
            i++;
        }
        if(g.size()>=n){

            int f=1;
            for(int j=1;j<=n;j++){
                if(mx[j]==-1||mi[j]==-1||mx[j]-mi[j]<tag){
                    f=0;
                    break;
                }
            }
            if(f)return true;
        }
    }
    return false;
}

int solve(int len){

    g.clear();
    int l=1,r=len,mid,ans=0;
    while(l<=r){

        mid=(r+l)>>1;
        if(ok(mid,len)){
            ans=mid;
            l=mid+1;
        }else{
            r=mid-1;
        }
    }
    return ans;
}

int main(){

    int l,T,tlen;
    scanf("%d",&T);
    while(T--){

        scanf("%d",&n);
        l=0;
        for(int i=1;i<=n;i++){

            scanf("%s",str);
            tlen=strlen(str);
            for(int j=0;j<tlen;j++){

                which[l]=i;
                ss[l++]=str[j]-'a'+1;
            }
            which[l]=i;
            ss[l++]=i+30;
        }
        l--;
        ss[l]=0;
        build_SA(ss,l+1,300);
        get_height(ss,l);
        printf("%d\n",solve(l));
    }
    return 0;
}
/**
100
2
adsf
asdf
2
fdsa
df
5
abcd
bcde
cdef
defg
efgh
5
nshiyedashai
goanshab
woshidiyi
firstorloser
iamthefirst
4
acaaaaaaaaaaaaa
aaaaa
aa
a
4
aa
a
d
d
2
pabcp
pertp
*/