POJ 题目2774 Long Long Message（后缀数组，求最长公共子串长度）

标签：

Description

Input

Output

Sample Input

yeshowmuchiloveyoumydearmotherreallyicannotbelieveit
yeaphowmuchiloveyoumydearmother

Sample Output

27

Source

POJ Monthly--2006.03.26,Zeyuan Zhu,"Dedicate to my great beloved mother."

ac代码

/*Problem: 2774  User: kxh1995 
Memory: 5024K  Time: 891MS 
Language: C++  Result: Accepted */

#include<stdio.h>   
#include<string.h>   
#include<algorithm>   
#include<iostream>  
#define min(a,b) (a>b?b:a) 
using namespace std;  
char str1[200100],str2[2000100];
int sa[220020],c[220020],t2[220020];
int t1[220020],s[200100];  
int rank[220020],height[220020]; 
int len1,len2; 
void build_sa(int s[],int n,int m)  
{  
    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]=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(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]]=y[sa[i-1]]==y[sa[i]]&&y[sa[i-1]+j]==y[sa[i]+j]?p-1:p++;  
        if(p>=n)  
            break;  
        m=p;  
    }  
}  
void getHeight(int s[],int n)  
{  
    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;  
    }  
}
int judge(int len,int k)
{
	int i;
	for(i=1;i<len;i++)
	{
		if(height[i]>=k)
		{
			if(sa[i]>len1&&sa[i-1]<=len1)
				return 1;
			if(sa[i-1]>len1&&sa[i]<=len1)
				return 1;
		}
	}
	return 0;
}
int main()
{
	while(scanf("%s%s",str1,str2)!=EOF)
	{
		int i;
		len1=strlen(str1);
		len2=strlen(str2);
		for(i=0;i<len1;i++)
		{
			s[i]=str1[i]-'a'+1;
		}
		s[len1]=27;
		int n=len1+1;
		for(i=0;i<len2;i++)
			s[n++]=str2[i]-'a'+1;
		s[n]=0;
		build_sa(s,len1+len2+3,28);
		getHeight(s,len1+len2+2);
		int l=0,r=min(len1,len2),ans=0;
		while(l<=r)
		{
			int mid=(l+r)>>1;
			if(judge(len1+len2+2,mid))
			{
				ans=mid;
				l=mid+1;
			}
			else
				r=mid-1;
		}
		printf("%d\n",ans);
	}
}

POJ 题目2774 Long Long Message（后缀数组，求最长公共子串长度）

标签：

原文地址：http://blog.csdn.net/yu_ch_sh/article/details/47755335