Введение
В статье присутствуют математические формулы, поэтому для чтения перейдите на сайт для их корректного отображения. Число \(\pi \) имеет богатую историю. Данная константа обозначает отношение длины окружности к ее диаметру.
В науке число \(\pi \) используют в любых расчетах, где есть окружности. Начиная от объема банки газировки, до орбит спутников. И не только окружности. Ведь в изучении кривых линий число \(\pi \) помогает понять периодические и колебательные системы. Например, электромагнитные волны и даже музыку.
В 1706 году в книге «Новое введение в математику» британского ученого Уильяма Джонса (1675-1749 гг.) для обозначения числа 3,141592… впервые была использована буква греческого алфавита \(\pi \). Это обозначение происходит от начальной буквы греческих слов περιϕερεια — окружность, периферия и περιµετρoς — периметр. Общепринятым обозначение стало после работ Леонарда Эйлера в 1737 году.
Геометрический период
Постоянство отношения длины любой окружности к её диаметру было замечено уже давно. Жители Междуречья применяли довольно грубое приближение числа \(\pi \). Как следует из древних задач, в своих расчетах они используют значение \(\pi ≈ 3 \).
Более точное значение для \(\pi \) использовали древние египтяне. В Лондоне и Нью-Йорке хранятся две части древнеегипетского папируса, который называют «папирус Ринда». Папирус был составлен писцом Армесом примерно между 2000-1700 гг. до н.э.. Армес в своем папирусе написал, что площадь круга с радиусом \(r\) равна площади квадрата со стороной, равной \(\frac{8}{9} \) от диаметра окружности \(\frac{8}{9} \cdot 2r \), то есть \(\frac{256}{81} \cdot r^2 = \pi r^2 \). Отсюда \(\pi = 3,16\).
Древнегреческий математик Архимед (287-212 гг. до н.э.) впервые поставил задачу измерения круга на научную почву. Он получил оценку \(3\frac{10}{71} < \pi < 3\frac{1}{7}\), рассмотрев отношение периметров вписанного и описанного 96-угольника к диаметру окружности. Архимед выразил приближение числа \(\pi \) в виде дроби \(\frac{22}{7}\), которое до сих называется архимедовым числом.
Метод достаточно простой, но при отсутствии готовых таблиц тригонометрических функций потребуется извлечение корней. Кроме этого, приближение сходится к \(\pi \) очень медленно: с каждой итерацией погрешность уменьшается лишь вчетверо.
Аналитический период
Несмотря на это, до середины 17 века все попытки европейских учёных вычислить число \(\pi \) сводились к увеличению сторон многоугольника. Так например, голландский математик Лудольф ван Цейлен (1540-1610 гг.) вычислил приближенное значение числа \(\pi \) с точностью до 20-ти десятичных цифр.
На вычисление ему понадобилось 10 лет. Удваивая по методу Архимеда число сторон вписанных и описанных многоугольников, он дошел до \(60 \cdot 2^{29} \) — угольника с целью вычисления \(\pi \) с 20 десятичными знаками.
После смерти в его рукописях были обнаружены ещё 15 точных цифр числа \(\pi \). Лудольф завещал, чтобы найденные им знаки были высечены на его надгробном камне. В честь него число \(\pi \) иногда называли «лудольфовым числом» или «константой Лудольфа».
Одним из первых, кто представил метод, отличный от метода Архимеда, был Франсуа Виет (1540-1603 гг.). Он пришел к результату , что круг, диаметр которого равен единице, имеет площадь:
\[\frac{1}{2 \sqrt{\frac{1}{2}} \cdot \sqrt{\frac{1}{2} + \frac{1}{2} \sqrt{\frac{1}{2}} } \cdot \sqrt{\frac{1}{2} + \frac{1}{2} \sqrt{\frac{1}{2} + \frac{1}{2} \sqrt{\frac{1}{2} \cdots }}}} \]
С другой стороны, площадь равна \(\frac{\pi}{4} \). Подставив и упростив выражение, можно получить следующую формулу бесконечного произведения для вычисления приближенного значения \(\frac{\pi}{2} \):
\[\frac{\pi}{2} = \frac{2}{\sqrt{2}} \cdot \frac{2}{\sqrt{2 + \sqrt{2}}} \cdot \frac{2}{\sqrt{2+ \sqrt{2 + \sqrt{2}}}} \cdots \]
Полученная формула представляет собой первое точное аналитическое выражение для числа \(\pi \). Кроме этой формулы, Виет, используя метод Архимеда, дал с помощью вписанных и описанных многоугольников, начиная с 6-угольника и заканчивая многоугольником с \(2^{16} \cdot 6 \) сторонами приближение числа \(\pi \) с 9 правильными знаками.
Английский математик Уильям Броункер (1620-1684 гг.), используя цепную дробь , получил следующие результаты вычисления \(\frac{\pi}{4}\):
\[\frac{4}{\pi} = 1 + \frac{1^2}{2 + \frac{3^2}{2 + \frac{5^2}{2 + \frac{7^2}{2 + \frac{9^2}{2 + \frac{11^2}{2 + \cdots }}}}}} \]
Данный метод вычисления приближения числа \(\frac{4}{\pi} \) требует довольно больших вычислений, чтобы получить хотя бы небольшое приближение.
Получаемые в результате подстановки значения то больше, то меньше числа \(\pi \), и каждый раз все ближе к истинному значению, но для получения значения 3,141592 потребуется совершить довольно большие вычисления.
Другой английский математик Джон Мэчин (1686-1751 гг.) в 1706 году для вычисления числа \(\pi \) со 100 десятичными знаками воспользовался формулой, выведенной Лейбницем в 1673 году, и применил её следующим образом:
\[\frac{\pi}{4} = 4 arctg\frac{1}{5} — arctg\frac{1}{239} \]
Ряд быстро сходится и с его помощью можно вычислить число \(\pi \) с большой точностью. Формулы подобного типа использовались для установки нескольких рекордов в эпоху компьютеров.
В XVII в. с началом периода математики переменной величины наступил новый этап в вычислении \(\pi \). Немецкий математик Готфрид Вильгельм Лейбниц (1646-1716 гг.) в 1673 году нашел разложение числа \(\pi \), в общем виде его можно записать следующим бесконечным рядом:
\[ \pi = 1 — 4(\frac{1}{3} + \frac{1}{5} — \frac{1}{7} + \frac{1}{9} — \frac{1}{11} + \cdots) \]
Ряд получается при подстановке x = 1 в \(arctg x = x — \frac{x^3}{3} + \frac{x^5}{5} — \frac{x^7}{7} + \frac{x^9}{9} — \cdots\)
Леонард Эйлер развивает идею Лейбница в своих работах, посвященных использованию рядов для arctg x при вычислении числа \(\pi \). В трактате «De variis modis circuli quadraturam numeris proxime exprimendi» (О различных методах выражения квадратуры круга приближенными числами), написанном в 1738 году, рассматриваются методы усовершенствования вычислений по формуле Лейбница.
Эйлер пишет о том, что ряд для арктангенса будет сходиться быстрее, если аргумент будет стремиться к нулю. Для \(x = 1\) сходимость ряда очень медленная: для вычисления с точностью до 100 цифр необходимо сложить \(10^{50}\) членов ряда. Ускорить вычисления можно, уменьшив значение аргумента. Если принять \(x = \frac{\sqrt{3}}{3}\), то получается ряд
\[ \frac{\pi}{6} = artctg\frac{\sqrt{3}}{3} = \frac{\sqrt{3}}{3}(1 — \frac{1}{3 \cdot 3} + \frac{1}{5 \cdot 3^2} — \frac{1}{7 \cdot 3^3} + \cdots) \]
По утверждению Эйлера, если мы возьмем 210 членов этого ряда, то получим 100 верных знаков числа. Полученный ряд неудобен, потому что необходимо знать достаточно точное значение иррационального числа \(\sqrt{3} \). Также Эйлер в своих вычислениях использовал разложения арктангенсов на сумму арктангенсов меньших аргументов :
\[где x = n + \frac{n^2-1}{m-n}, y = m + p, z = m + \frac{m^2+1}{p} \]
Далеко не все формулы для вычисления \(\pi \), которые использовал Эйлер в своих записных книжках, были опубликованы. В опубликованных работах и записных книжках он рассмотрел 3 различных ряда для вычисления арктангенса, а также привел множество утверждений, касающихся количества суммируемых членов, необходимых для получения приближенного значения \(\pi \) c заданной точностью.
В последующие годы уточнения значения числа \(\pi \) происходили все быстрее и быстрее. Так, например, в 1794 году Георг Вега (1754-1802 гг.) определил уже 140 знаков , из который только 136 оказались верными.
Период компьютерных вычислений
XX век ознаменован совершенно новым этапом в вычислении числа \(\pi \). Индийский математик Сриниваса Рамануджан (1887-1920 гг.) обнаружил множество новых формул для \(\pi \). В 1910 году он получил формулу для вычисления \(\pi \) через разложение арктангенса в ряд Тейлора:
\[\pi = \frac{9801}{2\sqrt{2} \sum\limits_{k=1}^{\infty} \frac{(1103+26390k) \cdot (4k)!}{(4\cdot99)^{4k} (k!)^2}} .\]
При k=100 достигается точность в 600 верных цифр числа \(\pi \).
Появление ЭВМ позволило существенно увеличить точность получаемых значений за более короткие сроки. В 1949 году всего за 70 часов с помощью ENIAC группа ученых под руководством Джона фон Неймана (1903-1957 гг.) получила 2037 знаков после запятой числа \(\pi \) . Давид и Грегорий Чудновские в 1987 году получили формулу, с помощью которой смогли установить несколько рекордов в вычислении \(\pi \):
\[\frac{1}{\pi} = \frac{1}{426880\sqrt{10005}} \sum\limits_{k=1}^{\infty} \frac{(6k)!(13591409+545140134k)}{(3k)!(k!)^3(-640320)^{3k}}.\]
Каждый член ряда дает по 14 цифр. В 1989 году было получено 1 011 196 691 цифр после запятой. Данная формула хорошо подходит для вычисления \(\pi \) на персональных компьютерах. На данный момент братья являются профессорами в политехническом институте Нью-Йоркского университета.
Важным событием недавнего времени стало открытие формулы в 1997 году Саймоном Плаффом . Она позволяет извлечь любую шестнадцатеричную цифру числа \(\pi \) без вычисления предыдущих. Формула носит название «Формула Бэйли — Боруэйна — Плаффа» в честь авторов статьи, где формула была впервые опубликована. Она имеет следующий вид:
\[\pi = \sum\limits_{k=1}^{\infty} \frac{1}{16^k} (\frac{4}{8k+1} — \frac{2}{8k+4} — \frac{1}{8k+5} — \frac{1}{8k+6}) .\]
В 2006 году Саймон, используя PSLQ, получил несколько красивых формул для вычисления \(\pi \). Например,
\[ \frac{\pi}{24} = \sum\limits_{n=1}^{\infty} \frac{1}{n} (\frac{3}{q^n — 1} — \frac{4}{q^{2n} -1} + \frac{1}{q^{4n} -1}), \]
\[ \frac{\pi^3}{180} = \sum\limits_{n=1}^{\infty} \frac{1}{n^3} (\frac{4}{q^{2n} — 1} — \frac{5}{q^{2n} -1} + \frac{1}{q^{4n} -1}), \]
где \(q = e^{\pi}\). В 2009 году японские ученые, используя суперкомпьютер T2K Tsukuba System, получили число \(\pi \) c 2 576 980 377 524 десятичными знаками после запятой. Вычисления заняли 73 часа 36 минут. Компьютер был оснащен 640-ка четырех ядерными процессорами AMD Opteron, что обеспечило производительность в 95 триллионов операций в секунду.
Следующее достижение в вычислении \(\pi \) принадлежит французскому программисту Фабрису Беллару , который в конце 2009 года на своем персональном компьютере под управлением Fedora 10 установил рекорд, вычислив 2 699 999 990 000 знаков после запятой числа \(\pi \). За последние 14 лет это первый мировой рекорд, который поставлен без использования суперкомпьютера. Для высокой производительности Фабрис использовал формулу братьев Чудновских. В общей сложности вычисление заняло 131 день (103 дня расчеты и 13 дней проверка результата). Достижение Беллара показало, что для таких вычислений не обязательно иметь суперкомпьютер.
Всего через полгода рекорд Франсуа был побит инженерами Александром Йи и Сингеру Кондо. Для установления рекорда в 5 триллионов знаков после запятой числа \(\pi \) был также использован персональный компьютер, но уже с более внушительными характеристиками: два процессора Intel Xeon X5680 по 3,33 ГГц, 96 ГБ оперативной памяти, 38 ТБ дисковой памяти и операционная система Windows Server 2008 R2 Enterprise x64. Для вычислений Александр и Сингеру использовали формулу братьев Чудновских. Процесс вычисления занял 90 дней и 22 ТБ дискового пространства. В 2011 году они установили еще один рекорд , вычислив 10 триллионов десятичных знаков числа \(\pi \). Вычисления происходили на том же компьютере, на котором был поставлен их предыдущий рекорд и занял в общей сложности 371 день. В конце 2013 года Александр и Сингеру улучшили рекорд до 12,1 триллиона цифр числа \(\pi \), вычисление которых заняло у них всего 94 дня. Такое улучшение в производительности достигнуто благодаря оптимизации производительности программного обеспечения, увеличения количества ядер процессора и значительного улучшения отказоустойчивости ПО.
Текущим рекордом является рекорд Александра Йи и Сингеру Кондо, который составляет 12,1 триллиона цифр после запятой числа \(\pi \).
Таким образом, мы рассмотрели методы вычисления числа \(\pi \), используемые в древние времена, аналитические методы, а также рассмотрели современные методы и рекорды по вычислению числа \(\pi \) на компьютерах.
Список источников
- Жуков А.В. Вездесущее число Пи – М.:Изд-во ЛКИ, 2007 – 216 с.
- Ф.Рудио. О квадратуре круга, с приложением истории вопроса, составленной Ф.Рудио. / Рудио Ф. – М.: ОНТИ НКТП СССР, 1936. – 235c.
- Arndt, J. Pi Unleashed / J. Arndt, C. Haenel. – Springer, 2001. – 270p.
- Шухман, Е.В. Приближенное вычисление числа Пи с помощью ряда для arctg x в опубликованных и неопубликованных работах Леонарда Эйлера / Е.В. Шухман. — История науки и техники, 2008 – №4. – С. 2-17.
- Euler, L. De variis modis circuli quadraturam numeris proxime exprimendi/ Commentarii academiae scientiarum Petropolitanae. 1744 – Vol.9 – 222-236p.
- Шумихин, С. Число Пи. История длиною в 4000 лет / С. Шумихин, А. Шумихина. — М.: Эксмо, 2011. — 192с.
- Борвейн, Дж.М. Рамануджан и число Пи. / Борвейн, Дж.М., Борвейн П.Б. В мире науки. 1988 – №4. – С. 58-66.
- Alex Yee. Number world. Access mode: numberworld.org
Понравилось?
Расскажи
Если сравнить окружности отличных друг от друга размеров, то можно заметить следующее: размеры разных окружностей пропорциональны. А это значит, что при увеличении диаметра окружности в некоторое количество раз, увеличивается и длина этой окружности в такое же количество раз. Математически это записать можно так:
| C 1 | C 2 | ||
| = | |||
| d 1 | d 2 | (1) |
где C1 и С2 – длины двух разных окружностей, а d1 и d2 – их диаметры.
Это соотношение работает при наличии коэффициента пропорциональности – уже знакомой нам константы π
. Из отношения (1) можно сделать вывод: длина окружности C равна произведению диаметра этой окружности на независящий от окружности коэффициент пропорциональности π
:
C = π d.
Также эту формулу можно записать в ином виде, выразив диаметр d через радиус R данной окружности:
С = 2π R.
Как раз эта формула и является проводником в мир окружностей для семиклассников.
Еще с древности люди пытались установить значение этой константы. Так, например, жители Месопотамии вычисляли площадь круга по формуле:
Откуда π = 3.
В древнем Египте значение для π было точнее. В 2000-1700 годах до нашей эры писец, именуемый Ахмесом, составил папирус, в котором мы находим рецепты разрешения различных практических задач. Так, например, для нахождения площади круга он использует формулу:
| 8 | 2 | |||||
| S | = | ( | d | ) | ||
| 9 |
Из каких соображений он получил эту формулу? – Неизвестно. Вероятно, на основе своих наблюдений, впрочем, как это делали и другие древние философы.
По стопам Архимеда
Какое из двух числе больше 22/7 или 3.14 ?
- Они равны.
- Почему?
- Каждое из них равно π
.
А. А. Власов. Из Экзаменационного билета.
Некоторы полагают, что дробь 22/7 и чисо π тождественно равны. Но это является заблуждением. Помимо вышеприведенного неверного ответа на экзамене (см. эпиграф) к этой группе можно также добавить одну весьма занимательную головоломку. Задание гласит: "переложите одну спичку так, чтобы равенство стало верным".
Решение будет таковым: нужно образовать "крышу" для двух вертикальных спичек слева, используя одну из вертикальных спичек в знаменателе справа. Получится визуальное изображение буквы π .
Многие знают, что приближение π = 22/7 определил древнегреческий математик Архимед. В честь этого часто такое приближение называют "Архимедовым" числом. Архимеду удалось не только установить приближенное значение для π, но также найти точность этого приближения, а именно – найти узкий числовой промежуток, которому принадлежит значение π . В одной из своих работ Архимед доказывает цепь неравенств, которая на современный лад выглядела бы так:
| 10 | 6336 | 14688 | 1 | |||||||||
| 3 | < | < | π | < | < | 3 | ||||||
| 71 | 1 | 1 | 7 | |||||||||
| 2017 | 4673 | |||||||||||
| 4 | 2 | |||||||||||
можно записать проще: 3,140 909 < π < 3,1 428 265...
Как видим из неравенств, Архимед нашел довольно-таки точное значение с точностью до 0,002. Самое удивительно то, что он нашел два первых знака после запятой: 3,14... Именно такое значение чаще всего мы используем в несложных расчетах.
Практическое применение
Едут двое в поезде:
− Вот смотри, рельсы прямые, колеса круглые.
Откуда же стук?
− Как откуда? Колеса-то круглые, а площадь
круга пи эр квадрат, вот квадрат-то и стучит!
Как правило, знакомятся с этим удивительным числом в 6-7 классе, но более основательно им занимаются к концу 8-го класса. В этой части статьи мы приведем основные и самые важные формулы, которые пригодятся вам в решении геометрических задач, только для начала условимся принимать π за 3,14 для удобства подсчета.
Пожалуй, самая известная формула среди школьников, в которой используется π , это – формула длины и площади окружности. Первая – формула площади круга – записывается так:
| π D 2 | |
| S=π R 2 = | |
| 4 |
где S – площадь окружности, R – ее радиус, D – диаметр окружности.
Длина окружности, или, как ее иногда называют, периметр окружности, вычисляют по формуле:
С = 2 π R = π d,
где C – длина окружности, R – радиус, d – диаметр окружности.
Понятно, что диаметр d равен двум радиусам R.
Из формулы длины окружности можно легко найти радиус окружности:
где D – диаметр, С – длина окружности, R – радиус окружности.
Это базовые формулы, знать которые должен каждый ученик. Также иногда приходится вычислять площадь не всей окружности, а только ее части – сектора. Поэтому представляем вам её – формулу для вычисления площади сектора окружности. Выглядит она так:
| α | |||
| S | = | π R 2 | |
| 360 ˚ |
где S – площадь сектора, R – радиус окружности, α – центральный угол в градусах.
Такое загадочное 3,14
И правда, оно загадочно. Потому что в честь этих магических цифр устраивают праздники, снимают фильмы, проводят общественные акции, пишут стихи и многое другое.
Например, в 1998 году вышел фильм американского режиссера Даррена Аронофски под названием "Пи". Фильм получил множество наград.
Каждый год 14 марта в 1:59:26 люди, интересующиеся математикой, празднуют "День числа Пи". К празднику люди подготавливают круглый торт, усаживаются за круглый стол и обсуждают число Пи, решают задачи и головоломки, связанные с Пи.
Вниманием это удивительное число не обошли и поэты, неизвестный написал:
Надо только постараться и запомнить всё как есть – три, четырнадцать, пятнадцать, девяносто два и шесть.
Давайте развлечемся!
Вашему вниманию предлагаются интересные ребусы с числом Пи. Разгадайте слова, какие зашифрованы ниже.
1. π р
2. π L
3. π k
Ответы: 1. Пир; 2. Надпил; 3. Писк.
Недавно на Хабре в одной статье упомянули про вопрос «Что было бы с миром, если бы число Пи равнялось 4?» Я решил слегка поразмышлять на эту тему, используя некоторые (пусть и не самые обширные) знания в соответствующих областях математики. Кому интересно – прошу под кат.
Чтобы представить такой мир, нужно математически реализовать пространство с иным соотношением длины окружности к ее диаметру. Это я и попытался сделать.
Попытка №1.
Оговорим сразу, что рассматривать я буду только двумерные пространства. Почему? Потому что окружность, собственно, определена в двумерном пространстве (если рассмотреть размерность n>2, то отношение меры (n-1)-мерной окружности к ее радиусу даже не будет константой).Так что для начала я попытался придумать хоть какое-то пространство, где Пи не равно 3.1415… Для этого я взял метрическое пространство с метрикой, в которой расстояние между двумя точками равно максимуму среди модулей разности координат (т.е. расстояние Чебышева).
Какой же вид будет иметь единичная окружность в этом пространстве? Возьмем точку с координатами (0,0) за центр этой окружности. Тогда множество точек, расстояние (в смысле заданной метрики) от которых до центра равно 1, есть 4 отрезка, параллельных осям координат, образующих квадрат со стороной 2 и с центром в нуле.
Да, в некоторой метрике это - окружность!
Посчитаем здесь Пи. Радиус равен 1, тогда диаметр, соответственно, равен 2. Можно также рассмотреть определение диаметра как наибольшего расстояния между двумя точками, но даже так оно равно 2. Осталось найти длину нашей «окружности» в данной метрике. Это сумма длин всех четырех отрезков, которые в данной метрике имеют длину max(0,2)=2. Значит, длина окружности равна 4*2=8. Ну а тогда Пи здесь равно 8/2=4. Получилось! Но нужно ли сильно радоваться? Результат этот практически бесполезен, ведь рассматриваемое пространство абсолютно абстрактно, в нем даже не определены углы и повороты. Вы можете представить себе мир, где по факту не определен поворот, и где окружностью является квадрат? Я пытался, честно, но у меня не хватило воображения.
Радиус равен 1, а вот с нахождением длины этой «окружности» есть некоторые сложности. После некоторых поисков информации в интернете, я пришел к выводу, что в псевдоевклидовом пространстве такое понятие как «число Пи» вообще не может быть определено, что, безусловно, плохо.
Если кто-нибудь в комментариях расскажет мне, как формально считать длину кривой в псевдоевклидовом пространстве, я буду очень рад, ибо моих познаний в дифференциальной геометрии, топологии (а также усердного гугления) для этого не хватило.
Выводы:
Не знаю, можно ли писать о выводах после таких не сильно продолжительных исследований, но кое-что сказать можно. Во-первых, попытавшись представить пространство с иным числом Пи, я понял, что оно будет слишком абстрактно, чтобы быть моделью реального мира. Во-вторых, когда если попытаться придумать более удачную модель (похожую на наш, реальный мир), выходит, что число Пи останется неизменным. Если принять за данность возможность отрицательного квадрата расстояния (что для обычного человека - просто абсурд), то Пи не будет определено вовсе! Все это и наводит на мысль, что, возможно, мира с другим числом Пи и вовсе быть не могло? Ведь не зря же Вселенная именно такая, какая она есть. А может быть, это и реально, только обычной математики, физики и человеческого воображения для этого недостаточно. А вы как считаете?Upd. Узнал точно. Длина кривой в псевдоевклидовом пространстве может быть определена только на каком-либо его евклидовом подпространстве. То есть, в частности, для получившейся в попытке N3 «окружности» вовсе не определено такое понятие как «длина». Соответственно, Пи там тоже посчитать нельзя.
С недавних пор существует элегантная формула для вычисления числа Пи, которую в 1995 году впервые опубликовали Дэвид Бэйли, Питер Борвайн и Саймон Плафф:
Казалось бы: что в ней особенного — формул для вычисления Пи великое множество: от школьного метода Монте-Карло до труднопостижимого интеграла Пуассона и формулы Франсуа Виета из позднего Средневековья. Но именно на эту формулу стоит обратить особое внимание — она позволяет вычислить n-й знак числа пи без нахождения предыдущих. За информацией о том, как это работает, а также за готовым кодом на языке C, вычисляющим 1 000 000-й знак, прошу под хабракат.
Как же работает алгоритм вычисления N-го знака Пи?
К примеру, если нам нужен 1000-й шестнадцатеричный знак числа Пи, мы домножаем всю формулу на 16^1000, тем самым обращая множитель, стоящий перед скобками, в 16^(1000-k). При возведении в степень мы используем двоичный алгоритм возведения в степень или, как будет показано в примере ниже, возведение в степень по модулю . После этого вычисляем сумму нескольких членов ряда. Причём необязательно вычислять много: по мере возрастания k 16^(N-k) быстро убывает, так что, последующие члены не будут оказывать влияния на значение искомых цифр). Вот и вся магия — гениальная и простая.
Формула Бэйли-Борвайна-Плаффа была найдена Саймоном Плаффом при помощи алгоритма PSLQ , который был в 2000 году включён в список Top 10 Algorithms of the Century . Сам же алгоритм PSLQ был в свою очередь разработан Бэйли. Вот такой мексиканский сериал про математиков.
Кстати, время работы алгоритма — O(N), использование памяти — O(log N), где N — порядковый номер искомого знака.
Думаю, уместно будет привести код на языке Си, написанный непосредственно автором алгоритма, Дэвидом Бэйли:
/*
This program implements the BBP algorithm to generate a few hexadecimal
digits beginning immediately after a given position id, or in other words
beginning at position id + 1. On most systems using IEEE 64-bit floating-
point arithmetic, this code works correctly so long as d is less than
approximately 1.18 x 10^7. If 80-bit arithmetic can be employed, this limit
is significantly higher. Whatever arithmetic is used, results for a given
position id can be checked by repeating with id-1 or id+1, and verifying
that the hex digits perfectly overlap with an offset of one, except possibly
for a few trailing digits. The resulting fractions are typically accurate
to at least 11 decimal digits, and to at least 9 hex digits.
*/
/* David H. Bailey 2006-09-08 */
#include
Какие возможности это даёт? Например: мы можем создать систему распределённых вычислений, рассчитывающую число Пи и поставить всем Хабром новый рекорд по точности вычисления (который сейчас, к слову, составляет 10 триллионов знаков после запятой). Согласно эмпирическим данным, дробная часть числа Пи представляет собой нормальную числовую последовательность (хотя доказать это достоверно ещё не удалось), а значит, последовательности цифр из него можно использовать в генерации паролей и просто случайных чисел, или в криптографических алгоритмах (например, в хэшировании). Способов применения можно найти великое множество - надо только включить фантазию.
Больше информации по теме вы можете найти в статье самого Дэвида Бэйли, где он подробно рассказывает про алгоритм и его имплементацию (pdf);
И, похоже, вы только что прочитали первую русскоязычную статью об этом алгоритме в рунете - других я найти не смог.
ЧИСЛО ПИ
Символ ПИ означает отношение длины окружности к ее диаметру. Впервые в этом смысле символ p был использован У. Джонсом в 1707, а Л. Эйлер, приняв это обозначение, ввел его в научный обиход. Еще в древности математикам было известно, что вычисление значения p и площади круга - задачи, тесно связанные между собой. Древние китайцы и древние евреи считали число p равным 3. Значение числа p, равное 3,1605, содержится в древнеегипетском папирусе писца Ахмеса (ок. 1650 до н. э.). Около 225 до н. э. Архимед, используя вписанный и описанный правильные 96-угольники, приближенно вычислил площадь круга с помощью метода, который привел к значению ПИ, заключенному между 31/7 и 310/71. Другое приближенное значение p, эквивалентное обычному десятичному представлению этого числа 3,1416, известно еще со 2 в.
Л. ван Цейлен (1540-1610) вычислил значение ПИ с 32 десятичными знаками. К концу 17 в. новые методы математического анализа позволили вычислять значение p множеством различных способов. В 1593 Ф. Виет (1540-1603) вывел формулу
В 1665 Дж. Валлис (1616-1703) доказал, что
В 1658 У. Броункер нашел представление числа p в виде непрерывной дроби
Г.Лейбниц в 1673 опубликовал ряд
Ряды позволяют вычислять значение p с любым числом десятичных знаков. В последние годы с появлением электронных вычислительных машин значение p было найдено более чем с 10 000 знаков. С десятью знаками значение ПИ равно 3,1415926536. Как число, ПИ обладает некоторыми интересными свойствами. Например, его нельзя представить в виде отношения двух целых чисел или периодической десятичной дроби; число ПИ трансцендентно, т.е. непредставимо в виде корня алгебраического уравнения с рациональными коэффициентами. Число ПИ входит во многие математические, физические и технические формулы, в том числе и не имеющие непосредственного отношения к площади круга или длине дуги окружности. Например, площадь эллипса A определяется формулой A = pab, где a и b - длины большой и малой полуосей.
Энциклопедия Кольера. - Открытое общество . 2000 .
Смотреть что такое "ЧИСЛО ПИ" в других словарях:
число - Прие моч ное Источник: ГОСТ 111 90: Стекло листовое. Технические условия оригинал документа Смотри также родственные термины: 109. Число бетатронных колебаний … Словарь-справочник терминов нормативно-технической документации
Сущ., с., употр. очень часто Морфология: (нет) чего? числа, чему? числу, (вижу) что? число, чем? числом, о чём? о числе; мн. что? числа, (нет) чего? чисел, чему? числам, (вижу) что? числа, чем? числами, о чём? о числах математика 1. Числом… … Толковый словарь Дмитриева
ЧИСЛО, числа, мн. числа, чисел, числам, ср. 1. Понятие, служащее выражением количества, то, при помощи чего производится счет предметов и явлений (мат.). Целое число. Дробное число. Именованное число. Простое число. (см. простой1 в 1 знач.).… … Толковый словарь Ушакова
Абстрактное, лишенное особенного содержания обозначение какоголибо члена некоторого ряда, в котором этому члену предшествует или следует за ним какой нибудь др. определенный член; абстрактный индивидуальный признак, отличающий одно множество от… … Философская энциклопедия
Число - Число грамматическая категория, выражающая количественные характеристики предметов мысли. Грамматическое число одно из проявлений более обшей языковой категории количества (см. Категория языковая) наряду с лексическим проявлением («лексическое… … Лингвистический энциклопедический словарь
Число, приближенно равное 2,718, которое часто встречается в математике и естественных науках. Например, при распаде радиоактивного вещества по истечении времени t от исходного количества вещества остается доля, равная e kt, где k число,… … Энциклопедия Кольера
А; мн. числа, сел, слам; ср. 1. Единица счёта, выражающая то или иное количество. Дробное, целое, простое ч. Чётное, нечётное ч. Считать круглыми числами (приблизительно, считая целыми единицами или десятками). Натуральное ч. (целое положительное … Энциклопедический словарь
Ср. количество, счетом, на вопрос: сколько? и самый знак, выражающий количество, цифра. Без числа; нет числа, без счету, многое множество. Поставь приборы, по числу гостей. Числа римские, арабские или церковные. Целое число, ·противоп. дробь.… … Толковый словарь Даля
ЧИСЛО, а, мн. числа, сел, слам, ср. 1. Основное понятие математики величина, при помощи к рой производится счёт. Целое ч. Дробное ч. Действительное ч. Комплексное ч. Натуральное ч. (целое положительное число). Простое ч. (натуральное число, не… … Толковый словарь Ожегова
ЧИСЛО «Е» (ЕХР), иррациональное число, служащее основанием натуральных ЛОГАРИФМОВ. Это действительное десятичное число, бесконечная дробь, равная 2,7182818284590...., является пределом выражения (1/) при п, стремящемся к бесконечности. По сути,… … Научно-технический энциклопедический словарь
Количество, наличность, состав, численность, контингент, сумма, цифра; день.. Ср. . См. день, количество. небольшое число, несть числа, расти числом... Словарь русских синонимов и сходных по смыслу выражений. под. ред. Н. Абрамова, М.: Русские… … Словарь синонимов
Книги
- Число имени. Тайны нумерологии. Выход из тела для ленивых. Учебник по экстрасенсорике (количество томов: 3)
- Число имени. Новый взгляд на числа. Нумерология - путь познания (количество томов: 3) , Лоуренс Ширли. Число имени. Тайны нумерологии. Книга Ширли Б. Лоуренс является всесторонним исследованием древней эзотерической системы – нумерологии. Чтобы научиться использовать вибрации чисел для…
