状压DP(超详细!!!)

1、定义

总述

状态压缩动态规划,就是咱们俗称的状压DP,是利用计算机二进制的性质来描述状态的一种DP方式。 算法

不少棋盘问题都运用到了状压,同时,状压也很常常和BFSDP连用。 spa

状压dp其实就是将状态压缩成2进制来保存 其特征就是看起来有点像搜索,每一个格子的状态只有1或0 ,是另外一类很是典型的动态规划blog

举个例子:有一个大小为n*n的农田,咱们能够在任意处种田,如今来描述一下某一行的某种状态: 搜索

设n = 9;遍历

有二进制数 100011011(九位),每一位表示该农田是否被占用,1表示用了,0表示没用,这样一种状态就被咱们表示出来了:见下表二进制

因此咱们最多只须要 2^(n + 1) - 1的十进制数就好(二进制形式是n个1) 如今咱们有了表示状态的方法,但内心也会有些不安:上面用十进制表示二进制的数,枚举了所有的状态,DP起来复杂度岂不是很大?没错,状压实际上是一种很暴力的算法,由于他须要遍历每一个状态,因此将会出现2^n的状况数量,不过这并不表明这种方法不适用:一些题目能够依照题意,排除不合法的方案,使一行的总方案数大大减小从而减小枚举方法

为了更好的理解状压dp,首先介绍位运算相关的知识。 im

  1. ’&’符号,x&y,会将两个十进制数在二进制下进行与运算(都1为1,其他为0) 而后返回其十进制下的值。例如3(11)&2(10)=2(10)。
  2. ’|’符号,x|y,会将两个十进制数在二进制下进行或运算(都0为0,其他为1) 而后返回其十进制下的值。例如3(11)|2(10)=3(11)。
  3. ’^’符号,x^y,会将两个十进制数在二进制下进行异或运算(不一样为1,其他 为0)而后返回其十进制下的值。例如3(11)^2(10)=1(01)。
  4. ’~’符号,~x,按位取反。例如~101=010。
  5. ’<<’符号,左移操做,x<<2,将x在二进制下的每一位向左移动两位,最右边用0填充,x<<2至关于让x乘以4。 ’>>’符号,是右移操做,x>>1至关于给x/2,去掉x二进制下的最右一位

1.判断一个数字x二进制下第i位是否是等于1。(最低第1位)数据

方法:if(((1<<(i−1))&x)>0) 将1左移i-1位,至关于制造了一个只有第i位 上是1,其余位上都是0的二进制数。而后与x作与运算,若是结果>0, 说明x第i位上是1,反之则是0。 img

2.将一个数字x二进制下第i位更改为1。

方法:x=x|(1<<(i−1)) 证实方法与1相似。

3.将一个数字x二进制下第i位更改为0。

方法:x=x&~(1<<(i−1))

4.把一个数字二进制下最靠右的第一个1去掉。

方法:x=x&(x−1)

2、典型例题

【例题1】骑士(P1896 [SCOI2005]互不侵犯)

题目描述

在 n×n(1<=n<=10) 的棋盘上放 k(0<=k<n×n)个国王,国王可攻击相邻的 8 个格子,求使它们没法互相攻击的方案总数。

输入格式

输入有多组方案,每组数据只有一行,包含两个整数 n 和 k。

输出格式

每组数据一行为方案总数,若不可以放置则输出 0。

输入样例

3 2

4 4

样例输出

16

79


实际状压dp顾名思义,就是采用位运算,来记录更多的必须记录的状态来作dp有了比较深的dp功底后只要对位运算有了解就能够解决问题。。。

考虑到每行每列之间都有互相的约束关系。所以,咱们能够用行和列做为另外一个状态的部分。用一个新的方法表示行和列的状态:数字。考虑任何一个十进制数均可以转化成一个二进制数,而一行的状态就能够表示成这样——例如:1010(2)

就表示:这一行的第一个格子没有国王,第二个格子放了国王,第三个格子没有放国王,第四个格子放了国王。而这个二进制下的数就能够转化成十进制: 10(10)

因而,咱们的三个状态就有了:第几行(用i表示)、此行放什么状态(用j表示)、包括这一行已经使用了的国王数(用s表示)。

考虑状态转移方程。咱们预先处理出每个状态(s[x])其中包含二进制下1的个数,及此状态下这一行放的国王个数(num[x]),因而就有:

f[i][j][s]=sum(f[i−1][k][s−num[j]]),f[i][j][s]就表示在只考虑前i行时,在前i行(包括第i行)有且仅有s个国王,且第i行国王的状况是编号为j的状态时状况的总数。而k就表明第i-1行的国王状况的状态编号

 

【例题2】牧场的安排(P1879 [USACO06NOV]玉米田Corn Fields)

Farmer John 新买了一块长方形的牧场,这块牧场被划分红 M列 N 行  (1≤M≤12;1≤N≤12),每一格都是一块正方形的土地。FJ 打算在牧场上的某几格土地里种上美味的草,供他的奶牛们享用。遗憾的是,有些土地至关的贫瘠,不能用来放牧。而且,奶牛们喜欢独占一块草地,因而 FJ 不会选择两块相邻的土地,即:没有哪两块草地有公共边。固然,FJ 尚未决定在哪些土地上种草。 做为一个好奇的农场主,FJ 想知道,若是不考虑草地的总块数,那么,一共有多少种种植方案可供他选择。固然,把新的牧场荒废,不在任何土地上种草,也算一种方案。请你帮 FJ 算一下这个总方案数。

输入格式

第 1行:两个正整数 M 和 N,用空格隔开; 第 2到 M+1行:每行包含 N 个用空格隔开的整数,描述了每块土地的状态。输入的第 i+1行描述了第 i行的土地。全部整数均为 0 或 1,1 表示这块土地足够肥沃,0 则表示这块地上不适合种草。

输出格式

第 1 行:输出一个整数,即牧场分配总方案数除以 10^8的余数。

样例输入

2 3

1 1 1

0 1 0

样例输出

9


 题目大意

给N*M的棋盘,每一个格子不是0就是1,1表明能够种草,不然不能。相邻两个格子不能同时种草,求种草的方案总数。

思路

状态压缩类动态规划,状压dp通常会有明显的数据范围特征,即n,m通常都在20之内。可将每一排的N个当作一个N位二进制,先预处理出每一行能够运行的状态,这样能够去掉不少无效状态(如110),而后DP处理,枚举当前有效状态和上一行有效状态的关系。

f[i][j] 表示第i行在状态j的时候的方案数,其中j咱们用一个二进制数来表示。

转移的时候只要判断与当前行和上一行是否冲突便可,若是不冲突,分f[i][j]=∑f[i−1][k]其中k为不冲突的状态。Ans=∑1≤i≤numf[n][i] 就是最后的答案(num为状态总数)。

初始条件:f[1][i]=1 (1<=i<=a[1].num).

相关文章
相关标签/搜索