您的位置:首页 > 其它

hdu 1080 Human Gene Functions(类最长公共子序列)

2011-08-14 13:17 309 查看
初始化一定要做好。

不然就如此地悲剧,

刚开始我还是严格按照最长公子序来写的,比较的时候是满足条件再赋值。

谁知道它就是三个一起比较。

准备用Map来映射,一直出错,不得不采用如此蹩脚的代码了。

最后,还想用滚动数组来写,发现……,不能正确初始化。

算了,反正100的长度也不长。

#include<iostream>
using namespace std;

#define MAX_M 101
#define MAX(a,b) (a)>(b)?(a):(b)
#define GetNum(c) G(c)
int value[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,-200} };

int GetNum(char c){
switch(c){
case 'A':return 0;
case 'C':return 1;
case 'G':return 2;
case 'T':return 3;
default:return 4;
}
}

void Init(int dp[][MAX_M],int n,int m,char s[],char t[]){
dp[0][0]=0;
for(int i=1;i<=m;i++)
dp[i][0]=dp[i-1][0]+value[G(s[i-1])][4];
for(int j=1;j<=n;j++)
dp[0][j]=dp[0][j-1]+value[4][G(t[j-1])];
}
int main(){

int T,n,m,dp[MAX_M][MAX_M];
char s[MAX_M],t[MAX_M];
cin>>T;
while(T--){
cin>>n>>s;
cin>>m>>t;
Init(dp,n,m,s,t);
for(int i=1;i<=n;i++)
for(int j=1;j<=m;j++)
dp[i][j]=MAX(MAX(dp[i-1][j]+value[G(s[i-1])][G('-')],
dp[i][j-1]+value[G('-')][G(t[j-1])]),
dp[i-1][j-1]+value[G(s[i-1])][G(t[j-1])]);
cout<<dp
[m]<<endl;
}
return 0;
}
内容来自用户分享和网络整理,不保证内容的准确性,如有侵权内容,可联系管理员处理 点击这里给我发消息
标签: