Реализация формулы Чудновского для Pi

@STMx32 · Регистрация: 17.09.2023

Author24 — интернет-сервис помощи студентам

Собственно у знакомого дите горит с курсовиком.
Собственно первую часть он не без помощи сделал сделал. вторая часть - посчитать по формуле Чудновского первые 10М знаков и вывести в текстовый файл.

помогите разобрать исходник питона (какже без инета) для воспроизведения на delphi.

Python

def chudnovsky2(n):
    pi = Decimal(13591409)
    ak = Decimal(1)
    k = 1
    while k < n:
        ak *= -Decimal((6*k-5)*(2*k-1)*(6*k-1))/Decimal(k*k*k*26680*640320*640320)
 
        val = ak*(13591409 + 545140134*k)
        
        d = Decimal((6*k-5)*(2*k-1)*(6*k-1))/Decimal(k*k*k*26680*640320*640320)
        
        pi += val
        k += 1
    pi = pi * Decimal(10005).sqrt()/4270934400
    pi = pi**(-1)
    return pi

без либы для "Длинной математики" не обойтись.
насколько я понял Decimal это тип данных из "Длинной математики" - верно?
как разкладывается pi = pi**(-1)?

Чет совсем туплю на выходных.

@AzAtom · 21.01.2024, 18:14

** это видимо, возведение в степень.

@STMx32 · 21.01.2024, 19:13 **[ТС]**

Даааааа...... А препод не без юмора оказался...
итак свои 16 разрядов на double получил. теперь дело за библиотекой длинной математики.

какую посоветуете? да и вобще что есть?

PS алгоритм правда с изюминкой - правильно считает только если n кратно 2.

@STMx32 · 22.01.2024, 11:17 **[ТС]**

Мдя.... а выбор то lib для расчетов не большой или GMP или BigDecimal от Rudy.
пока попробывал BigDecimal - если без извращений типа 1000 знаков. прекрасно работает даже мудрить нечего не надо. указал новый тип и уперед. на большую разрядность пока не получается..... гдето косячу. правда демка про Pi жутко медленная. медленне примера из инета, где ошибка с integer....

попробую прикрутить GMP для делфей....

@snake32 · 22.01.2024, 18:55

Есть ещё такая формула, с поверхностным пояснением:
https://www.overclockers.at/ne... ort-thread
Используют int64 и DoubleDouble для деления
Умудряются миллиарды знаков после запятой считать в параллель

@Cyborg Drone · 23.01.2024, 06:11

Не настолько "в лоб", как в сообщении #1, но, опять же, на питоне:
Pi - Chudnovsky
Может, пригодится.

@Joeymax · 26.01.2024, 14:48

и еще один интересный проект (в плане рекорда) от fabrice bellard

@AzAtom Модератор 3756 / 2260 / 783 Регистрация: 15.11.2015 Сообщений: 8,994
	21.01.2024, 18:14	2
	** это видимо, возведение в степень. 1

@STMx32 0 / 0 / 0 Регистрация: 17.09.2023 Сообщений: 22
	21.01.2024, 19:13 [ТС]	3
	Даааааа...... А препод не без юмора оказался... итак свои 16 разрядов на double получил. теперь дело за библиотекой длинной математики. какую посоветуете? да и вобще что есть? PS алгоритм правда с изюминкой - правильно считает только если n кратно 2. 0

@STMx32 0 / 0 / 0 Регистрация: 17.09.2023 Сообщений: 22
	22.01.2024, 11:17 [ТС]	4
	Мдя.... а выбор то lib для расчетов не большой или GMP или BigDecimal от Rudy. пока попробывал BigDecimal - если без извращений типа 1000 знаков. прекрасно работает даже мудрить нечего не надо. указал новый тип и уперед. на большую разрядность пока не получается..... гдето косячу. правда демка про Pi жутко медленная. медленне примера из инета, где ошибка с integer.... попробую прикрутить GMP для делфей.... 0

@snake32 3460 / 1648 / 236 Регистрация: 26.02.2009 Сообщений: 8,051 Записей в блоге: 5
	22.01.2024, 18:55	5
	Есть ещё такая формула, с поверхностным пояснением: https://www.overclockers.at/ne... ort-thread Используют int64 и DoubleDouble для деления Умудряются миллиарды знаков после запятой считать в параллель 0

@Cyborg Drone Модератор 10060 / 5402 / 3353 Регистрация: 17.08.2012 Сообщений: 16,522
	23.01.2024, 06:11	6
	Не настолько "в лоб", как в сообщении #1, но, опять же, на питоне: Pi - Chudnovsky Может, пригодится. 4

@Joeymax 1180 / 421 / 107 Регистрация: 31.03.2012 Сообщений: 1,150
	26.01.2024, 14:48	7
	и еще один интересный проект (в плане рекорда) от fabrice bellard 0