Python编程之欧拉计划33~34

欧拉计划

欧拉计划（Project Euler）。

33、消去数字的分数

分数 49/98 是一个奇怪的分数：当一个菜鸟数学家试图对其进行简化时，他可能会错误地可以认为通过将分子和分母上的

9 同时去除得到 49/98 = 4/8。但他得到的结果却是正确的。

我们将 30/50 = 3/5 这样的分数作为普通个例。

一共有四个这样的非普通分数，其值小于 1，并且包括分子和分母都包括 2 位数。

如果将这四个分数的乘积约分到最简式，分母是多少？
第一个拥有超过500个约数的三角形数是多少？

分析：AB/BC = A/C为真分数说明B>A,且C>A。

ans=1
for i in range(1,10):
    for j in range(i+1,10):
        for k in range(i+1,10):
            a,b=10*i+j,j*10+k
            if a*k==b*i:
                print(f'{a}/{b}={i}/{k}')
                ans *=k/i
print(f'最简式分母为{int(ans)}')

输出：
16/64=1/4
19/95=1/5
26/65=2/5
49/98=4/8
最简式分母为100

34、数字阶乘

145 是一个奇怪的数字, 因为 1! + 4! + 5! = 1 + 24 + 120 = 145.

找出所有等于各位数字阶乘之和的数字之和。

注意: 因为 1! = 1 和 2! = 2 不是和的形式，所以它们不算在内。

分析：7个9组成的数的各位数阶乘和为7位，8个9也是七位，可知符合条件的数不会超过7位数。实际上 $9!=362880\ast 7=2540160$

# 个位数的阶乘表
f=(1, 1, 2, 6, 24, 120, 720, 5040, 40320, 362880)
ans = []

for i in range(10, 2540160):
    if i == sum(map(lambda x:f[int(x)],list(str(i)))):
        ans.append(i)

print(f'符合条件的数字有{ans},和为{sum(ans)}')

输出：
符合条件的数字有[145, 40585],和为40730