feat(euler_37.py):修复combine_lists函数处理空列表的情况,并改进输出格式显示每个截断素数

📝 docs(0037.TruncatablePrimes):添加readme.md文档解释截断素数只有11个的原因
 feat(euler_38.py):新增全数字倍数问题的解决方案,包含优化算法和详细注释
This commit is contained in:
2026-01-08 14:38:53 +08:00
parent 5b1a9f0e4e
commit 745989cb32
3 changed files with 113 additions and 1 deletions

View File

@@ -26,6 +26,8 @@ def timer(func):
def combine_lists(a: list[int], b: list[int | None], c: list[int]) -> list[int]:
"""将三个列表的每个元素组合成数字"""
if b in [[],[None]] :
return [int(f"{x}{z}") for x, z in product(a, c)]
return [int(f"{x}{y}{z}") for x, y, z in product(a, b, c)]
@@ -80,7 +82,10 @@ def TruncatablePrime() -> list[int]:
@timer
def main():
print(sum(TruncatablePrime()))
nums = TruncatablePrime()
for num in nums:
print(num)
print(f"Answer: {sum(nums)}")
if __name__ == "__main__":

View File

@@ -0,0 +1,23 @@
# 截断素数
截断素数为什么只有11个
我看到很多人有这个疑问,我自己也是思考了很久,才恍然的。
首先向左向右分别截断一个素数得到的依然是素数这有多难呢最右的数位显然只能是3、7。
2、5会让这个向右截断的数编程可以合数并且截断到最后只有一位数时只能是2、3、5、7在排除前面说的2和5那么最右的数位只能是3、7。
相似的向左截断时最左的数位显然只能是2、3、5、7。
中间的数在向左截断时显然不能时双数那只能从1、3、5、7、9中选取而5结尾的数必然可以被5整除所以中间的数不能有5结尾。
这也就导致中间只有1、3、7、9这四个数可选。
到这我们了解了这个最基本的数字约束。
其次,由于截断的每一步都是一个素数,这就导致在高位数时可满足上述条件的素数会更为稀疏,
简单尝试验证就能发现从7位数开始满足上述约束的截断素数会更少且不会超过6位数的情况而6位数只有一个截断素数739397。
由于素数的分布是高度不均匀的,所以,而偏向连续位数的素数出现概率本身就会随素数本身的不均匀分布而更难出现。
这也能很好的解释,为什么没有五位数的截断素数。
总之,需要严格证明的话,我们可以根据 [【Angell & Godwin,1977】](https://www.ams.org/journals/mcom/1977-31-137/S0025-5718-1977-0427213-2/S0025-5718-1977-0427213-2.pdf) 简单扩展得到右截断素数在文章中以证明有83个且最大的为73939133。那么就可以简单枚举到这个最大的右截断素数中所有的左截断素数就是结果了。
-----
MiroMind的模型真的不错很好的总结了截断素数相关的论文研究 [【分享链接】](https://dr.miromind.ai/share/e9e297da-7d1b-49d3-b4a9-bdae9249c9aa) 。

View File

@@ -0,0 +1,84 @@
"""
Take the number 192 and multiply it by each of 1, 2, and 3:
192 * 1 = 192
192 * 2 = 384
192 * 3 = 576
By concatenating each product we get the 1 to 9 pandigital, 192384576.
We will call 192384576 the concatenated product of 192 and (1,2,3).
The same can be achieved by starting with 9 and multiplying by 1,2,3,4, and 5,
giving the pandigital, 918273645, which is the concatenated product of 9 and (1,2,3,4,5).
What is the largest 1 to 9 pandigital 9-digit number that can be formed as the concatenated product
of an integer with (1,2,...,n) where n > 1?
"""
import time
def timer(func):
def wrapper(*args, **kwargs):
start_time = time.time()
result = func(*args, **kwargs)
end_time = time.time()
print(f"Execution time: {end_time - start_time:.6f} seconds")
return result
return wrapper
# def max_pandiMulti():
# max_num = 0
# key = [1,2,3,4,5]
# res = (0, [])
# ok = False
# for i in range(9999, 9, -1):
# mul_keys = [i * j for j in key]
# for j in range(5) :
# numx = "".join(map(str, mul_keys[:j+1]))
# if len(numx) == 9 and set(numx) == set('123456789'):
# max_num = max(max_num, int(numx))
# if max_num == int(numx):
# res = (i, key[:j+1])
# ok = True
# break
# if ok:
# break
# return max_num, res
def max_pandiMulti_better():
max_num = 0
res = (0, [])
# 配置列表:(n, start, end),其中 start > end表示从大到小遍历
configs = [
(2, 9999, 5000), # n=2: i ∈ [5000, 9999]
(3, 333, 100), # n=3: i ∈ [100, 333]
(4, 33, 25), # n=4: i ∈ [25, 33]
(5, 9, 5) # n=5: i ∈ [5, 9]
]
for n, start, end in configs:
for i in range(start, end - 1, -1):
# 生成拼接字符串
s = ''.join(str(i * k) for k in range(1, n + 1))
# 检查是否为1-9 pandigital
if set(s) == set('123456789'):
num = int(s)
if num > max_num:
max_num = num
res = (i, list(range(1, n + 1)))
break # 当前n中i递减第一个满足条件的即为最大可终止内层循环
return max_num, res
@timer
def main():
maxnum, ress = max_pandiMulti_better()
print(f"{maxnum} : {ress}")
if __name__ == "__main__":
main()