有个位置,写错了一个下标,卡了半个小时。人裂开了
here
第一次过DP蓝题 虽然写了一个小时,但过了还是挺开心的。
// god with me
//#pragma GCC optimize(1)
//#pragma GCC optimize(2)
//#pragma GCC optimize(3,"Ofast","inline")
//#include <bits/stdc++.h>
#include <bits/stdc++.h>
#include <cmath>
#include <math.h>
#include<cstdio>
#include<string>
#include<cstring>
#include<iostream>
#include<algorithm>
#include<vector>
#define inf 0x7fffffff
#define ll long long
//#define int long long
//#define double long double
//#define double long long
#define re int
//#define i int i
//#define void inline void
#define eps 1e-8
//#define mod 1e9+7
#define ls(p) p<<1
#define rs(p) p<<1|1
#define pi acos(-1.0)
#define pb push_back
#define mk make_pair
#define P pair < int , int >
// typedef long long s64;
using namespace std;
const int mod=1e6+7;
const int M=45;
const int N=3e3+5;//?????????? 4e8
char A[N],B[N];
int sc[10][10];
int a1,b1;
int dp[N][N][3];
int n,m;
int dis(char x,char y)
{
int a;
if(x=='A')a=1;
if(x=='T')a=2;
if(x=='G')a=3;
if(x=='C')a=4;
int b;
if(y=='A')b=1;
if(y=='T')b=2;
if(y=='G')b=3;
if(y=='C')b=4;
return sc[a][b];
}
void godp()
{
memset(dp,-0x3f,sizeof(dp));
// for(int i=0;i<=max(n,m);i++)for(int j=0;j<=max(n,m);j++)for(int k=0;k<=2;k++)dp[i][j][k]=-1e8;
for(int i=1;i<=max(n,m);i++)dp[i][0][2]=dp[0][i][2]=dp[i][0][1]=dp[0][i][1]=-a1-b1*(i-1);
// for(int i=0;i<=max(n,m);i++)dp[i][0][1]=dp[i][0][0]=dp[0][i][1]=dp[0][i][0]=-1e8;
dp[0][0][0]=0;
for(int i=1;i<=n;i++)
for(int j=1;j<=m;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]))+dis(A[i],B[j]);
dp[i][j][1]=max(dp[i-1][j][0]-a1,max(dp[i-1][j][1]-b1,dp[i-1][j][2]-a1));
dp[i][j][2]=max(dp[i][j-1][0]-a1,max(dp[i][j-1][1]-a1,dp[i][j-1][2]-b1));
}
cout<<max(dp[n][m][0],max(dp[n][m][1],dp[n][m][2]))<<endl;
}
void solve()
{
scanf("%s",A+1);
scanf("%s",B+1);
n=strlen(A+1),m=strlen(B+1);
for(int i=1;i<=4;i++)
for(int j=1;j<=4;j++)scanf("%d",&sc[i][j]);
scanf("%d%d",&a1,&b1);
godp();
}
signed main()
{
// ios::sync_with_stdio(false);
int T=1;
// cin>>T;
for(int index=1;index<=T;index++)
{
// printf("Case %d:\n",index);
solve();
// puts("");
}
return 0;
}