【动态规划】最长公共上升子序列------------------------------dp

第二天叫醒我的不是闹钟,是梦想!

题目描述
研究发现,大猩猩的基因序列和人的基因序列只有1.3%的区别,更进一步,不仅仅离人最近的大猩猩和人的基因序列高度近似,就连以打洞为生的老鼠和人的基因序列也有高达95%的相同序列。于是有魔法师提出一个大胆设想,即改变人类的某些特定基因以期产生超级人类。

现在,他们要做的第一步是将两种不同生物的基因序列转换成两个整数序列,并试图确定他们的最大公共上升子序列的长度,例如有A序列为4 3 2 1 7 8 9,B序列为7 8 9 4 3 2 1,其最长公共子序列是4 3 2 1,而最长公共递增子序列应该是 7 8 9。

输入
输入每个序列由M个整数组成(1 ≤ M ≤500),M个整数范围在(-231 ≤Bi < 231)之间。
输出
输出最长公共上升子序列长度L。
样例输入 Copy
5
1 4 2 5 -12
4
-12 1 2 4
样例输出 Copy
2

解析:
设f[i][j]:表示所有再A[1…i]和B[1…j]都出现过,且以B[j]结尾的公共上升子序列集合。

分为两种集合:
第一种A[i]不包含在最长公共上升子序列中 。状态方程就是:f[i-1][j]

第二种A[i]包含在公共子序列中, 需要在、再划分子序列。这时候A[i]和B[j]相等枚举i-1的状态

#include<bits/stdc++.h>
using namespace std;
const int N=1000;
int dp[N][N];
int a[N];
int b[N];
int n,m;
int main()
{
  scanf("%d",&n);
  for(int i=1;i<=n;i++) cin>>a[i];
  scanf("%d",&m);
  for(int i=1;i<=m;i++) cin>>b[i];
  for(int i=1;i<=n;i++)
  {
    int maxv=1;
    for(int j=1;j<=m;j++)
    {
      dp[i][j]=dp[i-1][j];
      if(a[i]==b[j]) dp[i][j]=max(dp[i][j],maxv);
      if(b[j]<a[i]) maxv=max(maxv,dp[i][j]+1);
    }
  }
    int res=0;
    for(int i=1;i<=m;i++) res=max(dp[n][i],res);
    cout<<res<<endl;
  

}

  • 0
    点赞
  • 2
    收藏
    觉得还不错? 一键收藏
  • 0
    评论
评论
添加红包

请填写红包祝福语或标题

红包个数最小为10个

红包金额最低5元

当前余额3.43前往充值 >
需支付:10.00
成就一亿技术人!
领取后你会自动成为博主和红包主的粉丝 规则
hope_wisdom
发出的红包
实付
使用余额支付
点击重新获取
扫码支付
钱包余额 0

抵扣说明:

1.余额是钱包充值的虚拟货币,按照1:1的比例进行支付金额的抵扣。
2.余额无法直接购买下载,可以购买VIP、付费专栏及课程。

余额充值