bzoj3717 [PA2014]Pakowanie 状压dp

Description

---

你有n个物品和m个包。物品有重量，且不可被分割；包也有各自的容量。要把所有物品装入包中，至少需要几个包？

第一行两个整数n,m(1<=n<=24,1<=m<=100)，表示物品和包的数量。
第二行有n个整数a[1],a[2],…,a[n](1<=a[i]<=10^8)，分别表示物品的重量。
第三行有m个整数c[1],c[2],…,c[m](1<=c[i]<=10^8)，分别表示包的容量。

如果能够装下，输出一个整数表示最少使用包的数目。若不能全部装下，则输出NIE。

Solution

---

写了一个随机化过不去

考虑状压。我们设f[x]表示选取物品状态为x时最少包的数量，g[x]表示选取物品状态为x，答案为f[x]时最后一个包剩余的最大空间
考虑转移。我们枚举x中1的位置，然后把这个物品塞进包里面就可以了

Code

---

#include <stdio.h>
#include <string.h>
#include <algorithm>
#include <math.h>
#define rep(i,st,ed) for (register int i=st;i<=ed;++i)
#define drp(i,st,ed) for (register int i=st;i>=ed;--i)
#define fill(x,t) memset(x,t,sizeof(x))
#define lowbit(x) ((x)&(-(x)))

const int INF=10777216;
const int N=16778216;

int a[N],b[205],f[N],g[N],n,m;

int read() {
	int x=0,v=1; char ch=getchar();
	for (;ch<'0'||ch>'9';v=(ch=='-')?(-1):(v),ch=getchar());
	for (;ch<='9'&&ch>='0';x=x*10+ch-'0',ch=getchar());
	return x*v;
}

bool cmp(int x,int y) {
	return x>y;
}

int main(void) {
	// freopen("data.in","r",stdin);
	// freopen("myp.out","w",stdout);
	n=read(),m=read(); int lim=(1<<n)-1;
	rep(i,1,n) a[i]=read();
	drp(i,n,1) a[1<<(i-1)]=a[i];
	rep(i,1,m) b[i]=read();
	std:: sort(b+1,b+m+1,cmp);
	rep(i,1,lim) {
		f[i]=INF;
		for (int x=i;x;x-=lowbit(x)) {
			int tmp=lowbit(x);
			if (g[i-tmp]>=a[tmp]&&(f[i-tmp]<f[i]||(f[i-tmp]==f[i]&&g[i-tmp]-a[tmp]>g[i]))) {
				f[i]=f[i-tmp];
				g[i]=g[i-tmp]-a[tmp];
			} else if (b[f[i-tmp]+1]>=a[tmp]&&(f[i-tmp]+1<f[i]||(f[i-tmp]+1==f[i]&&b[f[i-tmp]+1]-a[tmp]>g[i]))) {
				f[i]=f[i-tmp]+1;
				g[i]=b[f[i]]-a[tmp];
			}
		}
	}
	if (f[lim]>m) puts("NIE");
	else printf("%d\n", f[lim]);
	return 0;
}