soj 3531 Number Pyramids(观察组合数 + 每种物品至少选一个的完全背包)

最新推荐文章于 2019-04-25 20:48:46 发布

羁绊残阳

最新推荐文章于 2019-04-25 20:48:46 发布

阅读量450

点赞数 3

分类专栏： ACM_动态规划

本文链接：https://blog.csdn.net/jibancanyang/article/details/49471949

版权

ACM_动态规划专栏收录该内容

121 篇文章 0 订阅

订阅专栏

题意
给你一个类似杨辉三角的堆积法,问在给定顶端的数和,底部长度的情况下,一共有多少种堆积方式?
分析
这个拿到找了半天分解子问题,记忆化搜索都没有办法…最后才发现了最底部的数,会被扩大,然后传递到顶部去,然后观察这个扩大倍数或发现恰恰是组合数!
假设最底部的数是: $a[1], a[2],....a[n]$ 那么最上面的数就是

C 0 n - 1 a [1] + C 1 n - 1 a [2] + . . . . . + C n n - 1 a [n]

$C_{n-1}^0a[1] + C_{n-1}^1a[2] +.....+C_{n-1}^na[n]$
然后就把这些组合数当做物品,他们的系数

a[i] $a[i]$ 当做选择的物品的个数来做完全背包就可以了.
然后这里比较特殊的是,每种物品都要至少选择一个,我们就先把背包总容量减去每个物品选一个所需要消耗的容量(

∑i=1...na[i] $\sum_{i =1 ...n}a[i]$ ),剩下的容量再来做完全背包.定义dp[i]为要达到这个容量有多少种选择方式,显然

dp[i]+=dp[i−Cin] $dp[i] += dp[i - C_n^i]$ ,

#include <cstdio>
#include <cstring>
#include <iostream>
#include <algorithm>
#include <set>
using namespace std;
const int mod = 1e9 + 9, maxn = 1e6 +9;
int n, m;
int dp[maxn];


int main(void) {
    while (~scanf("%d%d", &n, &m)) {
        if (n > 20) {
            puts("0");
            continue;
        }
        int x = 1 << (n - 1), y = n - 1;
        if (x > m) {
            puts("0");
            continue;
        }
        int c = 1;
        for (int i = 0; i <= y; i++) {
            m -= c;
            c = c * (y - i) / (i + 1);
        }
        fill(dp, dp + m + 9, 0);
        dp[0] = 1;
        for (int i = 0, c = 1; i <= y; i++) {
            for (int j = c; j <= m; j++) {
                dp[j] = (dp[j] + dp[j - c]) % mod;
            }
            c = c * (y - i) / (i + 1);
        }
        printf("%d\n", dp[m] % mod);
    }
    return 0;
}