码迷,mamicode.com
首页 > 其他好文 > 详细

LuoguP4059 找爸爸-序列DP

时间:2018-10-22 23:16:42      阅读:196      评论:0      收藏:0      [点我收藏+]

标签:str   简单   int   吐槽   reg   就是   name   line   class   

容我吐槽一波题目名字。。

首先我们考虑如果没有空格连续产生代价这一额外条件,那么是不是很好做。

因为我们注意到,不管怎么放空格,我们最终都是把两个串匹配完了。

同时,同一位置上下都放空格是完全没有意义的。

所以我们直接设f[i][j]表示第一个匹配到了i,第二个匹配到了j。

那么转移就是f[i][j]=max {f[i][j-1], f[j][i-1], f[i-1][j-1]+d[i][j]}。

然而这里多了一个条件。。那该怎么做呢??

同样的,我们仍然容易发现,同一位置上下都放空格不会是结果更优。

那么对于我们就可以把这个求额外贡献的式子拆成:

第一个0就减a,以后的连起来的0都减b。

在此基础上,我们考虑设f[i][j][0/1/2]表示第一个匹配到了i,第二个匹配到了j,且此时末尾(没有/第一个串/第二个串)出现了0。

那么的话,转移就很简单了:

       f[i][j][0]=max (f[i-1][j-1][0], max (f[i-1][j-1][1], f[i-1][j-1][2]))+d[ad[i]][bd[j]];
            f[i][j][1]=max (f[i][j-1][0], f[i][j-1][2])-a;
            f[i][j][1]=max (f[i][j][1], f[i][j-1][1]-b);
            f[i][j][2]=max (f[i-1][j][0], f[i-1][j][1])-a;
            f[i][j][2]=max (f[i][j][2], f[i-1][j][2]-b);

最后稍微注意下边界就好了。。

#include <bits/stdc++.h>
using namespace std;
inline int gi () {
    int x=0,w=0; char ch=0;
    while (!(ch>=0 && ch<=9)) {if (ch==-) w=1; ch=getchar ();}
    while (ch>=0 && ch<=9) x= (x<<3)+(x<<1)+(ch^48), ch=getchar ();
    return w?-x:x;
}

#define RG register

const int N=3010;

char ADNA[N],BDNA[N];
int a,b,alen,blen,ad[N],bd[N],d[5][5],f[N][N][3];

int main ()
{
    RG int i,j;
    scanf ("%s%s", ADNA+1, BDNA+1);
    for (i=1;i<=4;++i)
        for (j=1;j<=4;++j)
            d[i][j]=gi ();
    a=gi (), b=gi ();
    alen=strlen (ADNA+1), blen=strlen (BDNA+1);
    for (i=1;i<=alen;++i) 
        switch (ADNA[i]) {
            case A: ad[i]=1; break;
            case T: ad[i]=2; break;
            case G: ad[i]=3; break;
            case C: ad[i]=4; break;
        }
    for (i=1;i<=blen;++i) 
         switch (BDNA[i]) {
            case A: bd[i]=1; break;
            case T: bd[i]=2; break;
            case G: bd[i]=3; break;
            case C: bd[i]=4; break;
        }
    memset (f, 0xcf, sizeof (f));
    f[0][0][0]=0, f[0][1][1]=f[1][0][2]=-a;
    for (i=1;i<=alen;++i)
        for (j=1;j<=blen;++j) {
            f[i][j][0]=max (f[i-1][j-1][0], max (f[i-1][j-1][1], f[i-1][j-1][2]))+d[ad[i]][bd[j]];
            f[i][j][1]=max (f[i][j-1][0], f[i][j-1][2])-a;
            f[i][j][1]=max (f[i][j][1], f[i][j-1][1]-b);
            f[i][j][2]=max (f[i-1][j][0], f[i-1][j][1])-a;
            f[i][j][2]=max (f[i][j][2], f[i-1][j][2]-b);
        }
    printf ("%d\n", max (f[alen][blen][0], max (f[alen][blen][1], f[alen][blen][2])));
    return 0;
}

 

LuoguP4059 找爸爸-序列DP

标签:str   简单   int   吐槽   reg   就是   name   line   class   

原文地址:https://www.cnblogs.com/Bhllx/p/9833647.html

(0)
(0)
   
举报
评论 一句话评论(0
登录后才能评论!
© 2014 mamicode.com 版权所有  联系我们:gaon5@hotmail.com
迷上了代码!