URAL 1517. Freedom of choice(后缀数组:最长公共连续子串)

URAL 1517. Freedom of choice(后缀数组:最长公共连续子串)

题目链接：点击打开链接

题意:

        给你两个等长的串,求他们的最长公共连续子串.

分析: 

        罗穗骞《后缀数组——处理字符串的有力工具》例题.

        其实本题就是POJ2774的略微加强版:

       poj2774博客题解

        这里要求我们输出这个最长公共串了.(题目中要求有多解就输出所有解,但是其实只有1个解.)

        对于输出解,我们只要记录过程中的每一个可以作为解的后缀编号以及长度即可一一输出.

AC代码:

#include<iostream>
#include<cstdio>
#include<cstring>
#include<algorithm>
using namespace std;
const int maxn=200000+100;
const int maxm=200000+100;
struct SuffixArray
{
    char s[maxn];
    int sa[maxn],rank[maxn],height[maxn];
    int t1[maxn],t2[maxn],c[maxm],n;
    void build_sa(int m)
    {
        int i,*x=t1,*y=t2;
        for(i=0;i<m;i++) c[i]=0;
        for(i=0;i<n;i++) c[x[i]=s[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(int k=1;k<=n;k<<=1)
        {
            int p=0;
            for(i=n-k;i<n;i++) y[p++]=i;
            for(i=0;i<n;i++)if(sa[i]>=k) y[p++]=sa[i]-k;
            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]]= y[sa[i]]==y[sa[i-1]]&&y[sa[i]+k]==y[sa[i-1]+k]? p-1:p++;
            if(p>=n) break;
            m=p;
        }
    }
    void build_height()
    {
        int i,j,k=0;
        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(s[i+k]==s[j+k]) k++;
            height[rank[i]]=k;
        }
    }
}sa;
void solve(int len,char a[])
{
	int ans=0;
	int pos;
	for(int i=2;i<sa.n;i++)
	{
		int a1=sa.sa[i-1],a2=sa.sa[i];
		if(a1>a2)
			swap(a1,a2);
		if(a1>=0&&a1<=len-1&&a2>=len+1&&a2<=len*2)
			if(ans<sa.height[i])
			{
				ans=sa.height[i];
				pos=a1;
			}
	}
	for(int i=pos;i<=ans+pos-1;i++)
		printf("%c",a[i]);
	printf("\n");
}
int main()
{
    char a[maxn],b[maxn];
    int n;
    while(~scanf("%d",&n)&&n)
    {
    	scanf("%s%s",a,b);
    	for(int i=0;i<n;i++)
			sa.s[i]=a[i];
		sa.s[n]='@';
		for(int i=n+1;i<=2*n;i++)
			 sa.s[i]=b[i-n-1];
		sa.s[2*n+1]=0;
		sa.n=2*n+2;
		sa.build_sa(128);
		sa.build_height();
		solve(n,a);
    }
    return 0;
}