用 python 解 NOIP 竞赛题

以下为2017年全国 NOIP 提高组复赛的第1题:

 

 

怎么样,读完题是不是感觉特别懵。—— 我是谁,我在哪里?

接下来我们来慢慢解析这道让人摸不着头脑的竞赛题。

 

一、首先翻译一下题的意思:

  • 假设,现在人民币只有两个币种,面值你定
  • 对面值数字的要求是这两个数字不能有公约数(除1这个公约数外)
       例如:2和5
  • 问用这两个面值的钱,不能凑齐的最大的价钱是多少
  • 当然,前提是你有任意多的钱

是不是稍微要清晰一些了呢。

 

二、那我们接下来用python写一个程序来完成这道题:

整个过程分两步:

第一步:借助 python 找规律划范围

 
 1 # 先找出能凑出来的金额
 2 def myFunc(a, b):  # 输入a,b 两个互素的面值
 3     c = 1  # 从1开始找出能凑出的金额
 4     while True:  # 不断循环,电脑配置低的,请远离,前方危险
 5         for i in range(c):
 6             an01 = a * i
 7             for j in range(c):
 8                 an02 = b * j
 9                 if an01 + an02 == c:  # 一旦找到能凑出当前金额c的i和j,打印出来
10                     # print(c, "=", a, '*', i, '+', b, '*', j)
11                     print(c,end='  ')
12         c += 1  # 金额不断上涨,上不封顶
13 
14 
15 if __name__ == "__main__":
16     myFunc(3, 5)

结果为:

 

 

如果将面值设置为 7,5 呢

结果:

 

 

再如果换成 7,9

结果:

 

综上:

  •  我们可以发现,不可组合的面值均集中在靠前的位置,但有多靠前,具体又在哪个位置呢?
  •  我们姑且假定这个数字就在 两数的乘积 之内,而且事实也是这样的。大家可以多试几对数字,检验一下。

 

 

第二步:范围找到后,我们再来考虑用 python 找出范围内的不可组合的金额值:

  备注:上面的程序是一个死循环,需要手动结束程序,建议不懂操作的小伙伴谨慎运行(嘿嘿,你是不是已经入坑啦!)。但下面这个程序就不一样了,小伙伴们尽管去运行吧。

 

 1 # 找出两数乘积范围内的可组合数据
 2 def myFunc(a, b):
 3     c = a * b
 4     my_list = []  # 创建存放所有组合出来的金额值
 5 
 6     # 找寻过程 -- 不断对比
 7     for i in range(0, c):
 8         an01 = a * i
 9         for j in range(c):
10             an02 = b * j
11             if an01 + an02 <= c:  # 只找在乘积范围内的组合,节省运算次数
12                 my_list.append(an01 + an02)  # 将符合的金额添加进目标列表
13     return list(set(sorted(my_list)))  # 返回经过去重和排序的目标列表
14 
15 
16 # 找到最大的那个不能组合的金额
17 def getMax(a, b):
18     my_list = myFunc(a, b)  # 调用找可拼凑数据函数得到目标列表
19     my_list.sort(reverse=True)  # 将目标列表反序排列
20 
21     # 判断目标列表是否连续,并输出断点数中的最大值
22     y = my_list[0] + 1  # 创建对比参数
23     for x in my_list:
24         if x + 1 != y:
25             print(x, y)
26             break
27         y = x
28     return y - 1  # 返回最大断点值
29 
30 
31 if __name__ == "__main__":
32     print(getMax(16, 27))

 


 结果为:

 

 

三、总结

不知道大家有没有发现一个问题,这个最大不可组合数据似乎有一定的规律,规律为:

 

  c = a * b - a - b

 

 

( 其中的a 和 b 为你输入的两个互为素数的币种面值,c为它们不能组合的金额 )
大家可以多试几组数据,验证一下。

 

而且我要悄悄告诉你的就是,这个公式可是一个牛哄哄的定理,名字叫:赛瓦维斯特定理

  赛瓦维斯特定理:

    已知a,b为大于1的正整数,(a,b)=1,则使不定方程 ax+by=c 无负整数解的最大整数c=abab   

    其中的 (a,b) 表示a和b的最大公约数

 

 

 

 

 

怎么样,通过两个程序,我们就很容易的解决了这个看起来不那么友好的竞赛题。
此时,是不是觉得 python 很酷呢!

笔者会不定时的更新一些跟python相关又和数学相关的一些有趣的程序,喜欢就关注我吧。

 

特别警示:本文为作者原创作品,禁止不经过本人同意就将其转载用于商业用途,否则将予以追究。
处于学习分享转载请附上出处链接,谢谢!
posted @ 2019-01-11 19:55  天使猿  阅读(3348)  评论(0编辑  收藏  举报