一道字符串DP,然而不需要状压之类的玄学操作

题目大意:给你两个串,由'A','C','G','T'组成,现在你可以在这两个串中的某些位置插入'-',最终要使得它们的长度相等

给出两个字符匹配时的匹配度表格(题面中给出),让你求如何添加'-'使得总匹配度最大

很简单的DP,我们用f[i][j]表示使用了第一个串中前i个字符,第二个串中前j个字符(注意一定是使用了,即人为填上去的'-'不算)

然后对于每一个位置都有三种选择:

  • 让两个串相应位置上的字符直接匹配,即f[i][j]=max(f[i][j],f[i-1][j-1]+g[s1[i]][s2[j]]

  • 在第一个串中填入一个'-',即f[i][j]=max(f[i][j],f[i][j-1]+g['-'][s2[j]]

  • 在第二个串中填入一个'-',即f[i][j]=max(f[i][j],f[i-1][j]+g[s1[i]]['-']

最后的ans就是f[len1][len2]

CODE

#include<cstdio>
#include<cstring>
using namespace std;
const int N=105,g[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[N<<1][N<<1],t,len1,len2;
char s1[N],s2[N];
inline int get(char ch)
{
if (ch=='A') return 0;
if (ch=='C') return 1;
if (ch=='G') return 2;
if (ch=='T') return 3;
if (ch=='-') return 4;
}
inline int max(int a,int b)
{
return a>b?a:b;
}
int main()
{
register int i,j; scanf("%d",&t);
while (t--)
{
scanf("%d",&len1); scanf("%s",s1+1);
scanf("%d",&len2); scanf("%s",s2+1);
memset(f,167,sizeof(f)); f[0][0]=0;
for (i=1;i<=len1;++i)
f[i][0]=f[i-1][0]+g[get(s1[i])][get('-')];
for (i=1;i<=len2;++i)
f[0][i]=f[0][i-1]+g[get('-')][get(s2[i])];
for (i=1;i<=len1;++i)
for (j=1;j<=len2;++j)
{
f[i][j]=max(f[i][j],f[i-1][j]+g[get(s1[i])][get('-')]);
f[i][j]=max(f[i][j],f[i][j-1]+g[get('-')][get(s2[j])]);
f[i][j]=max(f[i][j],f[i-1][j-1]+g[get(s1[i])][get(s2[j])]);
}
printf("%d\n",f[len1][len2]);
}
return 0;
}

其实这个很显然也可以滚存,但我实在懒得写了

POJ1080的更多相关文章

  1. 【poj1080】 Human Gene Functions

    http://poj.org/problem?id=1080 (题目链接) 题意 给出两个只包含字母ACGT的字符串s1.s2,可以在两个字符串中插入字符“-”,使得s1与s2的相似度最大. Solu ...

  2. LCS poj1080

    题目链接:https://vjudge.net/problem/POJ-1080 参考博客:https://yq.aliyun.com/ziliao/372259 题意:给两个字符串,只含有'A',' ...

  3. POJ-1080 Human Gene Functions---类似LCS

    题目链接: https://cn.vjudge.net/problem/POJ-1080 题目大意: 给定两组序列,要你求出它们的最大相似度,每个字母与其他字母或自身和空格对应都有一个打分,求在这两个 ...

  4. poj1080 dp

    //Accepted 200 KB 0 ms //dp //dp[i][j]表示s1用前i个,s2用前j个字符能得到的最大分数 //dp[i][j]=max(dp[i-1][j]+score[s1[i ...

  5. POJ1080 Human Gene Functions(LCS)

    题目链接. 分析: 和 LCS 差不多. #include <iostream> #include <cstdio> #include <cstdlib> #inc ...

  6. POJ1080 Human Gene Functions 动态规划 LCS的变形

    题意读了半年,唉,给你两串字符,然后长度不同,你能够用'-'把它们补成同样长度,补在哪里取决于得分,它会给你一个得分表,问你最大得分 跟LCS非常像的DP数组 dp[i][j]表示第一个字符串取第i个 ...

  7. POJ1080(LCS变形)

    Human Gene Functions Time Limit: 2000/1000 MS (Java/Others)    Memory Limit: 65536/32768 K (Java/Oth ...

  8. poj1080 - Human Gene Functions (dp)

    题面 It is well known that a human gene can be considered as a sequence, consisting of four nucleotide ...

  9. poj分类 很好很有层次感。

    初期: 一.基本算法:      (1)枚举. (poj1753,poj2965)      (2)贪心(poj1328,poj2109,poj2586)      (3)递归和分治法.      ( ...

随机推荐

  1. C#和C++互相操作 结构体数组的传递

    C++中结构体定义: typedef struct // 平面 { double time;  float normal[3]; float center[3];  } plane; C++中方法声明 ...

  2. OleDbConnection SqlConnection DB2Connection 区别

    OleDbConnection适合于连接任何类型的数据库(如Oracle,SQL   Server,ACCESS等),其命名空间为:using System.Data.OleDb;.而SqlConne ...

  3. 脚本设置IP bat 命令行设置自动获取IP和固定IP

    由于办公室网络需要固定IP和DNS才能上网, 在连接公共网络或者家里又需要自动获取IP和DNS才能上网. 频繁手动切换很麻烦,就搞了两个脚本一键设置. 1.新建文本文件, 命名为固定IP.bat 复制 ...

  4. 网工最实用最常用的网络命令之一——Ping 命令详解(一)

    Ping是Windows.Unix和Linux系统下的一个命令.ping也属于一个通信协议,是TCP/IP协议的一部分.利用“ping”命令可以检查网络是否连通,可以很好地帮助我们分析和判定网络故障. ...

  5. NPOI 导出Excel 数据方式

    使用NPOI的库进行Excel导出操作 公共帮助类: using NPOI.HSSF.UserModel; using NPOI.SS.UserModel; using System; using S ...

  6. 全球首款完全开源的堡垒机,符合 4A 的专业运维审计系统Jumpserver

    Jumpserver是全球首款完全开源的堡垒机,是符合 4A 的专业运维审计系统. http://www.jumpserver.org https://github.com/jumpserver/ju ...

  7. linux正则表达式(基础正则表达式+扩展正则表达式)

    正则表达式应用非常广泛,例如:php,Python,java等,但在linux中最常用的正则表达式的命令就是grep(egrep),sed,awk等,换句话 说linux三剑客要想能工作的更高效,就一 ...

  8. 解决iPhone滑动时滑到另一个层级导致卡顿问题

    问题概览: 两个div都可以滑动时,会造成滑动顶层div时,底层div也会跟着滑动.如图示. 解决方法: 添加CSS即可. 代码如下 * { -webkit-overflow-scrolling: t ...

  9. Netty入门(四)ByteBuf 字节级别的操作

     Netty 中使用 ByteBuf 代替 Java NIO 提供的 ByteBuffer 作为字节的容器. 一.索引 ByteBuf 提供两个指针变量支持读和写操作,读操作使用 readerInde ...

  10. NOIP模拟赛-2018.11.7

    NOIP模拟赛 如果用命令行编译程序可以发现没加头文件之类的错误. 如果用命令行编译程序可以发现没加头文件之类的错误. 如果用命令行编译程序可以发现没加头文件之类的错误. 编译之前另存一份,听说如果敲 ...