خوارزمية مايسل-ليمر

خوارزمية Meissel -Lehmer (بعد Ernst Meissel و Derrick Henry Lehmer ) هي خوارزمية تحسب القيم الدقيقة لوظيفة العد الأولي . [ 1 ] [ 2 ]

وصف

تعود مشكلة حساب العدد الدقيق للأعداد الأولية الأقل من أو تساوي x ، دون سردها جميعًا، إلى ليجاندر . وقد لاحظ من خلال كتاب "المنخل" لإراتوستينس أن

π(x)-π(x1/2)+1=x-أناx/صأنا+أنا<جx/صأناصج-...{\displaystyle \pi (x)-\pi (x^{1/2})+1=\lfloor x\rfloor -\sum _{i}\lfloor x/p_{i}\rfloor +\sum _{i<j}\lfloor x/p_{i}p_{j}\rfloor -\ldots }

حيث x هي دالة الجزء الصحيح ، والتي تشير إلى أكبر عدد صحيح أصغر من أو يساوي و p i تمتد على جميع الأعداد الأولية x . [ 1 ] [ 2 ]

بما أن حساب صيغة المجموع هذه يصبح أكثر تعقيدًا وإرباكًا مع ازدياد قيمة x ، فقد حاول مايسل تبسيط عملية عدّ الأعداد في غربال إراتوستينس. ولذلك، قدّم هو وليمر دوال غربلة معينة، سيتم تفصيلها أدناه.

الوظائف الرئيسية

ليكن p₁ , p₂ , , pₙ الأعداد الأولية n الأولى . بالنسبة لعدد طبيعي a 1 ، عرّف

φ(x،أ):=|{نx:ص|نص>صأ}|،{\displaystyle \varphi (x,a):=\left|\left\{n\leq x:p|n\implies p>p_{a}\right\}\right|,}

والذي يحسب الأعداد الطبيعية التي لا تزيد عن x والتي جميع عواملها الأولية أكبر من p . عرّف أيضًا للعدد الطبيعي k ،

Pك(x،أ):=|{نx:ن=q1q2qك، مع q1،...،qك>صأ}|،{\displaystyle P_{k}(x,a):=\left|\left\{n\leq x:n=q_{1}q_{2}\cdots q_{k},~{\text{مع}}~q_{1},\ldots ,q_{k}>p_{a}\right\}\right|,}

والتي تحسب الأعداد الطبيعية التي لا تزيد عن x والتي لها k عوامل أولية بالضبط ، وكلها أكبر من p . باستخدام هذه الأعداد، لدينا

φ(x،أ)=ك=0Pك(x،أ)،{\displaystyle \varphi (x,a)=\sum _{k=0}^{\infty }P_{k}(x,a),}

حيث أن المجموع يحتوي فقط على عدد محدود من الحدود غير الصفرية لأن P <sub>k</sub> ( x , a ) = 0 عندما p <sub>ka </sub> > x . باستخدام حقيقة أن P <sub>0</sub> ( x , a ) = 1 و P <sub>1</sub> ( x , a ) = π ( x ) - a ، نحصل على

π(x)=φ(x،أ)+أ-1-ك=2Pك(x،أ)،{\displaystyle \pi (x)=\varphi (x,a)+a-1-\sum _{k=2}^{\infty }P_{k}(x,a),}

وهذا يثبت أنه يمكن حساب π ( x ) عن طريق حساب φ ( x , a ) و Pk ( x , a ) لـ k 2. وهذا ما تفعله خوارزمية مايسل-ليمر.

صيغة P k ( x , a )

بالنسبة لـ k = 2 ، نحصل على الصيغة التالية لـ P k ( x , a ) :

P2(x،أ)=|{ن:نx، ن=صبصج، مع أ<بج}|=ب=أ+1π(x1/2)|{ن:نx، ن=صبصج، مع بجπ(xصب)}|=ب=أ+1π(x1/2)(π(xصب)-(ب-1))=(أ2)-(π(x1/2)2)+ب=أ+1π(x1/2)π(xصب).\displaystyle \begin{aligned}P_{2}(x,a)&=\left|\left\{n:n\leq x,~n=p_{b}p_{c},~{\text{مع}}~a<b\leq c\right\}\right|\\&=\sum _{b=a+1}^{\pi (x^{1/2})}\left|\left\{n:n\leq x,~n=p_{b}p_{c},~{\text{مع}}~b\leq c\leq \pi \left({\frac {x}{p_{b}}}\right)\right\}\right|\\&=\sum _{b=a+1}^{\pi (x^{1/2})}\left(\pi \left({\frac \beginom {x}{p_{b}}}\right)-(b-1)\right)\\&={\binom {a}{2}}-{\binom {\pi (x^{1/2})}{2}}+\sum _{b=a+1}^{\pi (x^{1/2})}\pi \left({\frac {x}{p_{b}}}\right).\end{aligned}}}

