首页 > 其他 > 详细

【一天一DP计划】状压DP

时间:2019-10-08 13:09:07      阅读:72      评论:0      收藏:0      [点我收藏+]

##P1896 互不侵犯【状压dp】

  • 用01串表示每一行的可行解
  • 列与列之间的限制条件在转移的时候continue
  • x&=(-x)可以找到x的二进制位上有多少个1
  • 关于位运算的优先级!不确定就无脑加括号哦哦
reference:
    https://www.luogu.org/blog/virus2017/1896dnqec
Date:
    2019.10.04
*/
#include<bits/stdc++.h>
using namespace std;
#define int long long
#define rep(i,a,b) for(int i=a;i<=b;++i)
#define dwn(i,a,b) for(int i=a;i>=b;--i) 
template <typename T> inline void rd(T &x){x=0;char c=getchar();int f=0;while(!isdigit(c)){f|=c=='-';c=getchar();}while(isdigit(c)){x=(x<<1)+(x<<3)+(c^48);c=getchar();}x=f?-x:x;} 
#define mem(a,b) memset(a,b,sizeof(a))

int f[9+5][512+5][81+5];
//f[i][j][k]表示第i行的状态选的是第j个状态,有k个国王的方案数
//上个状态:f[i-1][jj][k-cnt[j]]; 
int n,K,tot;
int st[512+5],cnt[512+5];
//st表示当前可行的状态(01串 
//cnt表示st[]的状态有多少个1(为了后面利用背包的思想转移) 

inline void get_num(int x){
    int now=0;
    while(x){
        x&=(x-1);
        now++;
    }
    cnt[tot]=now;
}

inline void pre_work(){
    for(int i=0;i<(1<<n);++i)
        if(!(i&(i<<1)))//遇到位运算最好加括号,比如不能写成i&(-i)==0(==优先级比&高) 
            st[++tot]=i,get_num(i);
}

#undef int
int main(){
#define int long long
    #ifdef WIN64 
    freopen("qinfan.txt","r",stdin);
    #endif 
    rd(n),rd(K);
    pre_work();
    rep(i,1,tot)f[1][i][cnt[i]]=1;
    rep(i,2,n)
        rep(j,1,tot)
            rep(jj,1,tot){
                if(st[j]&st[jj] || (st[j]<<1)&st[jj] || (st[j]>>1)&st[jj])continue;
                dwn(k,K,cnt[j])
                    f[i][j][k]+=f[i-1][jj][k-cnt[j]];
            } 
    int ans=0;
    rep(i,1,tot)
        ans+=f[n][i][K];
    printf("%lld\n",ans);
    return 0;
}

【一天一DP计划】状压DP

原文:https://www.cnblogs.com/sjsjsj-minus-Si/p/11634651.html

(0)
(0)
   
举报
评论 一句话评论(0
关于我们 - 联系我们 - 留言反馈 - 联系我们:wmxa8@hotmail.com
© 2014 bubuko.com 版权所有
打开技术之扣,分享程序人生!