SPOJ

传送门：SPOJ - PHRASES（后缀数组+二分）

题意：给你n个字符串，找出一个最长的子串，他必须在每次字符串中都出现至少两次。

题解：被自己蠢哭...记录一下自己憨憨的操作，还一度质疑评测鸡（哭... 

首先是多个字符串的常规操作（目前写的题少，见到的都是这样）连成一个字符串，中间用不同的且没出现过的字符隔开。然后后缀数组求出sa数组和height数组，二分子串长度，用mx数组记录该串中的这个子串的起点最大值，mi数组记录最小值，如果所有串中的mx-mi都>=k说明这个长度可以。

#include<cstdio>
#include<algorithm>
#include<queue>
#include<iostream>
#include<cmath>
#include<cstring>
using namespace std;

//sa:字典序中排第i位的起始位置在str中第sa[i]  sa[1~n]为有效值

//rnk:就是str第i个位置的后缀是在字典序排第几 rnk[0~n-1]为有效值

//height：字典序排i和i-1的后缀的最长公共前缀  height[2~n]为有效值，第二个到最后一个

const int INF=0x3f3f3f3f;
const int maxn = 1e5+10;
int wa[maxn],wb[maxn],wsf[maxn],wv[maxn],sa[maxn];
int rnk[maxn],height[maxn],be[maxn];
char s[maxn];
int r[maxn],n,len,mx[maxn],mi[maxn];

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

void getsa(int *r,int *sa,int n,int m)//n为添加0后的总长
{
    int i,j,p,*x=wa,*y=wb,*t;
    for(i=0; i<m; i++)  wsf[i]=0;
    for(i=0; i<=n; i++)  wsf[x[i]=r[i]]++;
    for(i=1; i<m; i++)  wsf[i]+=wsf[i-1];
    for(i=n; i>=0; i--)  sa[--wsf[x[i]]]=i;
    p=1;
    j=1;
    for(; p<=n; j*=2,m=p){
        for(p=0,i=n+1-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<=n; i++)  wv[i]=x[y[i]];
        for(i=0; i<m; i++)  wsf[i]=0;
        for(i=0; i<=n; i++)  wsf[wv[i]]++;
        for(i=1; i<m; i++)  wsf[i]+=wsf[i-1];
        for(i=n; i>=0; i--)  sa[--wsf[wv[i]]]=y[i];
        swap(x,y);
        x[sa[0]]=0;
        for(p=1,i=1; i<=n; i++)
            x[sa[i]]=cmp(y,sa[i-1],sa[i],j)? p-1:p++;
    }
}

void getheight(int *r,int n)//n为添加0后的总长
{
    int i,j,k=0;
    for(i=1; i<=n; i++)  rnk[sa[i]]=i;
    for(i=0; i<n; i++){
        if(k)
            k--;
        else
            k=0;
        j=sa[rnk[i]-1];
        while(r[i+k]==r[j+k])
            k++;
        height[rnk[i]]=k;
    }
}

bool check(int k)
{
    for(int j=0;j<n;j++) mx[j]=-INF,mi[j]=INF;      //一开始忘记初始化
    for(int i=2;i<=len;i++){
        if(height[i]>=k){
            mx[be[sa[i]]]=max(mx[be[sa[i]]],sa[i]);
            mi[be[sa[i]]]=min(mi[be[sa[i]]],sa[i]);
            mx[be[sa[i-1]]]=max(mx[be[sa[i-1]]],sa[i-1]);
            mi[be[sa[i-1]]]=min(mi[be[sa[i-1]]],sa[i-1]);       //没眼睛写成了sa[i]
        }
        else{
            int j=0;
            for(j=0;j<n;j++) mx[j]=-INF,mi[j]=INF;
            mx[be[sa[i]]]=sa[i],mi[be[sa[i]]]=sa[i];
        }
        int j=0;
        for(j=0;j<n;j++){
            if(mx[j]-mi[j]<k) break;
        }
        if(j==n) return true;
    }
    return false;
}

int main()
{
    ios::sync_with_stdio(false);
    cin.tie(0);
    cout.tie(0);
    int t;
    cin>>t;
    while(t--){
        cin>>n;
        len=0;
        for(int i=0;i<n;i++){
            cin>>s+len;
            int p=strlen(s+len);
            for(int j=0;j<p;j++){
                r[j+len]=s[j+len]-'a'+1;
                be[len+j]=i;
            }
            len+=p;
            if(i!=n-1){
                be[len]=i;
                r[len++]=i+100;
            }
        }
        r[len]=0;
        getsa(r,sa,len,200);
        getheight(r,len);
        int l=0,r=min(len,10000);
        int ans=0;
        while(l<=r){
            int mid=l+r>>1;
            if(check(mid)){
                ans=mid;
                l=mid+1;
            }
            else r=mid-1;
        }
        cout<<ans<<endl;
    }
    return 0;
}

SPOJ - PHRASES