コード例 #1
0
ファイル: P80.py プロジェクト: gubenkoved/project-euler
def sqrt_fraction_part_with_custom_persicion(x, target_presicion):
    terms_amount = 1
    
    while True:
        terms_amount += 5
        expansion = sqrt_expansion.continued_fraction(x, terms_amount)
        evaluated = expansion.evaluate()
        denom_len = math.floor(math.log10(evaluated.denominator))

        f_part = int ( evaluated * Fraction( 10 ** denom_len ) )

        if f_part > 0 and math.log10(f_part) > target_presicion:
            return f_part
コード例 #2
0
ファイル: P66.py プロジェクト: gubenkoved/project-euler
from project_euler import *
import math
import sqrt_expansion

solution = { 'answer': -1, 'D': -1 }

for D in range(2, 1001):
    if is_square(D): continue

    n = 0
    while True:
        n += 1
        
        convergent = sqrt_expansion.continued_fraction(D, n).evaluate()
        y = convergent.denominator

        x_sq = 1 + D * y ** 2

        if is_square(x_sq):
            x = maximal_root(x_sq)
            
            if x > solution['answer']:
                solution['answer'] = x
                solution['D'] = D

            correct = (x ** 2 - D * y ** 2 == 1)

            print("{3:6} {0:40d}^2 - {1:4d}*{2:40d}^2 = 1".format(x, D, y, 'OK' if correct else 'FAIL'))
            break

print(solution)