BZOJ 4565 状压DP

思路：

f[i][j][S]表示从i到j压成S状态

j-m是k-1的倍数

$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]);$

为了区分001 01 和1

更新的时候要新开一个数组记录 最后再更新

//By SiriusRen
#include <cstdio>
#include <cstring>
#include <algorithm>
using namespace std;
#define int long long
int n,k,f[][][],g[],ans,c[],w[];
char s[];
signed main(){
    scanf("%lld%lld",&n,&k);
    scanf("%s",s+);
    for(int i=;i<(<<k);i++)scanf("%lld%lld",&c[i],&w[i]);
    memset(f,0xcf,sizeof(f));
    for(int i=;i<=n;i++)s[i]-='',f[i][i][s[i]]=;
    for(int l=;l<=n;l++)
        for(int i=;i<=n-l+;i++){
            int j=i+l-,len=j-i;
            while(len>=k)len-=k-;
            for(int m=j;m>i;m-=k-)
                for(int S=;S<(<<len);S++)
                    f[i][j][S<<]=max(f[i][j][S<<],f[i][m-][S]+f[m][j][]),
                    f[i][j][S<<|]=max(f[i][j][S<<|],f[i][m-][S]+f[m][j][]);
            if(len==k-){
                g[]=g[]=-100000000000000ll;
                for(int S=;S<(<<k);S++)
                    g[c[S]]=max(g[c[S]],f[i][j][S]+w[S]);
                f[i][j][]=g[],f[i][j][]=g[];
            }
        }
    for(int i=;i<(<<k);i++)ans=max(ans,f[][n][i]);
    printf("%lld\n",ans);
}