鱼C论坛

 找回密码
 立即注册
查看: 5076|回复: 10

题目60:找出一个5个质数的集合,该集合中每两个质数相连接都能产生一个素数

[复制链接]
发表于 2022-10-29 15:45:16 | 显示全部楼层
  1. import time as t
  2. import math

  3. start = t.perf_counter()


  4. def check_prime(a_num):
  5.     is_prime = True
  6.     for prime_i in range(2, int(math.sqrt(a_num) + 1)):
  7.         if a_num % prime_i == 0:
  8.             is_prime = False
  9.     return is_prime


  10. def is_exchangeable(prime_list):
  11.     length_list = len(prime_list)
  12.     for each_prime1 in range(length_list - 1):
  13.         for each_prime2 in range(each_prime1 + 1, length_list):
  14.             new_num_1 = int(str(prime_list[each_prime1]) + str(prime_list[each_prime2]))
  15.             new_num_2 = int(str(prime_list[each_prime2]) + str(prime_list[each_prime1]))
  16.             if not check_prime(new_num_1) or not check_prime(new_num_2):
  17.                 return False
  18.     return True


  19. def p60():
  20.     primes = []
  21.     prime_range = 10000
  22.     for numbers in range(2, prime_range):
  23.         if check_prime(numbers):
  24.             primes.append(numbers)
  25.     primes.remove(2)
  26.     primes.remove(5)
  27.     test_range = len(primes)
  28.     for prime1 in range(test_range):
  29.         for prime2 in range(prime1 + 1, test_range):
  30.             cur_plist = [primes[prime1], primes[prime2]]
  31.             if is_exchangeable(cur_plist):
  32.                 for prime3 in range(prime2 + 1, test_range):
  33.                     cur_plist = [primes[prime1], primes[prime2], primes[prime3]]
  34.                     if is_exchangeable(cur_plist):
  35.                         for prime4 in range(prime3 + 1, test_range):
  36.                             cur_plist = [primes[prime1], primes[prime2], primes[prime3], primes[prime4]]
  37.                             if is_exchangeable(cur_plist):
  38.                                 for prime5 in range(prime4 + 1, test_range):
  39.                                     cur_plist = [primes[prime1], primes[prime2], primes[prime3],
  40.                                                  primes[prime4], primes[prime5]]
  41.                                     if is_exchangeable(cur_plist):
  42.                                         return cur_plist


  43. answer = p60()
  44. print(answer)
  45. print(sum(answer))
  46. print("It costs %f s" % (t.perf_counter() - start))
复制代码



[13, 5197, 5701, 6733, 8389]
26033
It costs 51.725813 s
小甲鱼最新课程 -> https://ilovefishc.com
回复 支持 反对

使用道具 举报

您需要登录后才可以回帖 登录 | 立即注册

本版积分规则

小黑屋|手机版|Archiver|鱼C工作室 ( 粤ICP备18085999号-1 | 粤公网安备 44051102000585号)

GMT+8, 2025-10-10 05:14

Powered by Discuz! X3.4

© 2001-2023 Discuz! Team.

快速回复 返回顶部 返回列表