POJ 1080 Gene

本题收获：
1 利用字符ascii值来替代数组下标，使代码更清晰
2 注意dp数组初始化时不是统一初始化为0，是依次错开所有位置（与-对应的情况）而且对于这些情况需要累加

#include <iostream>
#include <cstdio>
#include <algorithm>
#include <cstring>
using namespace std;#define MAX_L 100
char ca[MAX_L+10], cb[MAX_L+10];
string a, b;
int dp[MAX_L+MAX_L+10][MAX_L+MAX_L+10], score[MAX_L+100][MAX_L+100];//ascii z = 90 enoughint T;int main()
{score['A']['A'] = 5;score['A']['C'] = -1;score['A']['G'] = -2;score['A']['T'] = -1;score['A']['-'] = -3; //注意这种技巧score['C']['A'] = -1;score['C']['C'] = 5;score['C']['G'] = -3;score['C']['T'] = -2;score['C']['-'] = -4;score['G']['A'] = -2;score['G']['C'] = -3;score['G']['G'] = 5;score['G']['T'] = -2;score['G']['-'] = -2;score['T']['A'] = -1;score['T']['C'] = -2;score['T']['G'] = -2;score['T']['T'] = 5;score['T']['-'] = -1;score['-']['A'] = -3;score['-']['C'] = -4;score['-']['T'] = -2;score['-']['G'] = -1;scanf("%d", &T);while(T){T--;int lena, lenb;cin >> lena >> a;/*for(int i = 1;i <= lena;i++)ca[i] = a[i-1];*/cin >> lenb >> b;/*for(int i = 1;i <= lenb;i++)cb[i] = b[i-1];*/memset(dp, 0, sizeof(dp));//initializationfor(int i = 1;i <= lena;i++)dp[i][0] = dp[i-1][0] + score[a[i-1]]['-'];for(int j = 1;j <= lenb;j++)dp[0][j] = dp[0][j-1] + score['-'][b[j-1]];//极限情况就是完全错开for(int i = 1;i <= lena;i++)for(int j = 1;j <= lenb;j++){if(dp[i-1][j-1]+score[a[i-1]][b[j-1]] > dp[i][j-1]+score[a[i-1]]['-']){if(dp[i-1][j-1]+score[a[i-1]][b[j-1]] > dp[i-1][j]+score['-'][b[j-1]])dp[i][j] = dp[i-1][j-1]+score[a[i-1]][b[j-1]];elsedp[i][j] = dp[i-1][j]+score['-'][b[j-1]];}else{if(dp[i][j-1]+score[a[i-1]]['-'] < dp[i-1][j]+score['-'][b[j-1]])dp[i][j] =dp[i-1][j]+score['-'][b[j-1]];elsedp[i][j] = dp[i][j-1]+score[a[i-1]]['-'];}//dp[i][j] = max(dp[i-1][j-1]+score[a[i-1]][b[j-1]], max(dp[i-1][j]+score['-'][b[j-1]], dp[i][j-1]+score[a[i-1]]['-']));//printf("%d ", dp[i][j]);}printf("\n");printf("%d\n", dp[lena][lenb]);}return 0;
}

POJ 1080 Gene相关推荐

POJ 1080 Human Gene Functions(DP:LCS)
题意: 人类基因有两条,现在给出两个字符串,有长有短,可以在字符串中插入空格,要求对照表使两个字符串得出的相似值最大. 要点: 是LCS的变形,用dp[i][j]存储第一个字符串前i个字符与第二个字符 ...
POJ - 1080 Human Gene Functions解题报告（求最长相同子序列）
题目大意: 太长的英语,我就直接读了测试样例和题目中一些可以看的....好像是关于什么人类基因工程什么的, 腺嘌呤脱氧核糖核苷酸~ 每组测试实例给你两串由AGCT四个字母组成的字符串(每个字符串100 ...
poj 1080 Human Gene Functions (最长公共子序列变形)
题意:有两个代表基因序列的字符串s1和s2,在两个基因序列中通过添加"-"来使得两个序列等长:其中每对基因匹配时会形成题中图片所示匹配值,求所能得到的总的最大匹配值. 题解:这题运 ...
动态规划总结与题目分类
源博客链接:http://blog.csdn.net/cc_again/article/details/25866971 动态规划一直是ACM竞赛中的重点,同时又是难点,因为该算法时间效率高,代码量少 ...
『ACM-算法-动态规划』初识DP动态规划算法
一.多阶段决策过程的最优化问题在现实生活中,有类活动的过程,由于它的特殊性,可将过程分成若干个互相阶段.在它的每一阶段都需要作出决策,从而使整个过程达到最好的活动效果.当阶段决策的选取不是任意确 ...
（转）dp动态规划分类详解
dp动态规划分类详解转自:http://blog.csdn.NET/cc_again/article/details/25866971 动态规划一直是ACM竞赛中的重点,同时又是难点,因为该算法时间 ...
《动态规划》— 动态规划分类
动态规划(英语:Dynamic programming,DP)是一种在数学.计算机科学和经济学中使用的,通过把原问题分解为相对简单的子问题的方式求解复杂问题的方法. 动态规划常常适用于有重叠子问题和最 ...
转：动态规划题目分类
https://blog.csdn.net/cc_again/article/details/25866971 一.简单基础dp 这类dp主要是一些状态比较容易表示,转移方程比较好想,问题比较基本常见 ...
【DP专辑】ACM动态规划总结
转载请注明出处,谢谢. http://blog.csdn.net/cc_again?viewmode=list ---------- Accagain 2014年5月15日 ...
ACM比赛经验、刷题记录及模板库总结（更新中）
前言本文所提及的部分题目代码,可以在我的Github上找到第一部分经验分享及感受第二部分刷题记录一.基础算法&程序语言 //strlen()函数的复杂度是O(n)要小心 //截取字 ...

POJ 1080 Gene

POJ 1080 Gene相关推荐

最新文章

热门文章