Поиск:

Ответ в темуСоздание новой темы Создание опроса
> Пи, Алгоритмы вычисления 
:(
    Опции темы
__vi
  Дата 23.1.2004, 11:28 (ссылка) | (нет голосов) Загрузка ... Загрузка ... Быстрая цитата Цитата


Опытный
**


Профиль
Группа: Участник
Сообщений: 301
Регистрация: 21.1.2004

Репутация: нет
Всего: -1



Извините может быть за глупый вопрос, но мне интересны алгоритмы вычесления пи.
Вот есть такой который позволяет вычислять пи с любой точностью?
PM MAIL   Вверх
MuToGeN
Дата 23.1.2004, 11:38 (ссылка) | (нет голосов) Загрузка ... Загрузка ... Быстрая цитата Цитата


Лесник
****


Профиль
Группа: Модератор
Сообщений: 4379
Регистрация: 15.8.2002
Где: Москва

Репутация: нет
Всего: 32



Хм... Да, тоже интересно


--------------------
Three pings for the token rings,
Five pings for the UNIX machines,
Hundred pings for the broken links,
One special ping to check them all
Through Simple Network Management Protocol!
PM MAIL ICQ   Вверх
shedon
Дата 23.1.2004, 12:11 (ссылка) | (нет голосов) Загрузка ... Загрузка ... Быстрая цитата Цитата


Эксперт
***


Профиль
Группа: Экс. модератор
Сообщений: 1209
Регистрация: 17.1.2003
Где: Нижнiй Новгородъ

Репутация: нет
Всего: 11



http://mathc.chat.ru/hist/pihist.htm

Это сообщение отредактировал(а) shedon - 23.1.2004, 12:13


--------------------
Programming is like sex: One mistake and you have to support it your lifetime
PM MAIL WWW ICQ   Вверх
__vi
Дата 23.1.2004, 13:27 (ссылка) | (нет голосов) Загрузка ... Загрузка ... Быстрая цитата Цитата


Опытный
**


Профиль
Группа: Участник
Сообщений: 301
Регистрация: 21.1.2004

Репутация: нет
Всего: -1



Спасибо!
PM MAIL   Вверх
PROme
Дата 23.1.2004, 14:11 (ссылка) | (нет голосов) Загрузка ... Загрузка ... Быстрая цитата Цитата


Опытный
**


Профиль
Группа: Участник
Сообщений: 346
Регистрация: 23.7.2003
Где: Kiev

Репутация: нет
Всего: 1



Очень-очень-очень советую разыскать №3 за Март 2002 журнал Hard'n'Soft (UA)

Там на 88-й странице начинается отличная статья "Пиратская тропинка к Пи"... в ней описыватся не только алгоритм поиска, но и возможности (что с помощью данного эдиного числа при достаточных вычислительных можностях можно закодировать (заархивировать) любую информацию любого размера)... вобще советую всем кто интересуются всяким таким программно-математическим почитать выпуски тех времен - было очень много отличнейших статей...

А вобще, в двух словах, насколько помню, ищется вписыванием многоугольника в круг - чем больше сторон у многоугольника - тем точнее получим значение Пи


--------------------
SEO-мастер
PM MAIL WWW   Вверх
MuToGeN
Дата 23.1.2004, 14:14 (ссылка) | (нет голосов) Загрузка ... Загрузка ... Быстрая цитата Цитата


Лесник
****


Профиль
Группа: Модератор
Сообщений: 4379
Регистрация: 15.8.2002
Где: Москва

Репутация: нет
Всего: 32



PROme
Это дело через пределы прогнать случайно нельзя? Изи я чего-то непонимаю...


--------------------
Three pings for the token rings,
Five pings for the UNIX machines,
Hundred pings for the broken links,
One special ping to check them all
Through Simple Network Management Protocol!
PM MAIL ICQ   Вверх
shedon
Дата 23.1.2004, 14:16 (ссылка) | (нет голосов) Загрузка ... Загрузка ... Быстрая цитата Цитата


Эксперт
***


Профиль
Группа: Экс. модератор
Сообщений: 1209
Регистрация: 17.1.2003
Где: Нижнiй Новгородъ

Репутация: нет
Всего: 11



Цитата
Очень-очень-очень советую разыскать №3 за Март 2002 журнал Hard'n'Soft (UA)

http://arbuz.ferghana.ru/z_pi.html


--------------------
Programming is like sex: One mistake and you have to support it your lifetime
PM MAIL WWW ICQ   Вверх
Chingachguk
Дата 23.1.2004, 15:04 (ссылка) | (нет голосов) Загрузка ... Загрузка ... Быстрая цитата Цитата


Эксперт
***


Профиль
Группа: Участник Клуба
Сообщений: 1232
Регистрация: 25.3.2002
Где: Москва

Репутация: 1
Всего: 18



Цитата

Это дело через пределы прогнать случайно нельзя? Изи я чего-то непонимаю...


Ну вот число e можно записать: limit (N->беск) (1+(1/N))^N.

А теперь запиши через предел отношение длины окружности (сумму сторон вписанного или описанного N-угольника) к радиусу.



--------------------
I don't like the drugs (but the drugs like me). M.Manson.
PM MAIL ICQ   Вверх
asmworm
Дата 24.1.2004, 14:29 (ссылка) | (нет голосов) Загрузка ... Загрузка ... Быстрая цитата Цитата


Шустрый
*


Профиль
Группа: Участник
Сообщений: 55
Регистрация: 20.1.2004

Репутация: 1
Всего: 1



Вот что-то откопал:

=====================================
THE TEN BILLIONTH HEXADECIMAL DIGIT of Pi is 9

By: Simon Plouffe, Peter Borwein and David Bailey.

The following is part of a paper titled "On The Rapid Computation of
Various Polylogarithmic Constants". The full text, as well as
Fortran code, is available in

http://www.cecm.sfu.ca/~pborwein/

as P123 under the link "Computing Pi and Related Matters".

ABSTRACT:

We give algorithms for the computation of the d-th digit of certain
transcendental numbers in various bases. These algorithms can be
easily implemented (multiple precision arithmetic is not needed),
require virtually no memory, and feature run times that scale nearly
linearly with the order of the digit desired. They make it feasible to
compute, for example, the billionth binary digit of log(2) or pi on a
modest work station in a few days run time.

Indeed we computed the 10 billionth hexadeximal digit of pi as well as
the billionth hexadecimal digits of pi^2, log(2) and log^2(2), the
billionth decimal digit of log (9/10) and the five billionth decimal
digit of log(1 - 10^{-96}).


These calculations rest on three observations. First, the d-th digit
of 1/n is "easy" to compute. Secondly, this scheme extends to
certain polylogarithm and arctangent series. Thirdly, very special
types of identities exist for certain numbers like pi, pi^2, log(2) and
log^2(2). These are essentially polylogarithmic ladders in an integer
base. A number of these identities that we derive in this work appear
to be new, for example the critical identity for the binary digits of
pi is:

Код

       infinity
        -----
        \         -n    /     4          2          1          1     \
pi =    )      16      |  ------- - -------  - ------- - -------   |
        /               \  8 n + 1   8  n + 4   8 n  + 5   8 n +  6 /
        -----
        n = 0

