利用滚动数组解决Longest common subsequence(LCS)问题

    Longest common subsequence问题:已知两序列,求这两个序列的最长公共子序列(不一定要连续)的长度。设二维数组dp[i][j]表示长度分别为i和j的序列A和B的LCS的最大长度,有状态转移方程:A[i]=B[j]时,dp[i][j]=dp[i-1][j-1]+1;A[i]≠B[j]时,dp[i][j]=max{dp[i-1][j],dp[i][j-1]}。
    通过观察发现,状态dp[i][j]只与当前行(dp[i][j-1])和上一行(dp[i-1][j],dp[i-1][j-1])的状态有关。这样,我们可以设置2个一维数组c1[],c2[],其中c1[]保存上一行的状态信息,c2[]更新当前行的状态信息,更新完毕之后再将c2[]的状态信息复制到c1[]中,循环结束后c1[length_B]就是所求LCS的最大长度,从而降低了空间需求。
 1 #include <iostream>
 2 using namespace std;
 3 
 4 const int MAXN = 5001;
 5 char str1[MAXN],str2[MAXN];
 6 int i,j,len1,len2,c1[MAXN],c2[MAXN];
 7 
 8 int main(){
 9     while(scanf("%s",str1)!=EOF){
10         scanf("%s",str2);
11         len1=strlen(str1),len2=strlen(str2);
12         for(i=0;i<=len2;i++) c1[i]=0;
13         for(i=0;i<=len1;i++){
14             for(j=0;j<=len2;j++){
15                 if(str1[i]==str2[j]) c2[j+1]=c1[j]+1;
16                 else c2[j+1]=c2[j]>c1[j+1? c2[j]:c1[j+1];
17             }
18             for(j=1;j<=len2;j++) c1[j]=c2[j];
19         }
20         printf("%d\n",c1[len2]);
21     }
22     return 0;
23 }

posted on 2009-04-22 16:08 极限定律 阅读(1152) 评论(4)  编辑 收藏 引用 所属分类: ACM/ICPC

评论

# re: 利用滚动数组解决Longest common subsequence(LCS)问题 2009-07-30 15:09 zc1991211

偶像!我以后一定多来拜访~  回复  更多评论   

# re: 利用滚动数组解决Longest common subsequence(LCS)问题 2010-10-24 16:22 rothchild

这样子效率是不是有点低...如果两个字符串都长度都上万...复制的时候就占用了好长时间...  回复  更多评论   

# re: 利用滚动数组解决Longest common subsequence(LCS)问题[未登录] 2011-08-28 09:18 -

n*m的算法上万本身就会超时@rothchild
  回复  更多评论   

# re: 利用滚动数组解决Longest common subsequence(LCS)问题 2012-10-26 20:16 蓝天下

可以用memcpy,复制速度比for快好几倍。
好像更常用的滚动数组做法是在二维数组中递推,对i的奇偶进行讨论,计算f[i%2][..]什么什么的。。  回复  更多评论   


只有注册用户登录后才能发表评论。
网站导航: 博客园   IT新闻   BlogJava   知识库   博问   管理


<2009年4月>
2930311234
567891011
12131415161718
19202122232425
262728293012
3456789

导航

统计

常用链接

留言簿(10)

随笔分类

随笔档案

友情链接

搜索

最新评论

阅读排行榜

评论排行榜