hdu 1080 ;poj 1080 Human Gene Functions

最新推荐文章于 2019-04-13 00:13:05 发布

原创最新推荐文章于 2019-04-13 00:13:05 发布 · 574 阅读

1 ·

CC 4.0 BY-SA版权

文章标签：

#dp #最长公共子序列 #字符串

本文介绍了一种基于动态规划的字符串匹配算法，该算法通过计算两个字符串间的最大匹配得分来找出最优匹配路径。文章详细解释了如何初始化状态转移矩阵，并通过具体实现代码展示了如何递推地更新状态转移矩阵以获得最终的最大匹配得分。

摘要生成于 C知道，由 DeepSeek-R1 满血版支持，前往体验 >

求两串字符串的匹配。类似于求最长公共子序列。不同的是：上下两串不同的匹配有不同的权值。

可以在两串加入‘-’，开始题目理解有误，以为只能在长度较短的那串加入‘-’ 。说到底还是太水了。

这个题目的初始化需要注意一下

#include <iostream>
#include<algorithm>
#include<cstdio>
#include<cstring>
using namespace std;
int sorce[5][5]= {{5,-1,-2,-1,-3},{-1,5,-3,-2,-4},{-2,-3,5,-2,-2},{-1,-2,-2,5,-1},{-3,-4,-2,-1,-1000}};
char s[5]= {'A','C','G','T','-'};
int get_index(char stc[],int id,int len)
{
    if(id<0||id>=len)return 4;
    int i;
    for(i=0; i<5; i++)
    {
        if(stc[id]==s[i])return i;
    }
    return 4;
}
int dp[111][111];
int main()
{
    int T,n,m,i,j,len1,len2,num2;
    char s1[111],s2[111];
    scanf("%d",&T);
    while(T--)
    {
        scanf("%d %s %d %s",&len1,s1,&len2,s2);
        //初始化 除0 0 外全部初始化为负无穷
        for(i=0; i<=len1; i++)
            for(j=0; j<=len2; j++)
                dp[i][j]=-1111111;
        dp[0][0]=0;
        //dp[i][j] 表示匹配到s1的第i个位置，s2的第j个位置
        for(i=0; i<=len1; i++)
        {
            for(j=0; j<=len2; j++)
            {
                int temx=get_index(s1,i-1,len1);
                int temy=get_index(s2,j-1,len2);
                if(i==0&&j==0)continue;
                //如果 前面加的都是‘-’
                else if(i==0)dp[i][j]=dp[i][j-1]+sorce[4][temy];
                else if(j==0)dp[i][j]=dp[i-1][j]+sorce[temx][4];
                //否则
                else dp[i][j]=max(dp[i-1][j-1]+sorce[temx][temy],max(dp[i-1][j]+sorce[temx][4],dp[i][j-1]+sorce[4][temy]));
            }
        }
        printf("%d\n",dp[len1][len2]);
    }
    return 0;
}