python动态规划经典题目_历届试题 包子凑数Python实现(完全背包，动规)

资源限制

时间限制：1.0s 内存限制：256.0MB

问题描述

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

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

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

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

输入格式

第一行包含一个整数N。(1 <= N <= 100)

以下N行每行包含一个整数Ai。(1 <= Ai <= 100)

输出格式

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

样例输入

2

4

5

样例输出

6

样例输入

2

4

6

样例输出

INF

样例说明

对于样例1，凑不出的数目包括：1, 2, 3, 6, 7, 11。

对于样例2，所有奇数都凑不出来，所以有无限多个。

数据规模和约定

峰值内存消耗(含虚拟机) < 256M

CPU消耗 < 1000ms

分析：

这个和之前做的“买不到的数目”很像，上一题是要求两个数字组合当中不能取到的最大的数字，而这一题，是要求，不能取到的数字有多少个。

首先分享一个不能AC的思路，一组数当中，不能取到的最大的数字最大是这组数字的最小公倍数，所以，我设一个数组长度是最小公倍数，然后再往数组里面排数字。本来是想固定数组长度，这样可以减少一定的复杂度，但是忽略了求一组数字的最小公倍数，这个循环起来是很可怕的，下面来看下代码：

不能AC的：

import math

while True:

try:

n = int(input())

s = []

for i in range(n):

s.append(int(input()))

#求最小公倍数

a,b = s[0],s[1]

a = a // math.gcd(a, b) * b // math.gcd(a, b)

if len(s)>2:

for i in range(2,len(s)):

b = s[i]

a = a//math.gcd(a,b) * b//math.gcd(a,b)

dp = [0 for i in range(a+1)]

dp[0] = 1

for x in range(n):

for y in range(min(s),len(dp)):

if dp[y-s[x]]==1:

dp[y] = 1

print(dp.count(0))

except:

break

但是完全背包的思想是对的，转换一下思路，省去求最小公倍数这个过程，我们看到，AiAiA_iAi​最大是100，那么最大的取不到的100∗100100∗100100*100100∗100肯定是够的，所以我们给动规的数组设置长度10000.

一组数字，如果他们的最大公约数不是1，那么肯定有无穷多个数字无法组合出来，所以我们要先求一下一组数字当中的最大公约数。

求完最大公约数之后我们来举个例子：以题目中给的两个数，4和5.

dp

0

1

2

3

4

5

6

7

8

9

10

11

12

…

1

0

0

0

0

0

0

0

0

0

0

0

0

…

动规数组建立如上。

然后看到4，将所有能用4组合起来的数组的位置置为1.

dp

0

1

2

3

4

5

6

7

8

9

10

11

12

…

1

0

0

0

1

0

0

0

1

0

0

0

1

…

再加入5，dp[j]=1dp[j]=1dp[j]=1dp[j]=1表示jjjj这个数字能够凑出来，反之dp[j]=0dp[j]=0dp[j]=0dp[j]=0表示凑不出来。

后面做一个判断：

if y >= s[i] and dp[y - s[i]] == 1:

dp[y] = 1

dp

0

1

2

3

4

5

6

7

8

9

10

11

12

…

1

0

0

0

1

1

0

0

1

1

1

0

1

…

后面的都能凑出来，所以不写了。

可以看出，数组中有6个零，所以凑不出来的个数是6个。

AC代码：

import math

while True:

try:

n = int(input())

s = []

g = 0

dp = [0 for i in range(10000)]

dp[0] = 1

for i in range(n):

s.append(int(input()))

if i == 0:

g = s[i]

else:

g = math.gcd(s[i], g) # 在输入的时候求出最大公约数

for y in range(1, len(dp)):

if y >= s[i] and dp[y - s[i]] == 1:

dp[y] = 1

if g == 1:

print(dp.count(0))

else:

print('INF')

except:

break

小白欢迎指点~

原文链接:https://blog.csdn.net/bianxia123456/article/details/106032268