【bzoj1087】互不侵犯king 状压dp

最新推荐文章于 2021-03-08 23:18:34 发布

Mininda

最新推荐文章于 2021-03-08 23:18:34 发布

阅读量372

点赞数

CC 4.0 BY-SA版权

分类专栏： bzoj 动态规划-状压dp 文章标签： bzoj 状压dp

本文链接：https://blog.youkuaiyun.com/Dadatu_Zhao/article/details/81139968

bzoj 同时被 2 个专栏收录

34 篇文章

订阅专栏

动态规划-状压dp

4 篇文章

订阅专栏

本文介绍了一道经典的状压动态规划入门题目，并提供了完整的代码实现。通过该例题，读者可以了解到如何利用位运算进行状态压缩，以及如何在状态转移过程中进行有效的剪枝。

摘要生成于 C知道，由 DeepSeek-R1 满血版支持，前往体验 >

状压dp入门题！太经典了就不多bb了。

#include<iostream>
#include<cstdio>
#include<cstring>
#define N 601
#define ll long long
using namespace std;
ll dp[85][N][85],ans;
int f[N][N],a[N],b[N],n,m,ToT;
bool judge(int x) {return !(x&(x<<1));}
bool check(int x,int y) {return !( (x&y) || (x&(y<<1)) || (x&(y>>1)) );}
int cal(int x) {
    int ret=0;
    for (;x;x>>=1)
        if (x&1) ret++;
    return ret; 
}
int main() {
    scanf("%d%d",&n,&m);
    for (int i=1;i<(1<<n);i++)
        if (judge(i)) a[++ToT]=i,b[ToT]=cal(i);
    for (int i=0;i<=ToT;i++)
        for (int j=0;j<=ToT;j++) 
            if (check(a[i],a[j])) f[i][j]=1;
    dp[0][0][0]=1;
    for (int i=1;i<=n;i++)
        for (int j=0;j<=ToT;j++)
            for (int k=0;k<=ToT;k++)
                if (f[k][j]) 
                    for  (int h=(b[k]+b[j]);h<=m;h++)
                        dp[i][j][h]+=dp[i-1][k][h-b[j]];
    for (int i=0;i<=ToT;i++) ans+=dp[n][i][m];
    printf("%lld\n",ans);   
    return 0;
}