[BZOJ]2287 消失之物 线段树分治|背包

2287: 【POJ Challenge】消失之物

Time Limit: 10 Sec   Memory Limit: 128 MB
Submit: 875   Solved: 499
[ Submit][ Status][ Discuss]

Description

ftiasch 有 N 个物品, 体积分别是 W₁, W₂, ..., W_N。 由于她的疏忽, 第 i 个物品丢失了。 “要使用剩下的 N - 1 物品装满容积为 x 的背包，有几种方法呢？” -- 这是经典的问题了。她把答案记为 Count(i, x) ，想要得到所有1 <= i <= N, 1 <= x <= M的 Count(i, x) 表格。

Input

 

第1行：两个整数 N (1 ≤ N ≤ 2 × 10³) 和 M (1 ≤ M ≤ 2 × 10³)，物品的数量和最大的容积。

第2行： N 个整数 W₁, W₂, ..., W_N, 物品的体积。

Output

 

一个 N × M 的矩阵， Count(i, x)的末位数字。

Sample Input

3 2
1 1 2

Sample Output

11
11
21

HINT

如果物品3丢失的话，只有一种方法装满容量是2的背包，即选择物品1和物品2。

Source

[ Submit][ Status][ Discuss] 
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;
}