خوارزمية بورواين

ابتكر جوناثان وبيتر بورواين خوارزمية بورواين لحساب قيمة1/π{\displaystyle 1/\pi }يمكن العثور على هذه الخوارزمية وغيرها في كتاب "باي والخوارزمية الحسابية متعددة الأعداد - دراسة في نظرية الأعداد التحليلية والتعقيد الحسابي" . [ 1 ]

سلسلة رامانوجان-ساتو

هذان مثالان على متسلسلة رامانوجان-ساتو . وتستخدم خوارزمية تشودنوفسكي ذات الصلة مُمَيِّزًا برقم فئة 1.

الصف رقم 2 (1989)

ابدأ بضبط [ 2 ]

أ=21217571091261+1657145277365ب=1377398089267261+107578229802750ج=(5280(236674+3030361))3{\displaystyle {\begin{aligned}A&=212175710912{\sqrt {61}}+1657145277365\\B&=13773980892672{\sqrt {61}}+107578229802750\\C&=\left(5280\left(236674+30303{\sqrt {61}}\right)\right)^{3}\end{aligned}}}

ثم

1π=12ن=0(-1)ن(6ن)!(أ+نب)(ن!)3(3ن)!جن+12{\displaystyle {\frac {1}{\pi }}=12\sum _{n=0}^{\infty }{\frac {(-1)^{n}(6n)!\,(A+nB)}{(n!)^{3}(3n)!\,C^{n+{\frac {1}{2}}}}}}

ينتج عن كل حد إضافي من المجموع الجزئي حوالي 25 رقماً.

رقم الفصل 4 (1993)

ابدأ بضبط [ 3 ]

أ=63365028312971999585426220+283377021408008420468256005+3845(10891728551171178200467436212395209160385656017+48709290865788102250773385345416887213512550405)12ب=7849910453496627210289749000+35105866782609320289656064005+25159683110(6260208323789001636993322654444020882161+27996502730604442965772068907188251902355)12ج=-214772995063512240-960494033386480325-12965(10985234579463550323713318473+49127462536923627546073959125)12{\displaystyle {\begin{aligned}A={}&63365028312971999585426220\\&{}+28337702140800842046825600{\sqrt {5}}\\&{}+384{\sqrt {5}}{\big (}10891728551171178200467436212395209160385656017\\&{}+\left.4870929086578810225077338534541688721351255040{\sqrt {5}}\right)^{\frac {1}{2}}\\B={}&7849910453496627210289749000\\&{}+3510586678260932028965606400{\sqrt {5}}\\&{}+2515968{\sqrt {3110}}{\big (}6260208323789001636993322654444020882161\\&{}+\left.2799650273060444296577206890718825190235{\sqrt {5}}\right)^{\frac {1}{2}}\\C={}&-214772995063512240\\&{}-96049403338648032{\sqrt {5}}\\&{}-1296{\sqrt {5}}{\big (}10985234579463550323713318473\\&{}+\left.4912746253692362754607395912{\sqrt {5}}\right)^{\frac {1}{2}}\end{aligned}}}

ثم

-ج3π=ن=0(6ن)!(3ن)!(ن!)3أ+نبج3ن{\displaystyle {\frac {\sqrt {-C^{3}}}{\pi }}=\sum _{n=0}^{\infty }{{\frac {(6n)!}{(3n)!(n!)^{3}}}{\frac {A+nB}{C^{3n}}}}}

ينتج عن كل حد إضافي من المتسلسلة ما يقرب من 50 رقمًا.

الخوارزميات التكرارية

التقارب التربيعي (1984)

ابدأ بضبط [ 4 ]

أ0=2ب0=0ص0=2+2{\displaystyle {\begin{aligned}a_{0}&={\sqrt {2}}\\b_{0}&=0\\p_{0}&=2+{\sqrt {2}}\end{aligned}}}

ثم كرر العملية

أن+1=أن+1أن2بن+1=(1+بن)أنأن+بنصن+1=(1+أن+1)صنبن+11+بن+1{\displaystyle {\begin{aligned}a_{n+1}&={\frac {{\sqrt {a_{n}}}+{\frac {1}{\sqrt {a_{n}}}}}{2}}\\b_{n+1}&={\frac {(1+b_{n}){\sqrt {a_{n}}}}{a_{n}+b_{n}}}\\p_{n+1}&={\frac {(1+a_{n+1})\,p_{n}b_{n+1}}{1+b_{n+1}}}\end{aligned}}}

