题目链接：http://www.lydsy.com/JudgeOnline/problem.php?id=5107
题意：见原题面。
发现空格匹配空格一定是更劣的，所以一个空格至少匹配了一个字母，总长不超过\(n+m\)
所以就可以 dp 了，见注释。
|
1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 |
#include <cmath> #include <queue> #include <cstdio> #include <iomanip> #include <cstdlib> #include <cstring> #include <iostream> #include <algorithm> #define N 3010 #define ll long long #define inf 1000000000 using namespace std; inline int read() { int x=0,f=1;char ch=getchar(); while(ch<'0'||ch>'9'){if(ch=='-')f=-1;ch=getchar();} while(ch>='0'&&ch<='9'){x=x*10+ch-'0';ch=getchar();} return x*f; } char stra[N],strb[N]; int a[N],b[N]; int dp[N][N][3];//A考虑到i B考虑到j 目前为 (字母 字母) (空格 字母) (字母 空格)的最大相似程度 int buc[100]; int d[4][4],A,B,n,m; void pre() { scanf("%s",stra+1),scanf("%s",strb+1); buc[(int)'A']=0,buc[(int)'T']=1,buc[(int)'G']=2,buc[(int)'C']=3; for(int i=0;i<=3;i++) for(int j=0;j<=3;j++)d[i][j]=read(); A=-read(),B=-read(); n=strlen(stra+1),m=strlen(strb+1); for(int i=1;i<=n;i++)a[i]=buc[(int)stra[i]]; for(int i=1;i<=m;i++)b[i]=buc[(int)strb[i]]; } int main() { pre(); for(int i=0;i<=n;i++) for(int j=0;j<=m;j++) for(int k=0;k<=2;k++) dp[i][j][k]=-inf; dp[0][0][0]=0; for(int i=0;i<=n;i++) for(int j=0;j<=m;j++) { if(i&&j)dp[i][j][0]=max(dp[i-1][j-1][0],max(dp[i-1][j-1][1],dp[i-1][j-1][2]))+d[a[i]][b[j]]; if(i) dp[i][j][1]=max(dp[i-1][j][1]+B,max(dp[i-1][j][0],dp[i-1][j][2])+A); if(j) dp[i][j][2]=max(dp[i][j-1][2]+B,max(dp[i][j-1][0],dp[i][j-1][1])+A); } printf("%d\n",max(dp[n][m][0],max(dp[n][m][1],dp[n][m][2]))); } |
Comments | NOTHING