CodeForces 427 D.Match & Catch（后缀数组）

Description
给出两个串s1和s2，要求找到一个最短的s1和s2的公共子串使得该子串在s1和s2中均只出现过一次
Input
两个字符串s1和s2(|s1|,|s2|<=5000)
Output
输出满足条件的最短子串的长度，如果不存在这样的子串则输出-1
Sample Input
apple
pepperoni
Sample Output
2
Solution
把两个串连起来用后缀数组得到按字典序排好的后缀，一个长度len是满足条件的当且仅当存在两个相邻的后缀sa[i]和sa[i+1]，使得这两个后缀分别输入s1和s2，且这两个后缀的height值大于等于len，而sa[i-1]和sa[i]以及sa[i+1]和sa[i+2]的height值小于len，从小到大枚举len然后扫一遍后缀即可
Code

#include<cstdio>
#include<iostream>
#include<cstring>
#include<algorithm>
#include<cmath>
#include<vector>
#include<queue>
#include<map>
#include<set>
#include<ctime>
using namespace std;
typedef long long ll;
#define INF 0x3f3f3f3f
#define maxn 11111
int t1[maxn],t2[maxn],sa[maxn],c[maxn],Rank[maxn],height[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 n,int m)
{
    n++;
    int i,j,p,*x=t1,*y=t2;
    for(i=0;i<m;i++)c[i]=0;
    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++;
            height[Rank[i]]=k;
        }
}
char s[maxn];
int a[maxn],len1,len2;
int belong(int x)
{
    if(x>=0&&x<len1)return 1;
    if(x>len1&&x<len2)return 2;
    return 0;
}
bool check(int n,int k)
{
    for(int i=2;i<=n;i++)
    {
        int x=belong(sa[i-1]),y=belong(sa[i]);
        if(x&&y&&x!=y)
        {
            if(i==2&&height[i]>=k&&height[i+1]<k)return 1;
            else if(i==n&&height[i]>=k&&height[i-1]<k)return 1;
            else if(height[i]>=k&&height[i-1]<k&&height[i+1]<k)return 1;
        }
    }
    return 0;
}
int main()
{
    while(~scanf("%s",s))
    {
        len1=strlen(s);
        s[len1]='1';
        scanf("%s",s+len1+1);
        len2=strlen(s);
        for(int i=0;i<len2;i++)a[i]=s[i];
        a[len2]=0;
        da(a,len2+1,128);
        int ans=-1;
        for(int i=1;i<=len1;i++)
            if(check(len2+1,i))
            {
                ans=i;
                break;
            }
        printf("%d\n",ans);
    }
    return 0;
}