beldmit: (Программизм)
Dmitry Belyavskiy ([personal profile] beldmit) wrote2009-12-18 10:40 pm

Занимательная математика - последнее про три девятки

У Перельмана приведено начало числа, вычисленное по какой-то сильно многозначной таблице логарифмов.

Так вот, ошибка - в последней приведенной у Перельмана цифре: 428124773175747048036987115 вместо 428124773175747048036987118 у Перельмана.
Общее число цифр, как я понимаю, 369693100 (у Перельмана указано 369693061).

Пытаясь загрузить этот файл, vim отожрал больше гига памяти и трудился почти час.
ext_613079: Default userpic (Default)

[identity profile] shaplov.livejournal.com 2009-12-18 08:31 pm (UTC)(link)
Вам, что, мегагерцы девать некуда?!
;-)

[identity profile] beldmit.livejournal.com 2009-12-18 08:39 pm (UTC)(link)
А куда?..

Коля, вычислить это число было одной из моих детских мечт. Я серьезно.
ext_613079: Default userpic (Default)

[identity profile] shaplov.livejournal.com 2009-12-18 08:42 pm (UTC)(link)
;-)
Понимаю...

Я лишь из любви к искусству ворчу... А на самом деле тайно завидую ;-)

[identity profile] beldmit.livejournal.com 2009-12-18 08:46 pm (UTC)(link)
Я уже года 3 подозревал, что современные компы это осилят за разумное время. Просто в этом году впервые получил инструмент, позволяющий не заморачиваться собственной реализацией длинной арифметики.

Теперь еще получил игрушку - llvm и pure, надо будет почитать, что они умеют еще, кроме как орехи колоть.
vitus_wagner: My photo 2005 (Default)

[personal profile] vitus_wagner 2009-12-18 09:40 pm (UTC)(link)
Ага. Несколько лет занимался разработкой такого инструмента, и не осознавал, что его получил.

[identity profile] beldmit.livejournal.com 2009-12-19 07:06 am (UTC)(link)
Не, на OpenSSL-ный BIGNUM я облизывался, но с ним же надо разбираться. А инструмент - это калькулятор.

[identity profile] slobin.livejournal.com 2009-12-19 01:15 pm (UTC)(link)
Ну да, считает-то там на самом деле GMP. Но читать её документацию и .h-файлы ломы, а тут простой и наглядный интерфейс. Я вот только не понял, автор pure в процессе перехода с GPL на LGPL интерфейс к GSL (которая там матрицы раньше считала) сделал необязательным или вообще невозможным? Если второе, то обидно: там, говорят, тоже хорошие алгоритмы.

... Estu tiel cxi ...

[identity profile] pashap.livejournal.com 2009-12-18 08:52 pm (UTC)(link)
>
Общее число цифр, как я понимаю, 369693100 (у Перельмана указано 369693061).

Забавно. Число цифр с точностью до 1 оценивается без проблем при наличии виндового калькулятора. Даже странно, что у Перельмана тут ошибка...
vitus_wagner: My photo 2005 (Default)

[personal profile] vitus_wagner 2009-12-18 09:41 pm (UTC)(link)
Ну не было у Перельмана виндового калькулятора.

[identity profile] pashap.livejournal.com 2009-12-18 09:55 pm (UTC)(link)
Девять в девятой можно и на бумажке посчитать. А таблицы логарифмов у него были. Судя по всему - недостаточно точные.

[identity profile] beldmit.livejournal.com 2009-12-19 07:08 am (UTC)(link)
200, что ли, знаков, самые длинные таблицы логарифмов... Я удивляюсь, что ошибка только в 27-й цифре.

Оценка с более высокой точностью (и более современная) мне тоже в какой-то из книг по занимательной математике попадалась.
vitus_wagner: My photo 2005 (Default)

[personal profile] vitus_wagner 2009-12-19 08:13 am (UTC)(link)
Попробовал я это посчитать через логарифмы в bc со scale=200. Получил гораздо больше точных цифр.

[identity profile] pashap.livejournal.com 2009-12-19 08:22 am (UTC)(link)
Собственно, это мне и странно - что в самом числе ошибка только в 27 знаке, а в числе цифр - всего лишь в 8-м.

[identity profile] kypexin.livejournal.com 2009-12-18 10:44 pm (UTC)(link)
"vim отожрал больше гига памяти и трудился почти час"

А ПЕРЕЛЬМАНУ КАКОВО БЫЛО?!

[identity profile] cmike.livejournal.com 2009-12-19 08:47 am (UTC)(link)
bc? А что считаем? 9**(9**9)?

[identity profile] beldmit.livejournal.com 2009-12-19 09:07 am (UTC)(link)
Ага.

[identity profile] slobin.livejournal.com 2009-12-19 01:19 pm (UTC)(link)
В первом посте Белявского на эту тему, в моих комментариях, идёт сравнение четырёх разных длинных арифметик (bc, J, clisp и pure). С bc мы оценивали время счёта порядка месяца. Поэтому Белявский так и обрадовался, когда я его на pure навёл. ;-) Ну то есть собственно считает, конечно, не pure, а библиотека GMP, но язык позволяет не изучать её документацию и .h-файлы, а просто набрать выражение и посчитать.

... В ответ они чем-то мигнули ...

[identity profile] maksa.livejournal.com 2009-12-19 11:51 am (UTC)(link)
Ничего не понимаю. Определить точно два десятка первых цифр и ошибиться в итоге на 39 порядков? Как это так?

[identity profile] beldmit.livejournal.com 2009-12-19 12:15 pm (UTC)(link)
Сам не очень понимаю. Разве что количество цифр он определял не по таблице логарифмов. Да и вполне может быть, что эта глава честно скомпилирована из нескольких источников.