P4059【动态规划】[Code+#1]找爸爸

最新推荐文章于 2022-11-17 15:49:21 发布

原创最新推荐文章于 2022-11-17 15:49:21 发布 · 176 阅读

1 ·

CC 4.0 BY-SA版权

动态规划专栏收录该内容

9 篇文章

订阅专栏

本文记录了一位开发者在解决一道动态规划问题时，由于一个小的下标错误导致卡壳，最终通过调整代码并完成蓝题的经历。重点讨论了代码优化和DP算法的应用。

有个位置，写错了一个下标，卡了半个小时。人裂开了
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;
}