pythonで99の95乗が計算できない

Question

pythonで多倍長を用いて99の95乗を計算したいと思っております。
そこで繰り返し二乗で計算する次のコードを書きました。
結果は
2762516676992086031162215252099373937751823498820663686038954381955181420661434999029402997104757596041430934979321430477936990632890966956832566192554952764119325614305369446826420895442039482384803123768496541126881538692265003099627365446602458705267201
となりまたが、wolfram alphaで計算した値3848960788934848611927795802824596789608451156087366034658627953530148126008534258032267383768627487094610968554286692697374726725853195657679460590239636893953692985541958490801973870359499
と一致しません。問題点・改善方法の分かる方いらっしゃったらお願いいたします。

def pow(a,n):
    ret = 1
    a=int(a)
    n=int(n)
    while n > 0:
        if n % 2 == 1:
            ret = ret * a
        n = int(n / 2)
        a = a * a
    print(a)
    return ret

def run():
    print(pow(99,95))

run()

Comments

[退会済みユーザー]

問題点はintの精度で、from decimal import *とした後に、a=Decimal(a)、n=Decimal(n)とすると誤差は減りますが、答えが「3.848960788934848611927795803E+189」となって大事な部分が欠けてしまいます。もっと詳しい方の意見が出ますように。

参考まで。

Answer 1

本日似たような質問がありました

n = int(n / 2)

の n / 2 が浮動小数演算で誤差を生じさせます。

n = n // 2

に変更すると正しい結果になります。