动态规划---状压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 ...
随机推荐
- 控制台——EventLog实现事件日志操作
我们应该如何通过写代码的方式向其中添加“日志”呢? 在操作之前,先明确几个概念: 1:事件日志名(logName):“事件查看器”中的每一项,如“应用程序”.“Internet Explorer”.“ ...
- inline-block兼容IE7
{ display:inline-block; *display:inblock; *zoom:1 }
- 为什么阿里规约手册要求谨慎使用Arrays.asList方法
前言 在开发中,有时候会碰到把多个参数,或者说把数组转成List的需求,通常我们会使用 Arrays.asList()方法.但是该方法在使用的过程中,稍有不慎就会出现严重的异常.有如下代码: @Tes ...
- EF入门
1.(安装EF)右键项目
- How To:配置Linux iSCSI客户端
1.安装客户端 [root@node01 Packages]# rpm -Uvh iscsi-initiator-utils-6.2.0.873-2.el6.x86_64.rpm warning: i ...
- [JS]window.location获取url各项参数详解
window.location方法后还还可以带href,search等参数,下面我们来看看获取url各项参数的办法. URL即:统一资源定位符 (Uniform Resource Locator, U ...
- Luogu P1041 [2003NOIP提高组]传染病控制
P1041 传染病控制 题目背景 近来,一种新的传染病肆虐全球.蓬莱国也发现了零星感染者,为防止该病在蓬莱国大范围流行,该国政府决定不惜一切代价控制传染病的蔓延.不幸的是,由于人们尚未完全认识这种传染 ...
- (C/C++学习)6.数组指针和指针数组
说明:int (*p)[4] 和 int *p[4](数组指针和指针数组),如果你是一个初学者,也许当你看到这两个名词的时候,已经懵了.其实,只要你理解了其中的含义.这两个名词对你来说会相当简单并且很 ...
- 23.match_phrase_prefix实现search-time搜索推荐
主要知识点: 搜索推荐的使用场景 用法 原理 一.搜索推荐的使用场景 搜索推荐,就是在你做搜索时,当你写出一部搜索词时,es会自提示接下来要写的词,比如当你在搜索hello w 时,如果es中有如下文 ...
- Python-----基本操作
Python是一种简单易学,且功能强大的编程语言.它是面向对象的编程语言. 对象的意义: 对象可以通过一个“.” 的方式来调用这个对象的方法. Python环境安装配置: Python是一种通用的计算 ...