===cut===

Теперь может кто-то кинет реализацию вышесказанного? Я что-то в матеиатике не очень.

Это сообщение отредактировал(а) asmworm - 24.1.2004, 14:33
PM MAIL   Вверх
acp
Дата 24.1.2004, 19:03 (ссылка) | (нет голосов) Загрузка ... Загрузка ... Быстрая цитата Цитата


Опытный
**


Профиль
Группа: Участник
Сообщений: 389
Регистрация: 4.2.2003
Где: Владимир

Репутация: нет
Всего: 2



А что тут сложного в реализации?
Обычный ряд. Вычисляется в цикле.
PM WWW ICQ   Вверх
asmworm
  Дата 24.1.2004, 21:33 (ссылка) | (нет голосов) Загрузка ... Загрузка ... Быстрая цитата Цитата


Шустрый
*


Профиль
Группа: Участник
Сообщений: 55
Регистрация: 20.1.2004

Репутация: 1
Всего: 1



Цитата(acp @ 24.1.2004, 19:03)
А что тут сложного в реализации?
Обычный ряд. Вычисляется в цикле.

Хотел бы я узнать как ты вычислишь 10000000000000000000000000000000 знак confused.gif . Думаю тебе нужно что-то помощнее чем long double или extended biggrin.gif
PM MAIL   Вверх
shedon
Дата 24.1.2004, 21:48 (ссылка) | (нет голосов) Загрузка ... Загрузка ... Быстрая цитата Цитата


