poj 1080 zoj 1027(最长公共子序列变种)

http://poj.org/problem?id=1080

http://acm.zju.edu.cn/onlinejudge/showProblem.do?problemId=27

/*
zoj 1027   poj 1080 思路： 
三种状态,取最大值：
s1[i]和s2[j]配 ：dp[i-1][j-1]+cost[my[s1[i]]][my[s2[j]]];
s1[i]和'-' 配： dp[i-1][j]+cost[my[s1[i]]][my['-']];
s2[j]和'-' 配： dp[i][j-1]+cost[my['-']][my[s2[j]]];
注意边界：
d[i][0]= cost[my[s1[i]]][my['-']]+dp[i-1][0]; 只能全部和'-'配，
同理： dp[0][j]=cost[my['-']][my[s2[j]]]+dp[0][j-1];
*/
#include <iostream>
#include <map>
#include <cstring>
#include <string>
using namespace std;
int cost[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,-100000005}};
int dp[1005][1005];
string s1,s2;
map<char,int> my;
int main(int argc, char *argv[])
{
	int t,i,n,j,m;
	my['A']=0; my['C']=1; my['G']=2; my['T']=3; my['-']=4;
	cin>>t;
	while(t--)
	{
		cin>>n>>s1>>m>>s2;
		s1='0'+s1; s2='0'+s2;
		dp[0][0]=0;
		for(i=1;i<=n;i++)
			dp[i][0]=cost[my[s1[i]]][my['-']]+dp[i-1][0];
		for(j=1;j<=m;j++)
			dp[0][j]=cost[my['-']][my[s2[j]]]+dp[0][j-1];
		for(i=1;i<=n;i++)
		for(j=1;j<=m;j++)
		{
			int t1,t2,t3;
			t1=dp[i-1][j-1]+cost[my[s1[i]]][my[s2[j]]];
			t2=dp[i-1][j]+cost[my[s1[i]]][my['-']];
			t3=dp[i][j-1]+cost[my['-']][my[s2[j]]];
			t1=max(t1,t2);  t1=max(t1,t3);
			dp[i][j]=t1;	
		}
		cout<<dp[n][m]<<endl;
	}
	return 0;
}