2287: 【POJ Challenge】消失之物
Time Limit: 10 Sec Memory Limit: 128 MBSubmit: 875 Solved: 499
[ Submit][ Status][ Discuss]
Description
ftiasch 有 N 个物品, 体积分别是 W1, W2, ..., WN。 由于她的疏忽, 第 i 个物品丢失了。 “要使用剩下的 N - 1 物品装满容积为 x 的背包,有几种方法呢?” -- 这是经典的问题了。她把答案记为 Count(i, x) ,想要得到所有1 <= i <= N, 1 <= x <= M的 Count(i, x) 表格。
Input
第1行:两个整数 N (1 ≤ N ≤ 2 × 103) 和 M (1 ≤ M ≤ 2 × 103),物品的数量和最大的容积。
第2行: N 个整数 W1, W2, ..., WN, 物品的体积。
Output
一个 N × M 的矩阵, Count(i, x)的末位数字。
Sample Input
1 1 2
Sample Output
11
21
HINT
如果物品3丢失的话,只有一种方法装满容量是2的背包,即选择物品1和物品2。
Source
HOME Back
新年开的第一道题... 话说csdn又改版了?? 怎么感觉就是像大白皮肤靠拢啊. markdown界面变小不过全屏模式很爽啊, xhtml编辑模式更新了好评.
这道题可以用线段树分治或者背包做.
首先我们考虑线段树分治. 暴力背包怎么做呢? for每个物品的时候把其他物品加进来, 这样是n^2v的, 显然过不了此题. 我们发现有些物品被重复加进去了很多次. 这是因为第i个物品被加进去的区间实际上就是[1, i) 和(i, n]. 那么我们就可以考虑线段树分治, 把for的顺序当做时间, 对时间进行线段树分治. 把每个点所存在的时间映射成两个区间. 在线段树分治的过程中就像线段树那样, 如果某个物品完整包括当前区间就把它加进背包, 否则分类讨论下方或分裂再下放区间即可(就是线段树里面L<=mid, R>mid的操作), 那么到叶子结点我们就把所有覆盖当前点的物品加了进来, 此时的背包也就是题目中关于i所要求的背包辣. 由于每个背包的存在区间最多2个(1, n只有1个), 那么由线段树的性质可以知道最多分成2 logn个区间, 那么每个物品最多被加入背包2logn次, 时间复杂度也就是O(nvlogn). 那么具体实现可以参考po姐的实现. 本蒟蒻就嘴巴Ac了
那么用背包怎么做呢? 我们看到这种限制很简单的就要想到补集转化辣. count(i, x)就是填满体积为j的背包的方案数 - 用了i号物品的填满体积为j的背包的方案数. 那么设c[i][j]表示不用i号物品填满体积为j的背包的方案数, f[j]为填满体积为j的背包的方案数. 那么用了i号物品的填满体积为j的背包的方案数怎么算呢? 对于背包问题我们其实可以在空间上限制来算. c[i][j - w[i]]实际上就表达了这种方案数. 因为我们只需要钦定最后一个物品为i即可. 为什么不是f[j - w[i]]? 因为有可能这种方案里面也用了i. 那么i就被用了两次, 这是不合法的. 当w[i] < j的时候, c[i][j] = f[j], 因为不可能用到i号物品. 那么我们就可以递推完成此题辣.
实现过程中发现枚举i的时候c数组只与当前i有关... 省掉一维XD.
#include<bits/stdc++.h>
using namespace std;
short n, m;
short a[2005], c[2005], f[2005];
int main() {
scanf("%d%d", &n, &m);
register short i, j;
for (i = 1; i <= n; ++ i) scanf("%d", &a[i]);
f[0] = 1;
for (i = 1; i <= n; ++ i)
for (j = m; j >= a[i]; -- j) {
f[j] += f[j - a[i]];
if (f[j] >= 10) f[j] -= 10;
}
for (i = 1; i <= n; ++ i, puts("")) {
for (j = 1; j <= m; ++ j) c[j] = 0;
c[0] = 1;
for (j = 1; j <= m; ++ j) {
if (j >= a[i]) c[j] = f[j] - c[j - a[i]];
else c[j] = f[j];
if (c[j] < 0) c[j] += 10;
}
for (j = 1; j <= m; ++ j) printf("%d", c[j]);
}
return 0;
}