HDU 1403 Longest Common Substring（最长公共前缀）

http://acm.hdu.edu.cn/showproblem.php?pid=1403

题意：
给出两个字符串，求最长公共子串的长度。

思路：

刚开始学后缀数组，确实感觉很难，但是这东西很强大，所以必须要学会它，推荐罗穗骞大牛的论文。

 1 #include<iostream>
 2 #include<algorithm>
 3 #include<cstring>
 4 #include<cstdio>
 5 #include<vector>
 6 #include<stack>
 7 #include<queue>
 8 #include<cmath>
 9 #include<map>
10 #include<set>
11 using namespace std;
12 typedef long long ll;
13 typedef pair<int,int> pll;
14 const int INF = 0x3f3f3f3f;
15 const int maxn=200000+5;
16
17 char s[maxn],s1[maxn];
18 int sa[maxn],t[maxn],t2[maxn],c[maxn],n;
19 int Rank[maxn],height[maxn];
20
21 void build_sa(int m)
22 {
23     int *x=t,*y=t2;
24     //基数排序
25     for(int i=0;i<m;i++)    c[i]=0;
26     for(int i=0;i<n;i++)    c[x[i]=s[i]]++;
27     for(int i=1;i<m;i++)    c[i]+=c[i-1];
28     for(int i=n-1;i>=0;i--) sa[--c[x[i]]]=i;
29     for(int k=1;k<=n;k<<=1)
30     {
31         int p=0;
32         //直接利用sa数组排序第二关键字
33         for(int i=n-k;i<n;i++)  y[p++]=i;
34         for(int i=0;i<n;i++)    if(sa[i]>=k)    y[p++]=sa[i]-k;
35         //基数排序第一关键字
36         for(int i=0;i<m;i++)    c[i]=0;
37         for(int i=0;i<n;i++)    c[x[y[i]]]++;
38         for(int i=1;i<m;i++)    c[i]+=c[i-1];
39         for(int i=n-1;i>=0;i--) sa[--c[x[y[i]]]]=y[i];
40         //根据sa和y计算新的x数组
41         swap(x,y);
42         p=1;
43         x[sa[0]]=0;
44         for(int i=1;i<n;i++)
45             x[sa[i]]=y[sa[i-1]]==y[sa[i]]&&y[sa[i-1]+k]==y[sa[i]+k]?p-1:p++;
46         if(p>=n)
47             break;
48         m=p;                //下次基数排序的最大值
49     }
50 }
51
52 void getHeight()
53 {
54     int i,j,k=0;
55     for(i=1;i<=n;i++)  Rank[sa[i]]=i;
56     for(i=0;i<n;i++)
57     {
58         if(k)  k--;
59         int j=sa[Rank[i]-1];
60         while(s[i+k]==s[j+k])  k++;
61         height[Rank[i]]=k;
62     }
63 }
64
65 int main()
66 {
67     //freopen("in.txt","r",stdin);
68     while(~scanf("%s",s))
69     {
70         int len = strlen(s);
71         s[len]=‘0‘;
72         scanf("%s",s+len+1);
73         n=strlen(s);
74         build_sa(128);
75         getHeight();
76         int ans=0;
77         for(int i=1;i<n;i++)
78         {
79             if(height[i]>ans)
80             {
81                 if((sa[i]>len && sa[i-1]<len) || (sa[i]<len && sa[i-1]>len))
82                     ans=height[i];
83             }
84         }
85         printf("%d\n",ans);
86     }
87     return 0;
88 }