P1140 相似基因 题解

题解
很明显这是一道动态规划的问题(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;
}