用筛法找出所有的素因子，再从中找出满足条件的n，结果是$n=8319823$（此时$\varphi(n)=8313928$）。
target=10**7
d=[[] if i%2 else [2] for i in range(target)]
p,q,n=target,1,3
while n < target:
if d[n]:
m=n
for k in d[n]:
m=m*(k-1)/k
if sorted(list(str(n)))==sorted(list(str(m))) and n*q < m*p:
p,q=n,m
else:
for k in range(n,target,n):
d[k].append(n)
n+=1
print p,q
|
