当前位置:   article > 正文

公共子序列的个数_公共子序列个数

公共子序列个数

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

给你两个数组,求公共子序列的个数。
比赛比了一会,就都有事撤了...现在开始补题,转移方程有点难想,千辛万苦想到了我还是WA了. 最后看了一下卿学姐@qscqesze的博客 ,才发现忘记处理如果dp[i][j]出现负的情况..
dp[i][j]代表A前i个数和B前j个数的公共子序列的个数。
那么dp[i][j]=dp[i][j-1]+dp[i-1][j]-dp[i-1][j-1];     因为dp[i][j-1] 为i ,j-1的公共子序列的个数,dp[i-1][j]为i-1,j的公共子序列的个数,所以两个加起来-重复的i-1,j-1的公共子序列的个数即为dp[i][j]的公共子序列个数,当a[i]==b[j] 时,还需要dp[i][j]+=dp[i-1][j-1]+1;即 当子序列都包含a[i]的时候那就是要多出i-1,j-1的个数,和仅有a[i],a[j]的子序列。
  1. #include <cstdio>
  2. #include <cstring>
  3. #include <iostream>
  4. using namespace std;
  5. #define Max(a,b) a>b?a:b
  6. int dp[1005][1005];
  7. int a[1005],b[1005];
  8. const int mod=1000000007;
  9. int main()
  10. {
  11. int n,m;
  12. while(~scanf("%d%d",&n,&m))
  13. {
  14. for(int i=1;i<=n;i++)
  15. scanf("%d",&a[i]);
  16. for(int i=1;i<=m;i++)
  17. scanf("%d",&b[i]);
  18. for(int i=0;i<=n;i++)
  19. for(int j=0;j<=m;j++)
  20. dp[i][j]=0;
  21. for(int i=1;i<=n;i++)
  22. for(int j=1;j<=m;j++)
  23. {
  24. dp[i][j]=dp[i-1][j]+dp[i][j-1]-dp[i-1][j-1];// 可能出现负
  25. if(a[i]==b[j]) dp[i][j]+=dp[i-1][j-1]+1;
  26. if(dp[i][j]<0) dp[i][j]+=mod;
  27. if(dp[i][j]>=mod) dp[i][j]%=mod;
  28. }
  29. printf("%d\n",dp[n][m]);
  30. }
  31. }



声明:本文内容由网友自发贡献,转载请注明出处:【wpsshop博客】
推荐阅读
相关标签
  

闽ICP备14008679号