POJ 1080

考查：LCS 动态规划

提交情况：1次AC

对状态转移方程又有了一个稍微深入的理解，但还差很远。对LCS的三种情况在这里的运用还没太懂。对于LCS，算法导论中的结论是，设两个字符串X={x0,x1...xm},Y={y0,y1...yn}和它们的最长公共子序列Z={z0,z1...zk},则：

(1)xm=yn，那么可得到zk=xm=yn且zk-1是Xm-1和Yn-1的一个最长公共子序列

(2)xm!=yn，那么由zk=xm可得到Z是Xm和Yn-1的一个最长公共子序列

(3)xm!=yn，那么由zk=yn可得到Z是Xm-1和Yn的一个最长公共子序列

因此若设c[i,j]为最长公共子序列的长度，则可得状态转移方程为

c[i,j]=0, i=0 或j=0

c[i,j]=c[i-1,j-1]+1, xi=yj

c[i,j]=max{c[i-1,j],c[i,j-1]}, xi!=yj

AC_CODE

#include <stdio.h> int table[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 score[105][105]; void ctoi(char *s,int n) { for(int i=0;i<n;i++) { switch(s[i]) { case 'A': s[i]='0'; break; case 'C': s[i]='1'; break; case 'G': s[i]='2'; break; case 'T': s[i]='3'; break; case '-': s[i]='4'; break; } } } int dp(int lena,int lenb,char *s1,char *s2) { ctoi(s1,lena); ctoi(s2,lenb); score[0][0]=0;// for(int i=1;i<=lena;i++) score[i][0]=score[i-1][0]+table[s1[i-1]-'0'][4]; for(int i=1;i<=lenb;i++) score[0][i]=score[0][i-1]+table[4][s2[i-1]-'0']; for(int i=1;i<=lena;i++)//从1开始 <=len for(int j=1;j<=lenb;j++) { score[i][j]=score[i-1][j-1]+table[s1[i-1]-'0'][s2[j-1]-'0']; if(score[i][j]<score[i-1][j]+table[s1[i-1]-'0'][4]) score[i][j]=score[i-1][j]+table[s1[i-1]-'0'][4]; if(score[i][j]<score[i][j-1]+table[4][s2[j-1]-'0']) score[i][j]=score[i][j-1]+table[4][s2[j-1]-'0']; } return score[lena][lenb]; } int main() { int cases,lena,lenb; char s1[105],s2[105]; scanf("%d",&cases); while(cases--) { scanf("%d%s%d%s",&lena,s1,&lenb,s2); printf("%d/n",dp(lena,lenb,s1,s2)); } return 0; }