包子凑数 (完全背包求数量问题)

  • Post author:
  • Post category:其他



1.题目引入:

小明几乎每天早晨都会在一家包子铺吃早餐。他发现这家包子铺有 NN 种蒸笼,其中第 ii 种蒸笼恰好能放 A_i 个包子。每种蒸笼都有非常多笼,可以认为是无限笼。

每当有顾客想买 X个包子,卖包子的大叔就会迅速选出若干笼包子来,使得这若干笼中恰好一共有 X 个包子。比如一共有 3 种蒸笼,分别能放 3、4 和 5 个包子。当顾客想买 11 个包子时,大叔就会选 2 笼 3 个的再加 1 笼 5 个的(也可能选出 1 笼 3 个的再加 2 笼 4 个的)。

当然有时包子大叔无论如何也凑不出顾客想买的数量。比如一共有 3 种蒸笼,分别能放 4、5 和 6 个包子。而顾客想买 7 个包子时,大叔就凑不出来了。

小明想知道一共有多少种数目是包子大叔凑不出来的。

输入描述

第一行包含一个整数 NN ( 100≤N≤100)。

以下 N 行每行包含一个整数 Ai​ (100≤Ai​≤100)。

输出描述

一个整数代表答案。如果凑不出的数目有无限多个,输出 INF。


2.样例输出:

示例 1

输入

2
4
5

输出

6

样例说明

凑不出的数目包括:1, 2, 3, 6, 7, 11。

示例 2

输入

2
4
6

输出

INF

思路:

这题是一个完全背包的问题:

数学结论:


两个数a,b(当gcd=1时)(gcd 应该知道吧,它是用来求最大公约数的即:欧几里德算法,相当于我们的辗转相除法)最大不能表示出来的数是:(a−1)(b−1)−1


,(大家可以验证一下,只要大于这个数的一定可以用a,b两个质数表示出来)根据数据范围,直接枚举到10000。

有些小伙伴可能就要问了,我的0< n< 200 那我们的范围岂不就变成了:200 ^200;有没有这样想的,如果你这样想了先点个赞我们继续向下看。只有当两个数时它的范围是最大的,显然两个数能表示的范围我多个数所表示的范围自然更大了。大脑又短路了不是。

这里的状态定义也挺特殊的,不过很容易想到:



f[i][j]表示从前i个数中选,目标总质量为j,属性是能否达到j。

状态转移:按照第i个物品选的数量来进行划分:



f[i][j] == f[i−1][j]|f[i−1][j−w[i]]|f[i−1][j−2∗w[i]]|…|f[i−1][j−k∗w[i]];



(注意:因为是满足一项就说明可以装满,这里用到的是异或运算符)



完全背包的优化:



f[i][j−w[i]] = f[i−1][j−w[i]]|f[i−1][j−2∗w[i]]|f[i−1][j−3∗w[i]]|…|f[i−1][j−k∗w[i]];



优化:



f[i][j] == f[i−1][j]|f[i][j−w[i]];



到一维:



f[j]|=f[j−w[i]];



到这里我们发现他其实就是一个完全背包的板子题,具体操作如下:


3.代码如下:

#include<iostream>
#include<cstring>
#include<algorithm>
using namespace std;

const int maxn=1e4+10,m=1e4;
int a[maxn],n;
bool f[maxn];
int gcd(int x,int y)
{
	return y?gcd(y,x%y):x;
 } 
int main()
{
	int q=0,ans=0;
	cin>>n;
	for(int i=1;i<=n;i++)
	{
		cin>>a[i];
		q=gcd(q,a[i]);
	}
	if(q!=1)
	{
		puts("INF");
		return 0;
	}
	f[0]=1;  // 确定边界,当背包没有空间时表示符合要求 
	for(int i=1;i<=n;i++)
	{
		for(int j=a[i];j<=m;j++)
		{
			f[j]|=f[j-a[i]];
		}
	}
	for(int i=0;i<=m;i++)
	{
		if(!f[i])
		{
			ans++;
		}
	}
	cout<<ans<<"\n";
	return 0;
}



版权声明:本文为weixin_52914088原创文章,遵循 CC 4.0 BY-SA 版权协议,转载请附上原文出处链接和本声明。