poj1080

实在是太菜太菜，被这道题折腾了不少时间，后面参看了不少解题报告才勉强算是懂了。

lcs算法的关键部分就是递推（动态规划）和边界条件。

#include<iostream>

#include<fstream>
using namespace std;
int max(int a,int b,int c)
{
return (a>b?a:b)>c?(a>b?a:b):c;
}
const int N=105;
char a[N];
char b[N];
int matrix[N][N];
int r[5][5]={
0,-3,-4,-2,-1,
-3,5,-1,-2,-1,
-4,-1,5,-3,-2,
-2,-2,-3,5,-2,
-1,-1,-2,-2,5
};
int match(char a)
{
switch(a)
{
case '-':return 0;
case 'A':return 1;
case 'C':return 2;
case 'G':return 3;
case 'T':return 4;
}
//return -1;
}
int main()
{
//fstream cin("C:\\Users\\wuyanyisb\\Desktop\\1.txt");
int cases;
cin>>cases;
while(cases--)
{
int m,n;
cin>>m;
cin>>a;
cin>>n;
cin>>b;
matrix[0][0]=0;
for(int i=1;i<=m;i++)
matrix[i][0]=r[match(a[i-1])][0]+matrix[i-1][0];
for(int i=1;i<=n;i++)
matrix[0][i]=r[0][match(b[i-1])]+matrix[0][i-1];
for(int i=1;i<=m;i++)
for(int j=1;j<=n;j++)
{
matrix[i][j]=max(
matrix[i-1][j-1]+r[match(a[i-1])][match(b[j-1])],
matrix[i-1][j]+r[match(a[i-1])][0],
matrix[i][j-1]+r[0][match(b[j-1])]);
}
cout<<matrix[m][n]<<endl;
}
system("pause");
return 0;
}