Suma schematów Bernoulliego oraz int too large to convert to float

Question

Muszę policzyć to:

import math

k = 1086
suma=0

for k in range(1086):
    a = math.comb(178322, k) * (1086/163861)**k * (1-1086/163861)**(178322-k)
    suma += a

print(suma)

Jest to suma schematów Bernoulliego i dostaję błąd OverflowError: int too large to convert to float. Próbowałem to przepisać używając "//", ale dostaję nieprawdopodobne wyniki - wszystkie składniki sumy poza ostatnim wynoszą zero. Jak to policzyć poprawnie?

Answer 1

Osiągasz granice w jakiej można zapisać wartość dla danego typu, dlatego masz taki błąd i to jest jasno napisane 

int too large to convert to float. 

musisz operować na innym rodzaju danych 

Float jest ograniczonej precyzji i JEST LIMITOWANY musisz zmienić rodzaj, najlepiej na decimal 

https://docs.python.org/3/library/decimal.html także pozostaje Ci zamienić dane na decimal 

tu masz nawet poczyniony artykuł o float i o problemie https://docs.python.org/3/tutorial/floatingpoint.html

 

Comments

Napisałem kod:

import math
import decimal

k = 1086
suma=0

for k in range(1086):
    a = math.comb(178322, k) * (decimal.Decimal(1106/163861))**k * (decimal.Decimal(1-1106/163861))**(178322-k)
    suma += a

print(suma)

Zwraca pewien wynik. Chyba działa poprawnie. Czy myślisz, że jest to napisane prawidłowo?