最长公共子串（后缀数组）



题目：http://poj.org/problem?id=2217

首先解释，DP中的最长公共子序列和此处的最长公共子串区别-------------------序列可以是不连续的，但是子串是连续的

其次，LCP，lcp[i]就是lcp[rank[i]]和lcp[rank[i]+1]的最长公共前缀，那么把两个字符串接起来，然后找最长的lcp，就是答案

思路还是比较清晰的

上代码：

[cpp] view plain copy print ?

1. /*******************************************************/  
2. //poj 2217 lcp+sa by Pilgrim  
3. //最长公共子串---注意与动态规划的最长公共子序列不同  
4. //2014.4.2  
5. /******************************************************/  
6.   
7. #include <string>  
8. #include <cstring>  
9. #include <cstdio>  
10. #include <algorithm>  
11. #include <iostream>  
12. #include <cstdlib>  
13. #include <vector>  
14.   
15. #define MAXN 10010  
16. #define INF 0x80000000  
17. //0x7fffffff  
18. //0x80000000  
19.   
20. using namespace std;  
21.   
22. int n,k;//n=strlen(s);  
23.   
24. int Rank[MAXN];  
25. int tmp[MAXN];  
26.   
27. /*使用Rank对sa排序*/  
28. bool cmpSa(int i, int j)  
29. {  
30.     if(Rank[i] != Rank[j])return Rank[i] < Rank[j];  
31.     else  
32.     {   /*下面的Rank[t]，已经是以t开头长度小于等于k/2的， 
33.         sa[i]的名次，只是以i开头的后缀，而长度不同*/  
34.         int ri = i+k <=n? Rank[i+k]:-1;  
35.         int rj = j+k <= n ? Rank[j+k]:-1;  
36.         return ri <rj;  
37.     }  
38. }  
39.   
40. /*计算SA*/  
41. void con_sa(char *s, int *sa)  
42. {  
43.     /*n=strlen(s);  必要时注明*/  
44.     /*初始化sa和rank保证两点 
45.         1、Rank[i]表示下标为i的是第几大，必须表示出相对大小，可以直接用字符代表其大小 
46.         2、sa[1...n]值为1..n*/  
47.     for(int i=0;i<=n;i++){  
48.         sa[i]=i;Rank[i] = i < n?s[i]:-1;  
49.     }  
50.   
51.     /*利用长度为k的字符串对长度为2*k的字符串排序*/  
52.     for(k=1;k<=n;k*=2)/*注意此代码中k是全局变量 别乱用,循环必须从1开始，因为0*2=0*/  
53.     {  
54.         sort(sa,sa+n+1,cmpSa);  
55.         tmp[sa[0]] = 0; /*此时tmp只是暂存rank*/  
56.         for(int i=1;i<=n;i++){  
57.             tmp[sa[i]] = tmp[sa[i-1]] +(cmpSa(sa[i-1],sa[i])?1:0);  
58.             /*这一句很关键，等号右侧的sa[i]在此循环里表示第i大的长度小于等于k/2的字符串， 
59.               从而求出第i大的长度小于等于k的字符串的sa[i]*/  
60.         }  
61.         for(int i=0;i<=n;i++){  
62.             Rank[i] = tmp[i];  
63.         }  
64.     }  
65. }  
66.   
67. void construct_lcp(char *s,int *sa,int *lcp)  
68. {  
69.     //n=strlen(s);  
70.     for(int i=0; i<=n; i++)Rank[sa[i]]=i;  
71.   
72.     int h=0;  
73.     lcp[0]=0;  
74.     for(int i=0;i<n;i++)  
75.     {  
76.         int j=sa[Rank[i]-1];  
77.   
78.         if(h>0)h--;  
79.         for(; j+h<n && i+h<n; h++)  
80.         {  
81.             if(s[j+h]!=s[i+h])break;  
82.         }  
83.         lcp[Rank[i]-1]=h;  
84.     }  
85. }  
86.   
87.   
88. int main()  
89. {  
90.     int sa[MAXN],lcp[MAXN];  
91.     char s[MAXN],t[MAXN];  
92.     char c;  
93.     int ncase,mmax,len,leng;  
94.   
95.     while(scanf("%d",&ncase)!=EOF)  
96.     {  
97.         while(ncase--)  
98.         {  
99.             mmax =0;  
100.             while((c=getchar())==' '|| c=='\n');  
101.             s[0]=c;  
102.             int tt=1;  
103.             while((c=getchar()) != '\n')  
104.             {  
105.                 s[tt++]=c;  
106.   
107.             }  
108.   
109.             s[tt]='\0';  
110.             while((c=getchar())==' '|| c=='\n');  
111.             t[0]=c;  
112.             tt=1;  
113.             while((c=getchar()) != '\n')  
114.             {  
115.                  t[tt++]=c;  
116.   
117.             }  
118.             t[tt]='\0';  
119.             len=strlen(s);  
120.             leng =len+1+strlen(t);  
121.             strcpy(s+len+1,t);  
122. ///  
123. //for(int i=0;i<leng;i++)  
124.  //   putchar(s[i]);  
125. //putchar('\n');  
126.             n=leng;  
127.             con_sa(s,sa);  
128.             construct_lcp(s,sa,lcp);  
129.   
130.            for(int i=0;i<leng;i++)  
131.             {  
132.               if((sa[i]<len) != (sa[i+1]<len))  
133.                   mmax=max(mmax,lcp[i]);  
134.             }  
135.             printf("Nejdelsi spolecny retezec ma delku %d.\n",mmax);  
136.         }  
137.     }  
138.   
139.     return 0;  
140. }