بالنسبة لـ k 3 ، يمكن اشتقاق متطابقات P k ( x , a ) بشكل مماثل. [ 1 ]

توسيع φ ( x , a )

مع حالة البداية

φ(x،0)=x،{\displaystyle \varphi (x,0)=\lfloor x\rfloor ,}

وتكرارها

φ(x،أ)=φ(x،أ-1)-φ(xصأ،أ-1)،{\displaystyle \varphi (x,a)=\varphi (x,a-1)-\varphi \left({\frac {x}{p_{a}}},a-1\right),}

يمكن حساب كل قيمة لـ φ ( x , a ) بشكل متكرر.

دمج المصطلحات

الشيء الوحيد المتبقي هو حساب قيمتي φ(x, a) و Pk( x , a ) لـ k2 ، وذلك لقيم معينة من x و a . ويمكن القيام بذلك عن طريق الغربلة المباشرة واستخدام الصيغ المذكورة أعلاه.

تاريخ

وجد مايسل بالفعل أنه بالنسبة لـ k 3 ، فإن P k ( x , a ) = 0 إذا كان a = π ( x 1/3 ) . وقد استخدم المعادلة الناتجة لحساب π ( x ) للقيم الكبيرة لـ x . [ 1 ]

قام مايسل بحساب π ( x ) لقيم x حتى 10 9 ، لكنه أخطأ النتيجة الصحيحة لأكبر قيمة لـ x . [ 1 ]

باستخدام طريقته وجهاز IBM 701 ، تمكن ليمر من حساب القيمة الصحيحة لـ π (10 9 ) وأخطأ في حساب القيمة الصحيحة لـ π (10 10 ) بمقدار 1. [ 1 ]

خوارزمية موسعة

نشر جيفري لاغارياس وفيكتور ميلر وأندرو أودليزكو تطبيقًا للخوارزمية التي تحسب π ( x ) في زمن O ( /³ + ε ) ومساحة O ( /³ + ε ) لأي قيمة ε > 0. [ 2 ] عند تعيين a = π ( ) ، تحتوي شجرة φ ( x , a ) على O ( ) من العقد الورقية. [ 2 ]

تحتاج خوارزمية Meissel-Lehmer الموسعة هذه إلى وقت حساب أقل من الخوارزمية التي طورها Meissel و Lehmer، وخاصة بالنسبة للقيم الكبيرة لـ x .

وقدّم كلٌّ من م. ديليجليس وج. ريفات في عام 1996، وإكس. جوردون في عام 2001 (غير منشور)، تحسينات إضافية على الخوارزمية. [ 3 ] [ 4 ]

مراجع

  1. 1 2 3 4 5 6 ليمر، ديريك هنري (1 أبريل 1958). "حول العدد الدقيق للأعداد الأولية الأقل من حد معين" . مجلة إلينوي للرياضيات . 3 (3): 381-388 . تم الاطلاع عليه في 1 فبراير 2017 .
  2. 1 2 3 4 لاغارياس، جيفري؛ ميلر، فيكتور؛ أودليزكو، أندرو (11 أبريل 1985). "الحوسبةπ(x){\displaystyle \pi (x)}: طريقة Meissel – Lehmer " (PDF) . رياضيات الحساب . 44 (170): 537–560 . دوى : 10.1090 / S0025-5718-1985-0777285-5 . تم الاسترجاع في 13 سبتمبر، 2016 .
  3. ديليغليز، مارك؛ ريفات، جويل (15 يناير 1996). "الحوسبةπ(x){\displaystyle \pi (x)}: طريقة Meissel، Lehmer، Lagarias، Miller، Odlyzko " . الرياضيات الحسابية . 65 (213): 235–245 . دوى : 10.1090 / S0025-5718-96-00674-6 .
  4. ^ أوليفيرا إي سيلفا ، توماس (1 مارس 2006). "الحوسبةπ(x){\displaystyle \pi (x)}"الطريقة التوافقية" (ملف PDF) . مجلة ديتوا . 4 (6): 759-768 . تاريخ الاسترجاع: 14 مارس 2023 .