『状压DP』[PA2014]Pakowanie

$\mathrm{P}\mathrm{r}\mathrm{o}\mathrm{b}\mathrm{l}\mathrm{e}\mathrm{m}$

你有n个物品和m个包。物品有重量，且不可被分割；包也有各自的容量。

要把所有物品装入包中，至少需要几个包？

$\mathrm{S}\mathrm{o}\mathrm{l}\mathrm{u}\mathrm{t}\mathrm{i}\mathrm{o}\mathrm{n}$

显然，先装大背包再装小背包一定是最优策略。因此我们需要对背包的体积排序。

而我们所需要安排的是背包里具体装的是哪些物品，我们考虑状压DP。

设$f[s]$表示状态$s$中的物品最少需要几个背包，$g[s]$表示放完$s$以后最后的背包所剩下的空间。

那么判断根据剩下体积能否装下，即$c[f[i]]-g[i]$和$a[j]$进行比较，可以得到两种转移。

• 若能够装下。
• 放一个新背包。
  

$\mathrm{C}\mathrm{o}\mathrm{d}\mathrm{e}$

#include <cstdio>
#include <cstring>
#include <iostream>
#include <algorithm>

using namespace std;
const int N = 24;

int n, m;
int c[N*N], a[N+9], f[(1<<N)+9], g[(1<<N)+9];

int read(void)
{
	int s = 0, w = 0; char c = getchar();
	while (c < '0' || c > '9') w |= c == '-', c = getchar();
	while (c >= '0' && c <= '9') s = s*10+c-48, c = getchar();
	return w ? -s : s;
}

void DP(void)
{
	f[0] = g[0] = 0;
	for (int i=0;i<(1<<n)-1;++i)
	{
		if (f[i] > m) continue;
		for (int j=1;j<=n;++j) if (!((i >> j-1) & 1))
		{
			if (c[f[i]] - g[i] >= a[j]) 
			{
				if (f[i] < f[i|(1<<j-1)]) f[i|(1<<j-1)] = f[i], g[i|(1<<j-1)] = g[i] + a[j];
				if (f[i] == f[i|(1<<j-1)]) g[i|(1<<j-1)] = min(g[i|(1<<j-1)],g[i]+a[j]);
			}
			if (c[f[i]] - g[i] < a[j])
			{
				if (f[i] >= m) continue;
				if (c[f[i]+1] < a[j]) continue;
				if (f[i] + 1 < f[i|(1<<j-1)]) f[i|(1<<j-1)] = f[i] + 1, g[i|(1<<j-1)] = a[j];
				if (f[i] + 1 == f[i|(1<<j-1)]) g[i|(1<<j-1)] = min(g[i|(1<<j-1)],a[j]);
			}
		}
	}
	if (f[(1<<n)-1] > m) puts("NIE");
	else printf("%d\n", f[(1<<n)-1]);
	return;
}

int main(void)
{
	freopen("pako.in","r",stdin);
	freopen("pako.out","w",stdout);
	n = read(), m = read();
	for (int i=1;i<=n;++i) a[i] = read();
	for (int i=1;i<=m;++i) c[i] = read();
	sort(c+1,c+m+1);
	reverse(c+1,c+m+1);
	memset(g,30,sizeof g);
	memset(f,30,sizeof f);
	DP();
	return 0;
}