ProjectEuler73.py

#!/usr/bin/env python
# coding:utf-8
"""
Project Euler Problem 73

nとdを正の整数として, 分数 n/d を考えよう. n<dかつHCF(n,d)=1のとき, 真既約分数と呼ぶ.
d ≤ 8について既約分数を大きさ順に並べると, 以下を得る:
1/8, 1/7, 1/6, 1/5, 1/4, 2/7, 1/3, 3/8, 2/5, 3/7, 1/2, 4/7, 3/5, 5/8, 2/3, 5/7, 3/4, 4/5, 5/6, 6/7, 7/8
1/3と1/2の間には3つの分数が存在することが分かる.
では, d ≤ 12,000 について真既約分数をソートした集合では, 1/3 と 1/2 の間に何個の分数があるか?
"""
import time
import Euler
t1 = time.time()

ans_list = []


def tagainiso(n, m):
    re = list(set(Euler.factoring(n)) & set(Euler.factoring(m)))
    if re == [1]:
        re = []
    return re

for d in range(1, 12001):
    if d % 50 == 0:
        print(d)
    for n in range(1, d):
        if 1.0 * n / d < 1.0 / 2 and 1.0 * n / d > 1.0 / 3:
            ans_list.append(1.0 * n / d)
    if d % 1000 == 0:
        ans_list = list(set(ans_list))

ans_list = list(set(ans_list))
answer = len(ans_list)
print(answer)
print(time.time() - t1, "seconds")