SDUT 1305 查找基因序列问题 dp

题目： http://acm.sdut.edu.cn/sdutoj/problem.php?action=showproblem&problemid=1305

 

这个题就是一个类似公共子串的dp，但是加了权值，其实还是很简单。

当匹配时相似程度是5，不匹配是-4，添加空隙是-7。那么dp[i][j]的值就来自于三个方式：

1、在s2[j]后添加空隙，那么相似程度就是dp[i][j-1] - 7

2、在s1[i]后添加空隙，那么相似程度就是dp[i-1][j] - 7

3、直接拿s1[i]和s2[j]匹配，如果匹配相似程度就是dp[i-1][j-1] + 5，如果不匹配相似程度就是dp[i-1][j-1] - 4

最后找一下字典序最小的，就是这样了。。。就是代码写的搓了一点。

 

#include <stdio.h>
#include <string.h>
#include <algorithm>

int dp[110][110];
char s1[110], s2[1010][110];

int main()
{
    int t, ans = 0, x;
    scanf("%s %d", s1, &t);
    int n = strlen(s1);
    for(int item = 0; item < t; item++)
    {
        scanf("%s", s2[item]);
        int m = strlen(s2[item]);
        dp[0][0] = 0;
        for(int i = 1; i <= n; i++)
            dp[i][0] = dp[i-1][0] - 7;
        for(int j = 1; j <= m; j++)
            dp[0][j] = dp[0][j-1] - 7;
        for(int i = 1; i <= n; i++)
        {
            for(int j = 1; j <= m; j++)
            {
                dp[i][j] = std::max(std::max(dp[i-1][j] - 7, dp[i][j-1] - 7), dp[i-1][j-1] + (s1[i-1] == s2[item][j-1] ? 5 : -4));
            }
        }
        if(dp[n][m] > ans || (dp[n][m] == ans && strcmp(s2[x], s2[item]) > 0))
        {
            ans = dp[n][m];
            x = item;
        }
    }
    printf("%d\n%s\n", ans, s2[x]);
    return 0;
}