bzoj 4565 字符合并 DP 解题报告

最新推荐文章于 2020-12-15 21:27:54 发布

原创最新推荐文章于 2020-12-15 21:27:54 发布 · 666 阅读

1 ·

本内容遵循CC 4.0 BY-SA版权协议

标签

#bzoj

————单个题目——— 同时被 3 个专栏收录

222 篇文章

订阅专栏

————DP————

50 篇文章

订阅专栏

DP——普通

20 篇文章

订阅专栏

本文介绍了一个关于01串合并的算法问题，通过状压动态规划方法求解最大得分。输入包括字符串长度、合并长度及具体规则，输出为最大可能得分。

Description

有一个长度为 n 的 01 串，你可以每次将相邻的 k 个字符合并，得到一个新的字符并获得一定分数。得到的新字
符和分数由这 k 个字符确定。你需要求出你能获得的最大分数。

Input

第一行两个整数n，k。接下来一行长度为n的01串，表示初始串。接下来2k行，每行一个字符ci和一个整数wi，ci
表示长度为k的01串连成二进制后按从小到大顺序得到的第i种合并方案得到的新字符,wi表示对应的第i种方案对应
获得的分数。1<=n<=300,0<=ci<=1,wi>=1,k<=8

Output

输出一个整数表示答案

Sample Input

3 2

101

1 10

0 20

1 30

Sample Output

//第3行到第6行表示长度为2的4种01串合并方案。00->1,得10分，01->1得10分，10->0得20分，11->1得30分。

思路

i,j,k都比较小，可以考虑状压dp。
令f[i][j][S]表示将i~j这一段消到S这个状态能获得的最大得分。
只考虑两种转移：
①类似区间dp的方法，枚举端点mid，枚举i~mid的状态，mid+1~j的状态为0或1。注意这里每次消除只能是k的整数倍，所以mid+1~j只能是k的整数倍。
即 $f[i][j][S<<1]=max(f[i][j][S<<1],f[i][m−1][S]+f[m][j][0])或f[i][j][S<<1|1]=max(f[i][j][S<<1|1],f[i][m−1][S]+f[m][j][1])$
刚开始没有这样写，而是枚举整个区间的状态之后，用S>>1和S&1转移。不过这样的话有一个问题，就是消除区间的时候由于是二进制表示，001,01,和1都表示成了一个状态没法区分。
②将长度正好可以消除的S直接向c[S]转移。但是这里要注意一个问题，就是要用一个临时数组记录一下，就是刚转移完的状态不能紧接着转移下一个状态，这里也存在我上面说的那个状态表示的问题。转移方程： $g[c[S]]=max(g[c[S]],f[i][j][S]+w[S])，f[i][j][0]=g[0],f[i][j][1]=g[1]。$
时间复杂度O(n3∗2k)，但是实际上有效的转移不会很多。注意转移来的状态必须是合法的状态，可以刚开始把f数组赋成负无穷来区分。

代码

#include<iostream>  
#include<cstdio>  
#include<cstring>  
#include<algorithm>  
#include<cmath>
#include<vector>
using namespace std;  
const int N=200+5;
const int inf=0x3f3f3f3f;
int n,k,c[N],a[N],dp[N][N][N],w[N];  
char s[N];   
int main()  
{  
    freopen("merge.in","r",stdin);
    freopen("merge.out","w",stdout);
    scanf("%d%d",&n,&k);
    scanf("%s",s+1);
    for (int i=1;i<=n;i++)
    a[i]=s[i]-'0';
    for (int i=0;i<=(1<<k)-1;i++)  
    {  
        char s1[2]; 
        scanf("%s",s1+1);  
        c[i]=s1[1]-'0';
        scanf("%d",&w[i]);
    }
    memset(dp,-inf,sizeof(dp));
    for (int i=n;i>=1;i--)
    for (int j=i;j<=n;j++)//枚举区间 
    {
        if (i==j) {dp[i][j][a[i]]=0;continue;}
        int len=j-i;
        while(len>=k) len-=k-1;
        for (int q=j;q>i;q-=k-1)//枚举断点 
        for (int t=(1<<len)-1;t>=0;t--)
        {
            int now=dp[i][q-1][t];
            if (now==-inf) continue;//如果上一次的都没有算过就不算了 
            if (dp[q][j][0]!=-inf)
            dp[i][j][t<<1]=max(dp[i][j][t<<1],dp[q][j][0]+now);
            if (dp[q][j][1]!=-inf)
            dp[i][j][t<<1|1]=max(dp[i][j][t<<1|1],dp[q][j][1]+now);
        }
        if (len==k-1)
        {
            for (int t=(1<<k)-1;t>=0;t--)
            if (dp[i][j][t]!=-inf) 
            dp[i][j][t]=max(dp[i][j][t],dp[i][j][t]+w[t]);
        }  
    }  
    int ans=0;  
    for (int i=(1<<k)-1;i>=0;i--)  
    ans=max(ans,dp[1][n][i]);  
    printf("%d\n",ans);
    return 0; 
}