在遇到某些动规题目时,发现状态太多,无从下手分析具体情况时,也许可以考虑一下状态压缩DP,一种很暴力但很好用的DP。
一、状态压缩的常规方法
一般来说,状态压缩是利用二进制数只有0和1的特点,将当前较为复杂存储状态表示成一个数储存起来。这样一来,繁多的状态就可以表示成一个数了。
比如背包问题,可以把放该物品记为1,不放记为0,那么三个物品如果只放第一个和第二个,就可以记作二进制的011(2)(默认为从二进制右侧开始记录状态),即十进制的3
二进制数位 | 0 | 1 | 2 |
---|---|---|---|
二进制数 | 1 | 1 | 0 |
状态 | √ | √ | × |
(注意顺序)
那么很容易得出,如果有n个物品,那么状态数最多有2n
种(0~2n-1)。
那么,位运算这个坎。。。就必须要过了。。
二、关于位运算
众所周知,电脑储存、计算数据都是利用二进制的,如果直接利用位运算进行计算,程序就能比普通的四则运算快,这也是为什么位运算很受欢迎。
1. 与(&):只有同为1时才是1,其他都是0;
eg:11001(2)& 10111(2)= 10001(2)
2. 或(|):只有同为0时才是0,其他都是1;
eg:11001(2)| 10111(2)= 11111(2)
3. 非(!):! 0 = 1,! 1 = 0;
4. 异或(^):相同为0,不同为1;
eg:11001(2)^ 10111(2)= 1110(2)
5. 左移(x<<n):将二进制的n所有数位向左移动n位(类似于十进制×2)
eg:11001(2)<<1 = 110010(2)
6. 右移(x>>n):将二进制的n所有数位向右移动n位(类似于十进制/2)
eg:11001(2)>>1 = 1100(2)
网上的一张图总结的不错,也贴在下面,不要求强记,大概知道就行。
三、一些例题~
magic_sheep大佬说我还要多做状压dp的题才能熟练掌握,所以下面的题目可能还会不定期更新~~~
1.Hihocoder P1846 ABC (来自vjudge)
问题描述
杂货铺老板一共有N件物品,每件物品具有ABC三种属性中的一种或多种。从杂货铺老板处购得一件物品需要支付相应的代价。
现在你需要计算出如何购买物品,可以使得ABC三种属性都在购买的物品中出现,并且支付的总代价最小。
输入
第一行包含一个整数N。
以下N行,每行包含一个整数C和一个只包含"ABC"的字符串,代表购得该物品的代价和其具有的属性。
对于50%的数据,1 ≤ N ≤ 20
对于100%的数据,1 ≤ N ≤ 1000 1 ≤ C ≤ 100000
输出
一个整数,代表最小的代价。如果无论如何凑不齐ABC三种属性,输出-1。
样例输入
5
10 A
9 BC
11 CA
4 A
5 B
样例输出
13
分析:
本题只有ABC三种物品属性,所以暴力一下多写几个if就可以过了,但用状压dp就会更清晰,代码量也会更少。
将ABC三种物品分别对应二进制数位的0,1,2位,这样用1表示该属性放了,0表示未放,则很容易列出:
f[j|a[i].w]=min(f[j|a[i].w],f[j]+a[i].v);
首先,f[j]表示在第 j 状态下当前的最小花费,a[i].w记录第i个物品的状态,则 j | a[i].w 的运算即表示在f[j]状态放入第i个物品后的最优解(或运算即可以模拟放入的操作)。
其实,这个状态转移方程式是不是和0/1背包很像呢:
f[j+w]=max(f[j+w],f[j]+v);
所以,这道题基本上就是个裸的背包+状压就行啦!
代码如下
#include <bits/stdc++.h>
using namespace std;
int f[8];
struct bag{
int v,w;
}a[1100];
int main(){
int n;
cin>>n;
for (int i=0;i<=7;i++) f[i]=99999999;
f[0]=0;//初始化一定要记得!
for (int m=1;m<=n;m++){
cin>>a[m].v;
string s;
cin>>s;
for (int i=0;i<s.size();i++){//按二进制存
if (s[i]=='A') a[m].w+=1;
if (s[i]=='B') a[m].w+=2;
if (s[i]=='C') a[m].w+=4;
}
for (int j=0;j<=7;j++){//dp过程
f[j|a[m].w]=min(f[j|a[m].w],f[j]+a[m].v);
}
}
if (f[7]!=99999999) cout<<f[7]<<endl;else cout<<-1<<endl;
return 0;
}
2.Hihocoder P1486 物品价值 (来自vjudge)
问题描述
小Hi现在有n个物品,每个物品都有一个价值。并且这n个物品总共有m个不同的属性,每个物品都具有其中若干属性。
小Ho要从中选出若干物品,满足每个属性都正好有奇数个物品拥有,且被选出的物品价值总和最大。你能帮助小Ho完成任务么?
输入
第一行一个数T(<=10),表示数据组数。对于每一组数据:
第一行两个数n,m(1<=n<=1000,m<=10)
接下来每两行描述一件物品。对于每一件物品:
第一行两个数v和s,表示其价值和所含属性数量(v<=100000,s<=m)
第二行s个数,表示该物品拥有的属性编号(1<=编号<=m)
输出
物品价值总和的最大值。
样例输入
1
3 2
2 1
1
2 1
2
5 2
1 2
样例输出
5
分析
这题明显比上一题稍微难点。。。
首先对于状态表示,题目说必须是奇数个,是不是可以想到,用1表示放了奇数个,0表示放了偶数个,那么通过异或操作就可完成放入过程了(可以自己手动推一下为什么异或可以表示放入状态)。
用f[i][j]表示在放入第i个物品,状态为j的时候的最优解,则可以列出和上题类似的状态转移方程:
f[i][j∧a[i].w]=max(f[i][j∧a[i].w],f[i-1][j]+a[i].v);
(具体的和上一题差不太多,理解就行)
代码如下
#include <bits/stdc++.h>
using namespace std;
struct cd{
int v,w;
}a[1010];
int f[2050][2050];
int main(){
int T,n,m;
cin>>T;
while (T--){
cin>>n>>m;
memset(f,-1,sizeof f);
memset(a,0,sizeof a);
f[0][0]=0;//初始化
int s;
for (int i=1;i<=n;i++){
cin>>a[i].v>>s;
for (int j=1;j<=s;j++){
int e;
cin>>e;
a[i].w|=1<<(e-1);//如此读入效率高
}
for (int j=0;j<(1<<m);j++){
if (f[i-1][j]!=-1){
f[i][j^a[i].w]=max(f[i][j^a[i].w],f[i-1][j]+a[i].v);//状态转移方程
f[i][j]=max(f[i][j],f[i-1][j]);
}
}
}
cout<<max(f[n][(1<<m)-1],0)<<endl;
}
return 0;
}
3.洛谷 P1896 互不侵犯
问题描述
在N×N的棋盘里面放K个国王,使他们互不攻击,共有多少种摆放方案。国王能攻击到它上下左右,以及左上左下右上右下八个方向上附近的各一个格子,共8个格子。
输入
只有一行,包含两个数N,K ( 1 <=N <=9, 0 <= K <= N * N)
输出
所得的方案数
样例输入
3 2
样例输出
16
分析
emmmmmmm,是不是无从下手。。。
暴力深搜?
想多了,老弟。。。(你要是能过就算我没说 )
考虑到状压,所以我们可以把每一行放或不放用1和0表示,这样一来,每行放置的状态就可以很轻易地表示出来。
首先我们要先预处理,判断一行中何种状态满足题中条件,这要用到与(&)运算。我们发现,判断二进制中1是否相邻两个数都是0,只需要判断 x&(x<<1) 和 x&(x>>1) 是否都为0即可(需要手推及理解,这里不细说)。摆放棋子的个数只需要计算一下二进制1的个数即可。
我们用f[i][j][k]表示当第i行是第j个状态并且放了k个棋子时的方案数。
那么,第一行的所有状态都可以预处理进数组了。接着,只需要当这行的状态与上一行不冲突,就更新方案数,差不多这题就OK了!
代码如下
#include <bits/stdc++.h>
using namespace std;
long long a[1000],f[10][1<<9][100];
int n,m;
bool g[1000];
int main(){
scanf("%d%d",&n,&m);
memset(f,0,sizeof f);
memset(g,false,sizeof g);
for (int x=0;x<1<<n;x++){
if (!(x&(x<<1))&&!(x&(x>>1))) g[x]=true;
int s=x;
while (s){
if (s&1) a[x]++;
s>>=1;
}
if (g[x]) f[1][x][a[x]]=1;
}
for (int i=2;i<=n;i++){
for (int j=0;j<1<<n;j++){
if (g[j]){
for (int k=0;k<1<<n;k++)
if (g[k]&&!(j&k)&&!(j&(k>>1))&&!(j&(k<<1)))
for (int l=0;l+a[k]<=m;l++)
f[i][k][l+a[k]]+=f[i-1][j][l];
}
}
}
long long ans=0;
for (int i=0;i<(1<<n);i++)
ans+=f[n][i][m];
printf("%lld\n",ans);
return 0;
}