HDU 1080 Human Gene Functions 带权的最长公共子序列

点击打开链接
题意：给定两个字符串 s1 和 s2 ，在两个串中都可以插入空格，使两个串的长度最后相等，然后开始匹配，怎样插入空格由匹配规则得到的值最大。

分析：LCS的变形

d[i][j]=max(d[i-1][j-1]+match[ s1[i] ][ s2[j] ] , d[i-1][j]+match[ s1[i] ][ '-' ] , d[i][j-1]+match[ '-' ][ s2[j] ] );

边界：d[0][0]=0; d[i][0]=d[i-1][0] + match[ s1[i] ][ '-' ] ; d[0][i]=d[0][i-1] + match[ '-' ][ s2[i] ];

 

d[i][j] = d[i-1][j]+match[ s1[i] ][ '-' ] 表示 s1[i-1] 与 s2[j] 匹配 于是 s1[i] 就只能和 ‘-’（空格）匹配了

#include<bits/stdc++.h>
#define INF 0x3f3f3f3f
#define M 1000
using namespace std;
int dp[M][M];
char s1[M],s2[M];
int vis[1000][1000];
int main()
{
    int T,n,m;
    scanf("%d",&T);
    vis['G']['G']=vis['T']['T']=vis['A']['A']=vis['C']['C']=5;

    vis['-']['A']=vis['A']['-']=-3;
    vis['-']['C']=vis['C']['-']=-4;
    vis['-']['T']=vis['T']['-']=-1;
    vis['-']['G']=vis['G']['-']=-2;

    vis['A']['C']=vis['C']['A']=-1;
    vis['A']['G']=vis['G']['A']=-2;
    vis['A']['T']=vis['T']['A']=-1;

    vis['G']['C']=vis['C']['G']=-3;
    vis['T']['C']=vis['C']['T']=-2;
    vis['T']['G']=vis['G']['T']=-2;
    while(T--)
    {
        scanf("%d%s%d%s",&n,s1+1,&m,s2+1);
        memset(dp,0,sizeof(dp));
        for(int i=1;i<=n;i++)
            dp[i][0]=dp[i-1][0]+vis[s1[i]]['-'];   ///初始化，为了使开头可以有‘-’与字符匹配的情况
        for(int i=1;i<=m;i++)
            dp[0][i]=dp[0][i-1]+vis['-'][s2[i]];

        for(int i=1;i<=n;i++)
        {
            for(int j=1;j<=m;j++)
            {                
                dp[i][j]=max(dp[i-1][j-1]+vis[s1[i]][s2[j]],max(dp[i-1][j]+vis[s1[i]]['-'],dp[i][j-1]+vis['-'][s2[j]]));
            }
        }
        printf("%d\n",dp[n][m]);
    }
    return 0;
}