题解
很明显这是一道动态规划的问题(dp)
思路:首先dp数组的定义,我是这样定义的:
dp[i][j]
代表的是,第一个碱基序列的前i位非空格脱氧核苷酸与第二个碱基的前j位脱氧核苷酸的最大值是多少(就是非空格)
例如
7 AGTGATG
5 GTTAG
dp[2][3]代表的就是AG与GTT相对应的最大值
当i=1,j=1时也就是A与G对应,很显然有三种状态:
很显然最大值为第一种状态,接着这个状态可以更新的状态分别有三个:
i=2,j=1 ; i=1,j=2 ; i=2,j=2
也就是dp[i][j]可以转移到dp[i+1][j] , dp[i][j+1] , dp[i+1][j+1]。
所以同理可以说明dp[i][j]可以由dp[i-1][j] , dp[i][j-1] , dp[i-1][j-1]转移过来
for(int i=1;i<=la;i++)
{
for(int j=1;j<=lb;j++)
{
dp[i][j]=max(dp[i][j],dp[i-1][j]+v[a[i]][5]);//第二个碱基序列的核苷酸匹配空格
dp[i][j]=max(dp[i][j],dp[i][j-1]+v[5][b[j]]);//第一个碱基序列的核苷酸匹配空格
dp[i][j]=max(dp[i][j],dp[i-1][j-1]+v[a[i]][b[j]]);//第一个碱基序列和第二个碱基序列的核苷酸匹配
}
}
有个细节要注意:
2. 要把dp[0][j]和dp[i][0]预处理
因为当一个碱基序列为0时,就是相当于另一个碱基序列全部和空格对应
for(int i=1;i<=la;i++)
{
dp[i][0]=dp[i-1][0]+v[a[i]][5];
}
for(int i=1;i<=lb;i++)
{
dp[0][i]=dp[0][i-1]+v[5][b[i]];
}
最后直接上代码吧
#include <bits/stdc++.h>
using namespace std;
const int N=100010;
int a[N],d1[N],d2[N],n;
int main()
{
while(cin>>a[++n]);n--; //输入
int len1=1,len2=1; //初始长度为1
d1[1]=a[1]; //用于求不上升序列长度
d2[1]=a[1]; //用于求上升序列长度
for(int i=2; i<=n; i++) { //从a[2]开始枚举每个数(a[1]已经加进去了)
if(d1[len1]>=a[i])
d1[++len1]=a[i];
else
{
int p1=upper_bound(d1+1,d1+1+len1,a[i],greater<int>())-d1;
d1[p1]=a[i];
}
if(d2[len2]<a[i])d2[++len2]=a[i]; //同上
else
{
int p2=lower_bound(d2+1,d2+1+len2,a[i])-d2;
d2[p2]=a[i];
}
}
cout<<len1<<endl<<len2; //输出
return 0;
}