HDU - 1403 - Longest Common Substring

先上题目：

Longest Common Substring

Time Limit: 8000/4000 MS (Java/Others)    Memory Limit: 65536/32768 K (Java/Others)
Total Submission(s): 4010    Accepted Submission(s): 1510

Problem Description

Given two strings, you have to tell the length of the Longest Common Substring of them.

For example:
str1 = banana
str2 = cianaic

So the Longest Common Substring is "ana", and the length is 3.

 

Input

The input contains several test cases. Each test case contains two strings, each string will have at most 100000 characters. All the characters are in lower-case.

Process to the end of file.

 

Output

For each test case, you have to tell the length of the Longest Common Substring of them.

 

Sample Input

banana

cianaic

 

Sample Output

3

 

　　题意：给出两个串，问你这两个串的最长公共子串的长度是多少。

　　后缀数组入门题。首先，不得不承认，现在我的水平只可以套一下模板，通过模板我们可以求出sa[],rank[],height[]三个数组。

　　对于这里的字符串，我们是从0~n-1。

　　sa[i]指的是字典序排第i的后缀的下标是什么，rank[i]指的是原串中第i个后缀(就是从第i个字符开始到末尾的字符串)在后缀数组中排第几。height[i]表示后缀数组中第i个后缀和第i-1一个后缀的最长公共前缀的长度是多少(其中height[0]=0)。

　　这里的做法是首先将两个字符串连接起来，在连接处加一个连接符(没在这两个字符串中出现过的字符即可)，然后求出height[]，再扫描height[]，寻找某个同时符合以下要求的值：①比最大值还要大，②suffex(sa[i])和suffex(sa[i-1])分属于两个不同的字符串。这里需要注意每个数组的长度足够。

 

上代码：

 

 #include <cstdio>
 #include <cstring>
 #include <algorithm>
 #define MAX 100002
 using namespace std;

 char s[(MAX<<)],b[MAX];
 int sa[MAX<<],rank[MAX<<],height[MAX<<],t[MAX<<],t2[MAX<<],c[MAX<<],n,li;
 int f[(MAX<<)];

 void build_sa(int m){
     int i,*x=t,*y=t2;
     for(i=;i<m;i++) c[i]=;
     for(i=;i<n;i++) c[x[i]=s[i]]++;
     for(i=;i<m;i++) c[i]+=c[i-];
     for(i=n-;i>=;i--) sa[--c[x[i]]]=i;
     for(int k=;k<=n;k<<=){
         int p=;
         for(i=n-k;i<n;i++) y[p++]=i;
         for(i=;i<n;i++) if(sa[i]>=k) y[p++]=sa[i]-k;
         for(i=;i<m;i++) c[i]=;
         for(i=;i<n;i++) c[x[y[i]]]++;
         for(i=;i<m;i++) c[i]+=c[i-];
         for(i=n-;i>=;i--) sa[--c[x[y[i]]]]=y[i];
         swap(x,y);
         p=;    x[sa[]]=;
         for(i=;i<n;i++){
             x[sa[i]]=y[sa[i-]]==y[sa[i]]&&y[sa[i-]+k]==y[sa[i]+k] ? p- : p++;
         }
         if(p>=n) break;
         m=p;
     }
     for(i=;i<n;i++) rank[sa[i]]=i;
 }

 void getHeight(){
     int i,j,k=;
     for(i=;i<n;i++){
         if(k) k--;
         j=sa[rank[i]-];
         while(s[i+k]==s[j+k]) k++;
         height[rank[i]]=k;
     }
 }

 int main()
 {
     int maxn;
     //freopen("data.txt","r",stdin);
     while(scanf("%s %s",s,b)!=EOF){
         strcat(s,"&");
         li=strlen(s);
         for(int i=;i<li-;i++) f[i]=;
         strcat(s,b);
         f[li-]=;
         n=strlen(s);
         for(int i=li;i<n;i++) f[i]=-;
         build_sa();
         getHeight();
         maxn=;
         for(int i=;i<n;i++){
             if(maxn<height[i] && f[sa[i-]]*f[sa[i]]<){
                 maxn=height[i];
             }
         }
         printf("%d\n",maxn);
     }
     return ;
 }

/*1403*/