hdu 1080 Human Gene Functions(DP)
题意:
人类基因由A、C、G、T组成。
有一张5*5的基因表。每格有一个值,叫相似度。例:A-C:-3。意思是如果A和C配对, 则它俩的相似度是-3【P.S.:-和-没有相似度,即-和-不能配对】
现在给两条基因片段。(长度不一定相等)
现在你要在两条基因片段中插入若干个-(空白基因),使得两个基因片段长度相等,且得到的整体相似度的值最大。【再次P.S.:-和-不能配对】
思路:
因为-和-不能匹配,所以插入的-的个数是有限的。
str1的第一个基因可以与str1的第一个或-配对。然后,,,,很快就看到了DP结构,,,
和最长公共子串一样滴的意思,,,
DP方程:dp[i][j]:str1的前i个和str2的前j个能获得的最大相似度值。
看代码,,
代码:
int l1,l2; char s1[105], s2[105]; int S1[105], S2[105]; int a[6][6]={{5,-1,-2,-1,-3},{-1,5,-3,-2,-4},{-2,-3,5,-2,-2},{-1,-2,-2,5,-1},{-3,-4,-2,-1,inf}}; int dp[105][105]; int maxx(int a,int b,int c){ int t=max(a,b); return max(t,c); } int main(){ int T; cin>>T; while(T--){ scanf("%d",&l1); scanf("%s",s1); scanf("%d",&l2); scanf("%s",s2); rep(i,0,l1-1){ if(s1[i]=='A') S1[i]=0; else if(s1[i]=='C') S1[i]=1; else if(s1[i]=='G') S1[i]=2; else S1[i]=3; } rep(i,0,l2-1){ if(s2[i]=='A') S2[i]=0; else if(s2[i]=='C') S2[i]=1; else if(s2[i]=='G') S2[i]=2; else S2[i]=3; } dp[0][0]=0; rep(i,0,l2-1){ dp[0][i+1]=dp[0][i]+a[4][S2[i]]; } rep(i,0,l1-1){ dp[i+1][0]=dp[i][0]+a[S1[i]][4]; } rep(i,1,l1){ rep(j,1,l2){ dp[i][j]=maxx( dp[i-1][j-1]+a[S1[i-1]][S2[j-1]],dp[i-1][j]+a[S1[i-1]][4],dp[i][j-1]+a[4][S2[j-1]] ); } } printf("%d\n",dp[l1][l2]); } return 0; }