动态规划---状压dp
状压dp,就是把动态规划之中的一个个状态用二进制表示,主要运用位运算。
这里有一道例题:蓝书P639猛兽军团1 [SCOI2005]互不侵犯
题目:
题目描述 在N×N的棋盘里面放K个国王,使他们互不攻击,共有多少种摆放方案。国王能攻击到它上下左右,以及左上左下右上右下八个方向上附近的各一个格子,共8个格子。 注:数据有加强(//)
输入输出格式
输入格式: 只有一行,包含两个数N,K ( <=N <=, <= K <= N * N) 输出格式: 所得的方案数 输入输出样例
输入样例#: 复制 输出样例#: 复制
直接上代码,注释很详细
#include<cstdio>
#include<iostream>
#include<cstring>
#define N 15
#define M 110
#define MAX 550
using namespace std;
/*
见蓝书641页
*/
int s[MAX]; // 记录一行可能的状态
int num[MAX]; //s数组对应每个状态放了多少个猛兽
int states;
long long f[N][M][MAX]; //f[i][j][k]第i行状态为k,放了j个猛兽
int n,m;
void init_state() //预处理s,num数组,代表一行之内所有的可能性
{
states = ;
for(int i = ; i < ( << n); i++) //注意,这里是枚举状态
{
if(i & (i << )) //处理一排上的冲突情况
continue;
int t = i;
num[states] = ;
while(t)
{
num[states] += (t & );
t = t >> ;
}
s[states++] = i; //保存状态
}
}
void dp()
{
int a,c,mm,b,cc;
long long ans;
memset(f,,sizeof(f));
//单独算第一行和最后一行
for(int i = ; i < states; i++)
{
int j = num[i];
if(j <= m) //不能超过总数
f[][j][i]++;
}
for(int i = ; i < n; i++) //2~n - 1行
{
for(int j = ; j <= m; j++) // 到第i行,一共放了j个猛兽
{
for(a = ; a < states; a++) //i行状态
{
c = num[a];
if(c > j)
continue;
mm = j - c;//前i - 1行的总数
for(int b = ; b < states; b++) //枚举i-1行
{
cc = num[b];
if(cc > mm)
continue;
if(s[a] & s[b]) //上下有攻击
continue;
if(s[a] & (s[b] << )) //对角有攻击
continue;
if(s[b] & s[a] << ) // 同上
continue;
f[i][j][a] += f[i - ][mm][b];
}
}
}
}
ans = ;
for(a = ; a < states; a++) //最后一行
{
c = num[a];
if(c > m)
continue;
int j = m - c;
for(int b = ; b < states; b++) //枚举n-1行
{
cc = num[b];
if(cc > j)
continue;
if(s[a] & s[b]) //上下有攻击
continue;
if(s[a] & (s[b] << )) //对角有攻击
continue;
if(s[b] & s[a] << ) // 同上
continue;
f[n][m][a] += f[n - ][j][b];
}
ans += f[n][m][a];
}
printf("%lld\n",ans);
}
int main()
{
scanf("%d%d",&n,&m);
init_state();
dp();
return ;
}
上面这个代码过于复杂,不好理解,我们换一种写法:
#include<bits/stdc++.h>
using namespace std;
const int M=<<;
long long g[M],h[M],f[][M][],n,k,tot=;
int main()
{
cin>>n>>k;
memset(f,,sizeof(f));
for(int x=; x < (<<n); x++) //第一排单独处理
{
if(!(x & (x>>)) && !(x&(x<<)))g[x] = ;//处理g数组,同一排左右不矛盾
int w = x;
while(w)
{
if(w % )h[x]++;
w /= ;
}
if(g[x])
f[][x][h[x]] = ;
}
for(int x = ; x <= n; x++)
{
for(int y = ; y < (<<n); y++)
{
if(g[y])
{
for(int z = ; z < (<<n); z++)
{
if(g[z] && !(y&z) && !(y & (z>>)) && !(y&(z<<))) //只用考虑该排与上一排
{
for(int w = ; w + h[z] <= k; w++)
f[x][z][w + h[z]] += f[x - ][y][w]; //w枚举总共放的国王的个数
}
}
}
}
}
for(int y=; y<(<<n); y++)tot+=f[n][y][k];
cout<<tot;
return ;
}
动态规划---状压dp的更多相关文章
- 状态压缩动态规划 状压DP
总述 状态压缩动态规划,就是我们俗称的状压DP,是利用计算机二进制的性质来描述状态的一种DP方式 很多棋盘问题都运用到了状压,同时,状压也很经常和BFS及DP连用,例题里会给出介绍 有了状态,DP就比 ...
- 状态压缩动态规划(状压DP)详解
0 引子 不要999,也不要888,只要288,只要288,状压DP带回家.你买不了上当,买不了欺骗.它可以当搜索,也可以卡常数,还可以装B,方式多样,随心搭配,自由多变,一定符合你的口味! 在计算机 ...
- 动态规划专题(一)——状压DP
前言 最近,决定好好恶补一下我最不擅长的\(DP\). 动态规划的种类还是很多的,我就从 状压\(DP\) 开始讲起吧. 简介 状压\(DP\)应该是一个比较玄学的东西. 由于它的时间复杂度是指数级的 ...
- 动态规划晋级——POJ 3254 Corn Fields【状压DP】
转载请注明出处:http://blog.csdn.net/a1dark 分析:刚开始学状压DP比较困难.多看看就发现其实也没有想象中那么难.这道题由于列数较小.所以将行压缩成二进制来看.首先处理第一行 ...
- HDU5117 Fluorescent 期望 计数 状压dp 动态规划
原文链接https://www.cnblogs.com/zhouzhendong/p/HDU5117.html 题目传送门 - HDU5117 题意 $T$ 组数据. 给你 $n$ 盏灯 ,$m$ 个 ...
- 动态规划:状压DP
状压DP可以用在NP问题的小规模求解中(不理解,感觉和可以搜索的题很类似) 如果状态是个网格,数据范围很小,基本锁定状压DP 例题是BZOJ1725 题意是这样的,给定一个黑白图,然后种田,要求田与田 ...
- [bzoj1879][Sdoi2009]Bill的挑战_动态规划_状压dp
Bill的挑战 bzoj-1879 Sdoi-2009 题目大意: 注释:$1\le t \le 5$,$1\le m \le 15$,$1\le length \le 50$. 想法: 又是一个看数 ...
- [bzoj3717][PA2014]Pakowanie_动态规划_状压dp
Pakowanie bzoj-3717 PA-2014 题目大意:给你n个物品m个包,物品有体积包有容量,问装下这些物品最少用几个包. 注释:$1\le n\le 24$,$1\le m\le 100 ...
- 【状压DP】bzoj1087 互不侵犯king
一.题目 Description 在N×N的棋盘里面放K个国王,使他们互不攻击,共有多少种摆放方案.国王能攻击到它上.下.左.右,以及左上.左下.右上.右下八个方向上附近的各一个格子,共8个格子. I ...
随机推荐
- Struts2框架实现简单的用户登入
Struts框架汲取了Struts的优点,以WebWork为核心,拦截器,可变和可重用的标签. 第一步:加载Struts2 类库: 第二步:配置web.xml <?xml version=&qu ...
- 【技术累积】【点】【java】【25】Orderd
基础概念 Orderd是spring core中定义的一个接口,使用它以及相关的Comparator和@Order注解,可以实现对元素的排序. @Order 直接先说下@Order注解吧,使用场景较多 ...
- chinason工作室-兄弟的工作室开张了,欢迎来访喔!
Chinason工作室,团队成员由多位多年从事软件开发及大型生产企业系统维护的工程师组成,借重传统国外协同软件的开发经验,结合国内企业实际需求,致力于本土企业工作流软件研发,workflow系统定制开 ...
- Centos6.7 编译安装 Apache PHP
Centos6.7 编译安装 Apache PHP 原文地址:http://www.cnblogs.com/caoguo/p/4968039.html ##### Apache 编译安装 #### [ ...
- =new、=null、.clear()、system.gc()的区别
开发经验告诉我 = new是指向另一个地址空间 =null对象被回收 .clear()对象被清空,但是仍然指向原来的地址空间 这三种方式都并没有真正的清理内存 只有system.gc()是直接清理,但 ...
- UpLoadify在IE下兼容问题
一.在IE9.IE10不能点击的问题解决 解决方法:进入uploadify的js文件中,搜索SWFUpload.prototype.getFlashHTML,找到它对应的语句,将方法全部替换为以下内容 ...
- 通过JS唤醒app(安卓+ios)
有需求说要通过页面按钮唤醒app,或者手机上没有这款app跳转到商店,然后刚开始也是查了资料的,结果发现一头雾水,不过最后还是捣鼓出来了,当然也参考了前人分享的经验,下面我就将方法整理一下: 首先明确 ...
- 数据结构与算法(6) -- heap
binary heap就是一种complete binary tree(完全二叉树).也就是说,整棵binary tree除了最底层的叶节点之外,都是满的.而最底层的叶节点由左至右又不得有空隙. 以上 ...
- Java Web项目实战第1篇之环境搭建
写在前面的话 从今天开始一个Java Web实战项目,参考自 http://blog.csdn.net/eson_15/article/details/51277324 这个博客(非常感谢博主的分享精 ...
- Accessoft-日期区间段查询示例,开始日期至截止日期区段查询
Accessoft-日期区间段查询示例,开始日期至截止日期区段查询 实现功能效果如下: 示例查询开始日期为2017年3月15日到2017年3月16日的内容: sql查询语句如下: SELECT Inf ...