ZOJ 1027

最新推荐文章于 2017-07-12 16:01:22 发布

xnature

最新推荐文章于 2017-07-12 16:01:22 发布

阅读量639

点赞数

CC 4.0 BY-SA版权

分类专栏： C++ 文章标签： c

本文链接：https://blog.youkuaiyun.com/xnature/article/details/2788181

C++ 专栏收录该内容

1 篇文章

订阅专栏

本文介绍了一种基于最长公共子序列（LCS）的加权版本算法实现，该算法用于寻找两个序列间的最大相似度得分。通过使用特定的权重矩阵来评估不同字符组合的价值，实现了对传统LCS算法的有效扩展。

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

又经典的DP题LCS改编

原来的LCS是求最长公共子序列

这里把长度换成了权重

#include <cstdlib> 
#include <iostream> 
#include <cstdio> 
using namespace std;
const int MAXN = 101;
const int MAX = 99999999;
const int point[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, 0}};
int f[MAXN][MAXN], data[2][MAXN], num[2];
int main()
{
    int cases;
    cin>>cases;
    while (cases--)
    {
        int i,j;
        for (i = 0; i < 2; i ++)
        {
            scanf("%d",num+i);getchar();
            char ch;
            for (j = 0; j < num[i]; j++)
            {
                ch = getchar();
                switch (ch)
                {
                case 'A':
                    data[i][j] = 0;
                    break;
                case 'C':
                    data[i][j] = 1;
                    break;
        
                case 'G':
                    data[i][j] = 2;
                    break;
               
                case 'T':
                    data[i][j] = 3;
                    break;
               
                }
            }
        }
            int u,v;
            for (u = 0; u < MAXN; u++ )
                for (v = 0; v < MAXN; v++)
                    f[u][v] = -MAX;
            f[0][0] = 0;
            for (u = 0; u <= num[0]; ++u)
            {
                for (v = 0; v <= num[1]; ++v)
                {
                    if (f[u+1][v+1] < f[u][v] + point[data[0][u]][data[1][v]])
                        f[u+1][v+1] = f[u][v] + point[data[0][u]][data[1][v]];
                    if (f[u][v+1] < f[u][v] + point[4][data[1][v]])
                        f[u][v+1] = f[u][v] + point[4][data[1][v]];
                    if (f[u+1][v] < f[u][v] + point[data[0][u]][4])
                        f[u+1][v] = f[u][v] + point[data[0][u]][4];
                }
            }
            cout<<f[num[0]][num[1]]<<endl;
        
    
    }
            
    return 0;
}