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

类似最长公共的DP
$f_{i,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;
}