动态规划---状压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 ...
随机推荐
- .net 大数据量,查找Where优化(List的Contains与Dictionary的ContainsKey的比较)
最近优化一个where查询条件,查询时间很慢,改为用Dictionary就很快了. 一.样例 假设:listPicsTemp 有100w条数据,pictures有1000w条数据. 使用第1段代码执 ...
- SQL基本操作——存储过程
存储过程类似于C#中的方法. --创建存储过程 create proc usp_TwoNumberAdd @num1 int, @num2 int as begin select @num1+@num ...
- Oracle基础理论笔记(一):模式概念
---oracle 10g 1.在oracle数据库中,数据对象是以模式为单位进行组织和管理的.模式是指一系列的逻辑数据结构或对象的集合. 2.模式与用户名相对应,一个模式只能对应一个用户,并且该模式 ...
- 【技术累积】【点】【java】【28】Map遍历
Map遍历 map的遍历一般有几种吧 遍历entrySet for(Map.Entry<String,String> entry : map.entrySet()){ } Iterator ...
- 【sqli-labs】 less60 GET -Challenge -Double Query -5 queries allowed -Variation3 (GET型 挑战 双查询 只允许5次查询 变化3)
http://192.168.136.128/sqli-labs-master/Less-60/?id=1")%23 http://192.168.136.128/sqli-labs-mas ...
- rrdtool 实践
rrdtool 实践 rrdtool 参数比较多,如果直接看文档,不知从何入手,直接从例子入手这样容易理解,模拟网卡流量 1. 创建数据库 rrdtool create Flow.rrd --star ...
- js 闭包 定时器
; !function (win) { ; //内部私有 , ; //内部私有 //test.prototype.tt1 = 0;//共有变量 var test = function () {}; t ...
- JS数组reduce()方法
1.语法 arr.reduce(callback,[initialValue]) reduce 为数组中的每一个元素依次执行回调函数,不包括数组中被删除或从未被赋值的元素,接受四个参数:初始值(或者上 ...
- linux 的sed命令解释 sed ':t;N;s/\n/,/;b t' 将换行符换成逗号
linux 的sed命令解释 sed ':t;N;s/\n/,/;b t' 将换行符换成逗号 实现的功能是吧换行符换成逗号了,自己试验过. 求解释,:t N b t 都是什么意思??? :t 定义la ...
- HDU 3152 Obstacle Course(优先队列,广搜)
题目 用优先队列优化普通的广搜就可以过了. #include<stdio.h> #include<string.h> #include<algorithm> usi ...