【HDU1085】母函数入门

1.题目链接。题目大意：给你一定数量的1元，2元，5元的钱，问你最下不能够凑出来的数目是哪一个？

2.额，显然是一个组合的问题了。这种选数据的题目，是母函数的最好的使用场所了。至于母函数是什么？我就不再一一推理说明了。其实很多人都只是会用，并不知道为什么。可能母函数的证明方法有很多种，但是我知道的一种证明使用了离散数学中的析取和合取这些基本的概念。所以要是证明又需要从盘古开天辟地说起，我就偷个懒得了，嘿嘿。

3.不多说了，直接给出代码，注释很写的很明白了。

#include<bits/stdc++.h>
#include<iostream>
#include<algorithm>
using namespace std;
const int MAX = 1e4 + 10;
int n1[MAX];//起点
int n2[MAX];//终点
int v[MAX];//权值
int P;//最大的组合
int a[MAX], b[MAX];
int n;//物品的种类
//母函数模板：init数组。作用：输入数量和权值，得到系数数组。
void init(int n,int n1[],int n2[],int v[], int a[])
{
	//最大值
	int P = 0;
	for (int i = 0; i < n; i++)
	{
		P += n2[i] * v[i];
	}
	/*如果n2是无穷，那么第二层循环条件j<=n2[i]可以去掉。*/
	int b[MAX];//中转数组
	for (int i = 0; i <n; i++)//循环每个因子
	{
		memset(b, 0, sizeof(b));
		for (int j = n1[i]; j <= n2[i] && j*v[i] <= P; j++)//循环每个因子的每一项
			for (int k = 0; k + j * v[i] <= P; k++)//循环a的每个项
				b[k + j * v[i]] += a[k];//把结果加到对应位
		memcpy(a, b, sizeof(b));//b赋值给a
	}
}
#pragma warning(disable:4996)
int main()
{
	v[0] = 1; v[1] = 2; v[2] = 5;//代表权值
	n = 3;//物品种类赋值
	while (scanf("%d%d%d", &n2[0], &n2[1], &n2[2]) && (n2[0] != 0 || n2[1] != 0 || n2[2] != 0))
	{
		memset(a, 0, sizeof(a));
		a[0] = 1;
		init(n, n1, n2, v,a);
		for (int i = 0; i <=MAX; i++)
		{
			if (!a[i])
			{
				cout << i << endl;
				break;
			}
		}
	}
	return 0;
}

后边的我的测试数据，我想说最后的判断条件，为什么是那个样子的。我们知道，a数组存的就是每一项的系数，如果这一项的系数是0，那么很不幸，说明凑不出来。而且，注意那里循环的上界，不是P，至少是P+1，因为从0到P可能都被赋值了，所以最小的不存在就是最大值加一了，开始没想好，在这里WA了一发。