动态规划(acm pku 1080 Human Gene Functions)

http://162.105.81.212/JudgeOnline/problem?id=1080

 

题目意思是说有两段基因序列，而每两个基因都有一个相似值，在任意匹配之下，求出最大相似值。

相似值表格如下，其中'-'代表无基因

DP的做法，类似于最长公共子序列，不过这里不是判定相同字符，而是最大的相似值。

初始化时注意要先找到超过长度的情况，即对应字符与'-'的相似值

 

 

#include <iostream> using namespace std; int dp[300][300] , len1 , len2; char ch1[300] , ch2[300]; int getMax( int a , int b , int c ) { return max( a,max(b,c) ); } int value( char a , char b ) { if ( a=='A' ) { if ( b=='A' ) return 5; if ( b=='C' ) return -1; if ( b=='G' ) return -2; if ( b=='T' ) return -1; if ( b=='-' ) return -3; } if ( a=='C' ) { if ( b=='A' ) return -1; if ( b=='C' ) return 5; if ( b=='G' ) return -3; if ( b=='T' ) return -2; if ( b=='-' ) return -4; } if ( a=='G' ) { if ( b=='A' ) return -2; if ( b=='C' ) return -3; if ( b=='G' ) return 5; if ( b=='T' ) return -2; if ( b=='-' ) return -2; } if ( a=='T' ) { if ( b=='A' ) return -1; if ( b=='C' ) return -2; if ( b=='G' ) return -2; if ( b=='T' ) return 5; if ( b=='-' ) return -1; } if ( a=='-' ) { if ( b=='A' ) return -3; if ( b=='C' ) return -4; if ( b=='G' ) return -2; if ( b=='T' ) return -1; if ( b=='-' ) return 0; } return 0; } int fun() { int i , j; dp[0][0] = 0; //初始化超过长度的情况 for (i=1 ; i<=len1 ; i++) dp[i][0] = dp[i-1][0] + value( ch1[i-1],'-' ); for (i=1 ; i<=len2 ; i++) dp[0][i] = dp[0][i-1] + value( '-',ch2[i-1] ); for (i=1 ; i<=len1 ; i++) for (j=1 ; j<=len2 ; j++) dp[i][j] = getMax( dp[i][j-1]+value(ch2[j-1],'-'),dp[i-1][j]+value(ch1[i-1],'-'),dp[i-1][j-1]+value(ch1[i-1],ch2[j-1])); return dp[len1][len2]; } int main() { int cases; cin >> cases; while (cases--) { memset(dp , 0 , sizeof dp); cin >> len1 >> ch1; cin >> len2 >> ch2; cout << fun() << endl; } return 0; }