【普通型母函数】HDU1085[Holding Bin-Laden Captive!]题解

题目概述

有 $num_1$ 个 $1$ 元硬币，$num_2$ 个 $2$ 元硬币，$num_3$ 个 $5$ 元硬币，求不能组成的最小价值。

解题报告

显然是二进制拆分做01背包，然而数据范围比较大，所以效率不是很好。

其实这道题可以用母函数轻松解决：（摘自维基和这里）

在数学中，某个序列的母函数是一种形式幂级数，其每一项的系数可以提供关于这个序列的信息。

也就是说母函数是针对某个序列的，它的外在表现形式是一种形式幂级数。比如说有这样一个序列 $a_0,a_1,\cdots,a_n$ ，构造一个函数 $f(x)=a_0+a_1x+a_2x^2+\cdots+a_nx^n$ 。则 $f(x)$ 是序列 $a_0,a_1,\cdots,a_n$ 的母函数。比如说最常见的 $(1+x)^n$ ，它是序列 ${n\choose0},{n\choose1},{n\choose2}\cdots{n\choose n}$ 的母函数。

由此可见母函数各项系数就是各项的出现次数。对于一种面值为 $c$ 的硬币，我们可以将 $kc$ 作为幂（ $k$ 表示选几次），则母函数为：

$$f(x)=(1+x+x^2+\cdots+x^{num_1})(1+x^2+x^4+\cdots+x^{2num_2})(1+x^5+x^{10}+\cdots+x^{5num_5})$$
只需要多项式乘法求出母函数各项系数，然后找到第一个系数为 0 <script type="math/tex" id="MathJax-Element-31">0</script> 的项就是答案。

示例程序

#include<cstdio>
using namespace std;
const int maxn=8000,val[3]={1,2,5};

int num[3],MAX,a[maxn+5],b[maxn+5];

int main(){
    freopen("program.in","r",stdin);
    freopen("program.out","w",stdout);
    while (~scanf("%d%d%d",&num[0],&num[1],&num[2])){
        if (!num[0]&&!num[1]&&!num[2]) break;
        num[1]*=val[1];num[2]*=val[2];MAX=num[0]+num[1]+num[2];
        for (int i=0;i<=MAX;i++) a[i]=i<=num[0],b[i]=0;
        for (int t=1;t<3;t++){
            for (int i=0;i<=MAX;i++)
                for (int j=0;i+j<=MAX&&j<=num[t];j+=val[t])
                    b[i+j]+=a[i];
            for (int i=0;i<=MAX;i++) a[i]=b[i],b[i]=0;
        }
        int ans=0;while (a[ans]&&ans<=MAX) ans++;printf("%d\n",ans);
    }
    return 0;
}