洛谷P1140相似基因

一、题目描述

https://www.luogu.com.cn/problem/P1140

二、思路

这道题感觉还是比较简单。相比刚刚做的那道题，状态是非常好找的，无非就是两个字符串找最长。
首先我们将碱基对的相似匹配使用二维数组存下来，由于输入的是字符，我们可以在输入时直接将其转化为整数。
dp[i][j]：第一个碱基长度为i ，第二个碱基长度为j时 ，最大相似度。
状态转移方程也比较容易想到，比较i和j的相似度时，就迭代来缩小长度比较。
迭代：i-1 比 j-1 或者 i-1 比空 或者 空比j-1 求三者最大。
即

4所对的即表示 为空

dp[i][j]=max(dp[i-1][j-1]+g[a[i]][b[j]],dp[i-1][j]+g[a[i][4]],dp[i][j-1]+g[4][b[j]]

初始化：
dp[0][0]：空对空为0
dp[0][i] ：第0行 即第二个碱基让其每一个都对空，第i个对空就等于前i-1个对空加上 第i个对空的最大相似
dp[i][0]：同理。

三、代码

#include <iostream>
#include <cstdio>
using namespace std;

int g[5][5]={
    //定义基因相似度表格
    //       A    C    G    T  -
    /*A*/ {5,-1,-2,-1,-3},
    /*C*/ {-1,5,-3,-2,-4},
    /*G*/ {-2,-3,5,-2,-2},
    /*T*/ {-1,-2,-2,5,-1},
    /*-*/ {-3,-4,-2,-1,0}
//不能两个空 碰头
};
int main(){
    int m,n,dp[105][105];
    char s;
    int a[105],b[105];
    scanf("%d",&m);
    //随读 随将字符转成表格下标
    for(int i=1;i<=m;i++){
        cin>>s;
        switch(s){
            case 'A':a[i]=0;break;
            case 'C':a[i]=1;break;
            case 'G':a[i]=2;break;
            case 'T':a[i]=3;break;
        }
    }
    scanf("%d",&n);
    for(int i=1;i<=n;i++){
        cin>>s;
        switch(s){
            case 'A':b[i]=0;break;
            case 'C':b[i]=1;break;
            case 'G':b[i]=2;break;
            case 'T':b[i]=3;break;
        }
    }
    //初始化 a的对空；空对b
    dp[0][0]=0;
    for(int i=1;i<=m;i++)
        dp[i][0]=dp[i-1][0]+g[a[i]][4];
    for(int i=1;i<=n;i++)
        dp[0][i]=dp[0][i-1]+g[4][b[i]];
    //dp:
    for(int i=1;i<=m;i++){
        for(int j=1;j<=n;j++){
        //状态转移方程：
            dp[i][j]=max(dp[i-1][j-1]+g[a[i]][b[j]],dp[i-1][j]+g[a[i]][4]);
            dp[i][j]=max(dp[i][j],dp[i][j-1]+g[4][b[j]]);
        }
    }
    cout<<dp[m][n];
    return 0;
}