poj_1080

源地址：http://poj.org/problem?id=1080

又是一道dp的题目，做dp没感觉啊。。- -

这是一道字符匹配的dp,感觉和最长公共子序列有点像.

最长公共子序列：

if(s[i]==s2[j]) dp[i][j] = dp[i-1][j-1]+1

else dp[i][j] = Max(dp[i-1][j],dp[i][j-1])

而对于这道题目，字母还可以 ’-‘来匹配。

我们用dp[i][j]来表示s1的前i个字母和s2前j个字母匹配后能得到的最大结果。

而事先，我们得做好初始化工作：

对于dp[i][0]来说，dp[i][0] = dp[i-1][0] + score[s1[i]]['-'],因为相当于前i个字母都和'-'匹配，那么匹配后的分数就是题中所给的那个表格中的。

同理，对于dp[0][j]来说，dp[0][j] = dp[0][j-1]+score['-'][s2[j]]。

而对于i>=1&&j>=1来说，就有三种状态：

要么s1[i]和s2[j]匹配，要么s1[i]和'-'匹配，要么'-'和s2[j]匹配，然后比较这三种情况下能得到的分数大小，取大的那一个，最后输出dp[len1][len2]就是s1前len1个字符和s2前len2个字符匹配后的最大结果。

#include<stdio.h>
#include<iostream>
#include<string>
#include<string.h>
#include<algorithm>
#include<iomanip>
#include<vector>
#include<time.h>
#include<queue>
#include<stack>
#include<iterator>
#include<math.h>
#include<stdlib.h>
#include<limits.h>
#include<set>
#include<map>
//#define ONLINE_JUDGE
#define eps 1e-8
#define INF 0x7fffffff
#define FOR(i,a) for((i)=0;i<(a);(i)++)
#define MEM(a) (memset((a),0,sizeof(a)))
#define sfs(a) scanf("%s",a)
#define sf(a) scanf("%d",&a)
#define sfI(a) scanf("%I64d",&a)
#define pf(a) printf("%d\n",a)
#define pfI(a) printf("%I64d\n",a)
#define pfs(a) printf("%s\n",a)
#define sfd(a,b) scanf("%d%d",&a,&b)
#define sft(a,b,c)scanf("%d%d%d",&a,&b,&c)
#define for1(i,a,b) for(int i=(a);i<b;i++)
#define for2(i,a,b) for(int i=(a);i<=b;i++)
#define for3(i,a,b)for(int i=(b);i>=a;i--)
#define MEM1(a) memset(a,0,sizeof(a))
#define MEM2(a) memset(a,-1,sizeof(a))
const double PI=acos(-1.0);
template<class T> T gcd(T a,T b){return b?gcd(b,a%b):a;}
template<class T> T lcm(T a,T b){return a/gcd(a,b)*b;}
template<class T> inline T Min(T a,T b){return a<b?a:b;}
template<class T> inline T Max(T a,T b){return a>b?a:b;}
using namespace std;
#define ll __int64
int n,m;
#define Mod 1000000007
#define N 510
#define M 1000100
int dp[N][N];	//dp[i][j]代表s1的前i个字符和s2的前j个字符之间匹配的最大值
int sc[10][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 getPos(char c){
	if(c == 'A') return 0;
	else if(c == 'C') return 1;
	else if(c == 'G') return 2;
	else if(c == 'T') return 3;
	else if(c == '-') return 4;
}
int main(){
#ifndef ONLINE_JUDGE
    freopen("in.txt","r",stdin);
//  freopen("out.txt","w",stdout);
#endif
    int t;
    sf(t);
    int len1,len2;
    char s1[N],s2[N];
    while(t--){
    	scanf("%d %s",&len1,s1+1);
    	scanf("%d %s",&len2,s2+1);
    	memset(dp,0,sizeof dp);
    	dp[0][0] = 0;
    	for(int i=1;i<=len1;i++){	//边界条件：即s1的前i个字母都和'-'进行匹配
    		dp[i][0] = dp[i-1][0]+sc[getPos(s1[i])][getPos('-')];
    	}
    	for(int i=1;i<=len2;i++){	//边界条件：即s2的前i个字母都和'-'匹配
    		dp[0][i] = dp[0][i-1]+sc[getPos('-')][getPos(s2[i])];
    	}
    	for(int i=1;i<=len1;i++){
    		for(int j=1;j<=len2;j++){
    			dp[i][j] = dp[i-1][j-1] + sc[getPos(s1[i])][getPos(s2[j])];<span style="white-space:pre">	</span>//三种情况
    			dp[i][j] = Max(dp[i][j],dp[i-1][j]+sc[getPos(s1[i])][getPos('-')]);
    			dp[i][j] = Max(dp[i][j],dp[i][j-1]+sc[getPos('-')][getPos(s2[j])]);
    		}
    	}
    	printf("%d\n",dp[len1][len2]);
    }
//    while(sf(n)!=EOF){
//    }
return 0;
}