ثم يتقارب p k تربيعيًا إلى π ؛ أي أن كل تكرار يُضاعف تقريبًا عدد الأرقام الصحيحة. الخوارزمية ليست ذاتية التصحيح؛ يجب تنفيذ كل تكرار بالعدد المطلوب من الأرقام الصحيحة للحصول على النتيجة النهائية π .

التقارب التكعيبي (1991)

ابدأ بالضبط

أ0=13s0=3-12{\displaystyle {\begin{aligned}a_{0}&={\frac {1}{3}}\\s_{0}&={\frac {{\sqrt {3}}-1}{2}}\end{aligned}}}

ثم كرر العملية

رك+1=31+2(1-sك3)13sك+1=رك+1-12أك+1=رك+12أك-3ك(رك+12-1){\displaystyle {\begin{aligned}r_{k+1}&={\frac {3}{1+2\left(1-s_{k}^{3}\right)^{\frac {1}{3}}}}\\s_{k+1}&={\frac {r_{k+1}-1}{2}}\\a_{k+1}&=r_{k+1}^{2}a_{k}-3^{k}\left(r_{k+1}^{2}-1\right)\end{aligned}}}

ثم يتقارب k بشكل مكعب إلى 1 / π ؛ أي أن كل تكرار يزيد عدد الأرقام الصحيحة بمقدار ثلاثة أضعاف تقريبًا .

التقارب الرباعي (1985)

ابدأ بضبط [ 5 ]

أ0=2(2-1)2y0=2-1{\displaystyle {\begin{aligned}a_{0}&=2\left({\sqrt {2}}-1\right)^{2}\\y_{0}&={\sqrt {2}}-1\end{aligned}}}

ثم كرر العملية

yك+1=1-(1-yك4)141+(1-yك4)14أك+1=أك(1+yك+1)4-22ك+3yك+1(1+yك+1+yك+12){\displaystyle {\begin{aligned}y_{k+1}&={\frac {1-\left(1-y_{k}^{4}\right)^{\frac {1}{4}}}{1+\left(1-y_{k}^{4}\right)^{\frac {1}{4}}}}\\a_{k+1}&=a_{k}\left(1+y_{k+1}\right)^{4}-2^{2k+3}y_{k+1}\left(1+y_{k+1}+y_{k+1}^{2}\right)\end{aligned}}}

ثم يتقارب k بشكل رباعي مع 1 / π ؛ أي أن كل تكرار يزيد عدد الأرقام الصحيحة بمقدار أربعة أضعاف تقريبًا. الخوارزمية ليست ذاتية التصحيح؛ يجب تنفيذ كل تكرار بالعدد المطلوب من الأرقام الصحيحة للحصول على النتيجة النهائية لـ π .

تُعادل دورة واحدة من هذه الخوارزمية دورتين من خوارزمية غاوس-ليجندر . ويمكن الاطلاع على برهان هذه الخوارزميات هنا: [ 6 ]

التقارب الخماسي

ابدأ بالضبط

أ0=12s0=5(5-2)=5ϕ3{\displaystyle {\begin{aligned}a_{0}&={\frac {1}{2}}\\s_{0}&=5\left({\sqrt {5}}-2\right)={\frac {5}{\phi ^{3}}}\end{aligned}}}

أينϕ=1+52{\displaystyle \phi ={\tfrac {1+{\sqrt {5}}}{2}}}هي النسبة الذهبية . ثم كرر العملية.

xن+1=5sن-1yن+1=(xن+1-1)2+7zن+1=(12xن+1(yن+1+yن+12-4xن+13))15أن+1=sن2أن-5ن(sن2-52+sن(sن2-2sن+5))sن+1=25(zن+1+xن+1zن+1+1)2sن{\displaystyle {\begin{aligned}x_{n+1}&={\frac {5}{s_{n}}}-1\\y_{n+1}&=\left(x_{n+1}-1\right)^{2}+7\\z_{n+1}&=\left({\frac {1}{2}}x_{n+1}\left(y_{n+1}+{\sqrt {y_{n+1}^{2}-4x_{n+1}^{3}}}\right)\right)^{\frac {1}{5}}\\a_{n+1}&=s_{n}^{2}a_{n}-5^{n}\left({\frac {s_{n}^{2}-5}{2}}+{\sqrt {s_{n}\left(s_{n}^{2}-2s_{n}+5\right)}}\right)\\s_{n+1}&={\frac {25}{\left(z_{n+1}+{\frac {x_{n+1}}{z_{n+1}}}+1\right)^{2}s_{n}}}\end{aligned}}}

