poj 1080 Human Gene Functions dp

这个题目也没什么可以说的了

dp[i][j]表示第一段长度为i，第二段长度为j的时候的最优值

状态转移方程

dp[i][j]=max(dp[i][j],dp[i-1][j-1]+mat[a[i]][b[j]]);

dp[i][j]=max(dp[i][j],dp[i-1][j]+mat[a[i]][0]);

dp[i][j]=max(dp[i][j],dp[i][j-1]+mat[0][b[j];

注意边界就ok了

#include

#include

#include

using namespace std;

const int maxn=111,inf=99999999;

char a[maxn],b[maxn];

int mat[501][501];

int dp[maxn][maxn];

const char A='A',G='G',T='T',C='C';

int max(int a,int b)

{

    if(a>b)return(a);

    return(b);

}

int main()

{

    int tcase;

   scanf("%d",&tcase);

   mat[A][A]=mat[C][C]=mat[T][T]=mat[G][G]=5;

   mat[A][C]=mat[C][A]=-1;

   mat[A][T]=mat[T][A]=-1;

   mat[A][G]=mat[G][A]=-2;

   mat[A][0]=mat[0][A]=-3;

   mat[C][G]=mat[G][C]=-3;

   mat[C][T]=mat[T][C]=-2;

   mat[C][0]=mat[0][C]=-4;

   mat[G][T]=mat[T][G]=-2;

   mat[G][0]=mat[0][G]=-2;

   mat[T][0]=mat[0][T]=-1;

    while(tcase--)

    {

       memset(a,0,sizeof(a));

       memset(b,0,sizeof(b));

       int n,m;

       scanf("%d %s %d%s",&n,&a[1],&m,&b[1]);

       dp[0][0]=0;

       for(int i=1;i<=n;i++)

       dp[i][0]=mat[a[i]][0]+dp[i-1][0];

       for(int i=1;i<=m;i++)

       dp[0][i]=mat[b[i]][0]+dp[0][i-1];

       dp[1][1]=mat[a[1]][b[1]];

       for(int i=1;i<=n;i++)

       for(int j=1;j<=m;j++)

       {

          dp[i][j]=-inf;

          dp[i][j]=max(dp[i][j],dp[i-1][j-1]+mat[a[i]][b[j]]);

          dp[i][j]=max(dp[i][j],dp[i-1][j]+mat[a[i]][0]);

          dp[i][j]=max(dp[i][j],dp[i][j-1]+mat[0][b[j]]);

       }

       printf("%d\n",dp[n][m]);

    }

    return 0;

}