[DP] 「CodePlus 2017 11 月赛」找爸爸

类似最长公共的DP
fi,j,k,s 表示匹配到A的第i位,B的第j位,k表示A的第i是之前是否有0,s表示B的第j位之前是否有0

瞎转移一下

#include <cstdio>
#include <algorithm>
#include <iostream>
#include <cstring>

using namespace std;

const int N=3010;

char a[N],b[N];

int c[10][10],A,B,n,m;
int f[N][N][2][2],g[N][N];

inline int code(char x){
  if(x=='A') return 1;
  if(x=='T') return 2;
  if(x=='G') return 3;
  return 4;
}

inline void fix(int &x,int y){
  if(x<y) x=y;
}

int main(){
  freopen("6250.in","r",stdin);
  freopen("6250.out","w",stdout);
  scanf("%s%s",a+1,b+1); n=strlen(a+1); m=strlen(b+1);
  memset(f,-0x7f>>1,sizeof(f));
  f[0][0][0][0]=0;
  for(int i=1;i<=4;i++)
    for(int j=1;j<=4;j++)
      scanf("%d",&c[i][j]);
  scanf("%d%d",&A,&B);
  for(int i=0;i<=n;i++)
    for(int j=0;j<=m;j++){
      if(i==j && i==0) continue;
      if(i && j){
    int cur=c[code(a[i])][code(b[j])];
    fix(f[i][j][0][0],f[i-1][j-1][0][0]+cur);
    fix(f[i][j][0][0],f[i-1][j-1][0][1]+cur);
    fix(f[i][j][0][0],f[i-1][j-1][1][0]+cur);
      }
      if(i){
    fix(f[i][j][0][1],f[i-1][j][0][0]-A);
    fix(f[i][j][0][1],f[i-1][j][1][0]-A);
    fix(f[i][j][0][1],f[i-1][j][0][1]-B);
      }
      if(j){
    fix(f[i][j][1][0],f[i][j-1][0][0]-A);
    fix(f[i][j][1][0],f[i][j-1][0][1]-A);
    fix(f[i][j][1][0],f[i][j-1][1][0]-B);
      }
    }
  cout<<max(max(f[n][m][0][0],f[n][m][0][1]),f[n][m][1][0])<<endl;
  return 0;
}
评论
添加红包

请填写红包祝福语或标题

红包个数最小为10个

红包金额最低5元

当前余额3.43前往充值 >
需支付:10.00
成就一亿技术人!
领取后你会自动成为博主和红包主的粉丝 规则
hope_wisdom
发出的红包
实付
使用余额支付
点击重新获取
扫码支付
钱包余额 0

抵扣说明:

1.余额是钱包充值的虚拟货币,按照1:1的比例进行支付金额的抵扣。
2.余额无法直接购买下载,可以购买VIP、付费专栏及课程。

余额充值