【poj 2774】Long Long Message 中文题意&题解&代码(C++)
2016-03-09 21:48
411 查看
题目链接:
http://poj.org/problem?id=2774
题意:
给出两个字符串,求这两个字符串的最长公共子串的长度。
题解:
利用后缀数组的性质,将两个字符串连接成一个长串,然后求出这个新串的height数组,因为height数组表示后缀排名为i和i-1的后缀的最长前缀,所以for循环一遍取height数组中的最大值即可,因为按字典序来排名,排名最近的字符串一定是最相近的字符串,所以这样做是对的。
不过要注意当循环到 i 位置时,要判断sa[i]和sa[i-1]是否分别属于两个字符串,这样求出的才是这两个字符串的最长公共子串。
代码:
http://poj.org/problem?id=2774
题意:
给出两个字符串,求这两个字符串的最长公共子串的长度。
题解:
利用后缀数组的性质,将两个字符串连接成一个长串,然后求出这个新串的height数组,因为height数组表示后缀排名为i和i-1的后缀的最长前缀,所以for循环一遍取height数组中的最大值即可,因为按字典序来排名,排名最近的字符串一定是最相近的字符串,所以这样做是对的。
不过要注意当循环到 i 位置时,要判断sa[i]和sa[i-1]是否分别属于两个字符串,这样求出的才是这两个字符串的最长公共子串。
代码:
#include<iostream> #include<algorithm> #include<stdio.h> #include<string.h> #define maxn (200005) using namespace std; int lim,ans=-1,sum,wa[maxn],wb[maxn],wv[maxn],R[maxn],SA[maxn],WS[maxn]; char s[maxn]; int cmp(int *r,int a,int b,int l) {return r[a]==r[b]&&r[a+l]==r[b+l]; } void da(int *r,int *sa,int n,int m) { int *x=wa,*y=wb,*t; int i,j,p; for (i=0;i<m;i++) WS[i]=0; for (i=0;i<n;i++) WS[x[i]=r[i]]++; for (i=1;i<m;i++) WS[i]+=WS[i-1]; for (i=n-1;i>=0;i--) sa[--WS[x[i]]]=i; for (j=1,p=1;p<n;j*=2,m=p) { for (p=0,i=n-j;i<n;i++) y[p++]=i; for (i=0;i<n;i++) if (sa[i]-j>=0) y[p++]=sa[i]-j; for (i=0;i<n;i++) wv[i]=x[y[i]]; for (i=0;i<m;i++) WS[i]=0; for (i=0;i<n;i++) WS[wv[i]]++; for (i=1;i<m;i++) WS[i]+=WS[i-1]; for (i=n-1;i>=0;i--) sa[--WS[wv[i]]]=y[i]; for (t=x,x=y,y=t,p=1,x[sa[0]]=0,i=1;i<n;i++) x[sa[i]]=cmp(y,sa[i-1],sa[i],j)?p-1:p++; } } int rank[maxn],ht[maxn]; void get_h(int *r,int *sa,int n) { int i,j,k=0; for (i=1;i<=n;i++) rank[sa[i]]=i; for (i=0;i<n;ht[rank[i++]]=k) for (k?k--:0,j=sa[rank[i]-1];r[i+k]==r[j+k];k++); return ; } int main() { sum=0; scanf("%s",s); int len=strlen(s); for (int i=0;i<len;i++) R[sum]=s[i]-'a'+1,sum++; R[sum]=0; sum++; lim=len; scanf("%s",s); len=strlen(s); for (int i=0;i<len;i++) R[sum]=s[i]-'a'+1,sum++; R[sum]=0;sum++; da(R,SA,sum,30); get_h(R,SA,sum-1); // for (int i=0;i<sum;i++) // cout<<ht[i]<<' '; // cout<<endl; for (int i=1;i<sum;i++) if ((SA[i]>lim&&SA[i-1]<lim)||(SA[i]<lim&&SA[i-1]>lim)) ans=max(ans,ht[i]); printf("%d\n",ans); }
相关文章推荐
- contor expansion
- 正态分布,泊松分布,指数分布的c/c++代码
- 2015年第六届蓝桥杯试题(C/C++本科B组)
- 第2周项目2-就拿胖子说事
- c++上机报告1-1
- C++的异常处理
- c语言:利用静态变量static,输出1到5的阶乘值
- C++ vector变量等导致内存泄露问题的解决方法
- CAF(C++ actor framework)(序列化之结构体,任意嵌套STL)(一)
- 判断任意数据是否奇偶数
- 装饰模式(Decorator)C++实现
- C++11 lambda表达式
- c++多线程编程(三)
- c++多线程编程(二)
- c++多线程编程(一)
- C++中友元函数
- c++11 可变参数模板
- C语言实现C++面向对象语言多态特性
- 第2周项目1
- c语言 谭浩强 一维数组内放10个学生成绩 全局变量写一个函数 最高分、最低分、平均分