Эксперт
***


Профиль
Группа: Экс. модератор
Сообщений: 1209
Регистрация: 17.1.2003
Где: Нижнiй Новгородъ

Репутация: нет
Всего: 11



Цитата
Хотел бы я узнать как ты вычислишь 10000000000000000000000000000000 знак  . Думаю тебе нужно что-то помощнее чем long double или extended 

Создай динамический массив и делай знак до каторого у тебя хватит памяти smile.gif


--------------------
Programming is like sex: One mistake and you have to support it your lifetime
PM MAIL WWW ICQ   Вверх
Chingachguk
Дата 24.1.2004, 22:35 (ссылка) | (нет голосов) Загрузка ... Загрузка ... Быстрая цитата Цитата


Эксперт
***


Профиль
Группа: Участник Клуба
Сообщений: 1232
Регистрация: 25.3.2002
Где: Москва

Репутация: 1
Всего: 18



По поводу приведеннной asmworm заметки. Я не очень хорошо знаю язык (позже со словарем посмотрю), но я не уверен, что эта формула описывает ПИ в смысле его "определения" или же результатом аналитических выкладок и скорее является чем-то вроде результата приближения ~"полиномами" вычисленных результатов или их части.


--------------------
I don't like the drugs (but the drugs like me). M.Manson.
PM MAIL ICQ   Вверх
__vi
Дата 25.1.2004, 12:30 (ссылка) | (нет голосов) Загрузка ... Загрузка ... Быстрая цитата Цитата


Опытный
**


Профиль
Группа: Участник
Сообщений: 301
Регистрация: 21.1.2004

Репутация: нет
Всего: -1



Я вот дома нашел сорс какой-то ужасно написанной проги (видно писал математик), но пи вычисляет, алгоритм понять не возможно, после некоторых модификаций смог заставить её вычеслять очень далеко (100 000 знаков после запятой за... 90 минут (у меня машина в 970bogomips)) очень медлено. Протестирую все вышепредложенные алгоритмы и скажу какой быстрее.

(Я так же хотел 10^6 после запятой но за 4 часа не было результата так я выключил)
Кстати если кто-нибудь может проверить достоверность вычисленний пишите мне личное сообщение, я пришлю.

Это сообщение отредактировал(а) __vi - 25.1.2004, 14:34
PM MAIL   Вверх
asmworm
  Дата 25.1.2004, 18:21 (ссылка) | (нет голосов) Загрузка ... Загрузка ... Быстрая цитата Цитата


Шустрый
*


Профиль
Группа: Участник
Сообщений: 55
Регистрация: 20.1.2004

Репутация: 1
Всего: 1



__vi
В принципе можно с помощью Mathematica 5 проверить, если у тебя её нет то пришли мне на мыло результат или лучше всего его CRC32(только скажи в какой системе счисления) или что-то в этом роде(чтоб не слать много кб)

Это сообщение отредактировал(а) asmworm - 25.1.2004, 18:23
PM MAIL   Вверх
Ответ в темуСоздание новой темы Создание опроса
Правила форума "Алгоритмы"

maxim1000

Форум "Алгоритмы" предназначен для обсуждения вопросов, связанных только с алгоритмами и структурами данных, без привязки к конкретному языку программирования и/или программному продукту.


Если Вам понравилась атмосфера форума, заходите к нам чаще! С уважением, maxim1000.

 
1 Пользователей читают эту тему (1 Гостей и 0 Скрытых Пользователей)
0 Пользователей:
« Предыдущая тема | Алгоритмы | Следующая тема »


 




[ Время генерации скрипта: 0.1284 ]   [ Использовано запросов: 20 ]   [ GZIP включён ]


Реклама на сайте     Информационное спонсорство

 
По вопросам размещения рекламы пишите на vladimir(sobaka)vingrad.ru
Отказ от ответственности     Powered by Invision Power Board(R) 1.3 © 2003  IPS, Inc.