#-*- coding:utf-8 -*- import math from datetime import datetime def Euler63(): count = 0 dim = 1 while (math.log10(9 ** dim) < dim) & (math.log10(9 ** dim) >= (dim -1)): for i in xrange(1,10): if (math.log10(i ** dim) < dim) & (math.log10(…

素因数の和が互いに等しい、差が１の自然数の組をルース・アーロン・ペアと呼ぶ。 ２００００以下には２６組しかない。#素因数分解を原始的な方法で解く def prime_desolve(n): List = [] i = 2 while n > 1: if n % i == 0: while n % i == 0: List.append(…

ある4桁の自然数は、２つの数の立方体（三乗）の和として表す方法が２通りある最小の数である。 def Pra_13(): List = [] for i in xrange(1, 22): for j in xrange(i, 22): total = i **3 + j ** 3 if total in List: return total break else: List.append…