ثم يتقارب k بشكل خماسي إلى 1 / π ( أي أن كل تكرار يضاعف عدد الأرقام الصحيحة تقريبًا خمس مرات)، ويتحقق الشرط التالي :

0<أن-1π<165نهـ-5نπ{\displaystyle 0<a_{n}-{\frac {1}{\pi }}<16\cdot 5^{n}\cdot e^{-5^{n}}\pi \,\!}

التقارب غير المتناظر

ابدأ بالضبط

أ0=13ر0=3-12s0=(1-ر03)13{\displaystyle {\begin{aligned}a_{0}&={\frac {1}{3}}\\r_{0}&={\frac {{\sqrt {3}}-1}{2}}\\s_{0}&=\left(1-r_{0}^{3}\right)^{\frac {1}{3}}\end{aligned}}}

ثم كرر العملية

تن+1=1+2رنuن+1=(9رن(1+رن+رن2))13vن+1=تن+12+تن+1uن+1+uن+12wن+1=27(1+sن+sن2)vن+1أن+1=wن+1أن+32ن-1(1-wن+1)sن+1=(1-رن)3(تن+1+2uن+1)vن+1رن+1=(1-sن+13)13{\displaystyle {\begin{aligned}t_{n+1}&=1+2r_{n}\\u_{n+1}&=\left(9r_{n}\left(1+r_{n}+r_{n}^{2}\right)\right)^{\frac {1}{3}}\\v_{n+1}&=t_{n+1}^{2}+t_{n+1}u_{n+1}+u_{n+1}^{2}\\w_{n+1}&={\frac {27\left(1+s_{n}+s_{n}^{2}\right)}{v_{n+1}}}\\a_{n+1}&=w_{n+1}a_{n}+3^{2n-1}\left(1-w_{n+1}\right)\\s_{n+1}&={\frac {\left(1-r_{n}\right)^{3}}{\left(t_{n+1}+2u_{n+1}\right)v_{n+1}}}\\r_{n+1}&=\left(1-s_{n+1}^{3}\right)^{\frac {1}{3}}\end{aligned}}}

ثم يتقارب k بشكل غير منتظم إلى 1 / π ؛ أي أن كل تكرار يضرب عدد الأرقام الصحيحة تقريبًا في تسعة. [ 7 ]

انظر أيضاً

مراجع

  1. جوناثان م. بورواين، بيتر ب. بورواين، باي والنموذج التراكمي المتوسط ​​- دراسة في نظرية الأعداد التحليلية والتعقيد الحسابي ، وايلي، نيويورك، 1987. العديد من نتائجهم متاحة في: يورغ أرندت، كريستوف هانيل، باي بلا حدود، سبرينغر، برلين، 2001، ISBN 3-540-66572-2
  2. بيلي، ديفيد هـ (2023-04-01). "بيتر بورواين: عالم رياضيات صاحب رؤية". إشعارات الجمعية الرياضية الأمريكية . 70 (4): 610-613 . doi : 10.1090/noti2675 . ISSN 0002-9920 . 
  3. بورواين، ج.م.؛ بورواين، ب.ب. (1993). "متسلسلة رامانوجان من النوع الثالث للفئة 1/π" . مجلة الرياضيات الحسابية والتطبيقية . 46 ( 1-2 ): 281-290 . doi : 10.1016/0377-0427(93)90302-R .
  4. ^ أرندت، يورج. هينيل ، كريستوف (1998). π أطلق العنان . سبرينغر-فيرلاغ. ص. 236. ردمك  3-540-66572-2.
  5. ماك، رونالد (2003). دليل مبرمجي جافا للحساب العددي . بيرسون التعليمية. ص 353. ISBN  0-13-046041-9.
  6. ميلا، لورنز (2019)، برهان سهل لثلاث خوارزميات π تكرارية ، arXiv : 1907.04110
  7. ^ هنريك فيسترمارك (4 نوفمبر 2016). "التنفيذ العملي لخوارزميات π" (PDF) . تم الاسترجاع في 29 نوفمبر 2020 .