خوارزمية بيرلكامب-رابين

إلوين ر. بيرلكامب في مؤتمر حول نظرية الألعاب التوافقية في محطة بانف الدولية للأبحاث

في نظرية الأعداد ، تُعد خوارزمية بيرلكامب لإيجاد الجذور ، والتي تُسمى أيضًا خوارزمية بيرلكامب-رابين ، الطريقة الاحتمالية لإيجاد جذور كثيرات الحدود على الحقل.Fص{\displaystyle \mathbb {F} _{p}}معص{\displaystyle p}العناصر. اكتشف إلوين بيرلكامب هذه الطريقة عام 1970 [ 1 ] كأداة مساعدة لخوارزمية تحليل كثيرات الحدود على الحقول المنتهية. ثم عدّل رابين الخوارزمية لاحقًا لتشمل الحقول المنتهية العشوائية عام 1979. [ 2 ] كما اكتشف باحثون آخرون هذه الطريقة بشكل مستقل قبل بيرلكامب. [ 3 ]

تاريخ

اقترح إلوين بيرلكامب هذه الطريقة في بحثه عام 1970 [ 1 ] حول تحليل كثيرات الحدود على الحقول المنتهية. افتقر عمله الأصلي إلى برهان رسمي على صحته [ 2 ] ، ثم قام مايكل رابين بتنقيحه وتعديله لاحقًا ليشمل الحقول المنتهية المختلفة . [ 2 ] وفي عام 1986، اقترح رينيه بيرالتا خوارزمية مشابهة [ 4 ] لإيجاد الجذور التربيعية فيFص{\displaystyle \mathbb {F} _{p}}[ 5 ] في عام 2000 تم تعميم طريقة بيرالتا للمعادلات التكعيبية . [ 6 ]

بيان المشكلة

يتركص{\displaystyle p}ليكن عددًا أوليًا فرديًا. لنعتبر كثيرة الحدودو(x)=أ0+أ1x++أنxن{\textstyle f(x)=a_{0}+a_{1}x+\cdots +a_{n}x^{n}}في الملعبFصZ/صZ{\displaystyle \mathbb {F} _{p}\simeq \mathbb {Z} /p\mathbb {Z} }الباقي moduloص{\displaystyle p}ينبغي أن تجد الخوارزمية جميعλ{\displaystyle \lambda }فيFص{\displaystyle \mathbb {F} _{p}}بحيثو(λ)=0{\textstyle f(\lambda )=0}فيFص{\displaystyle \mathbb {F} _{p}}[ 2 ] [ 7 ]

الخوارزمية

التوزيع العشوائي

يتركو(x)=(x-λ1)(x-λ2)(x-λن){\textstyle f(x)=(x-\lambda _{1})(x-\lambda _{2})\cdots (x-\lambda _{n})}إن إيجاد جميع جذور هذه المعادلة متعددة الحدود يكافئ إيجاد تحليلها إلى عوامل خطية. ولإيجاد هذا التحليل، يكفي تقسيم المعادلة إلى أي قاسمين غير تافهين وتحليلهما بشكل متكرر. وللقيام بذلك، لننظر إلى المعادلة متعددة الحدود التالية:وz(x)=و(x-z)=(x-λ1-z)(x-λ2-z)(x-λن-z){\textstyle f_{z}(x)=f(xz)=(x-\lambda _{1}-z)(x-\lambda _{2}-z)\cdots (x-\lambda _{n}-z)}أينz{\displaystyle z} هو عنصر منFص{\displaystyle \mathbb {F} _{p}}إذا أمكن تمثيل هذه متعددة الحدود كحاصل ضربوz(x)=ص0(x)ص1(x){\displaystyle f_{z}(x)=p_{0}(x)p_{1}(x)}ثم من حيث كثير الحدود الأولي، فهذا يعني أنو(x)=ص0(x+z)ص1(x+z){\displaystyle f(x)=p_{0}(x+z)p_{1}(x+z)}، مما يوفر التحليل اللازم لـو(x){\displaystyle f(x)}[ 1 ] [ 7 ]

تصنيفFص{\displaystyle \mathbb {F} _{p}}عناصر

بسبب معيار أويلر ، لكل حد أحادي(x-λ){\displaystyle (x-\lambda )}تتحقق إحدى الخصائص التالية فقط: [ 1 ]

  1. الحدّ الواحد يساويx{\displaystyle x}لوλ=0{\displaystyle \lambda =0}،
  2. يقسم الحد الأحاديز0(x)=(x(ص-1)/2-1){\textstyle g_{0}(x)=(x^{(p-1)/2}-1)}لوλ{\displaystyle \lambda } هو الباقي التربيعي moduloص{\displaystyle p}،
  3. يقسم الحد الأحاديز1(x)=(x(ص-1)/2+1){\textstyle g_{1}(x)=(x^{(p-1)/2}+1)}لوλ{\displaystyle \lambda } هو باقي تربيعي moduloص{\displaystyle p}.

وبالتالي إذاوz(x){\displaystyle f_{z}(x)}لا يقبل القسمة علىx{\displaystyle x}والتي يمكن التحقق منها بشكل منفصل، ثموz(x){\displaystyle f_{z}(x)}يساوي حاصل ضرب القواسم المشتركة الكبرىالقاسم المشترك الأكبر(وz(x)؛ز0(x)){\displaystyle \gcd(f_{z}(x);g_{0}(x))}والقاسم المشترك الأكبر(وz(x)؛ز1(x)){\displaystyle \gcd(f_{z}(x);g_{1}(x))}[ 7 ]

أسلوب بيرلكامب

تؤدي الخاصية المذكورة أعلاه إلى الخوارزمية التالية: [ 1 ]

  1. احسب معاملاتوz(x)=و(x-z){\displaystyle f_{z}(x)=f(xz)}،
  2. احسب باقي قسمةx،x2،x22،x23،x24،...،x2سجل2ص{\textstyle x,x^{2},x^{2^{2}},x^{2^{3}},x^{2^{4}},\ldots ,x^{2^{\lfloor \log _{2}p\rfloor }}}moduloوz(x){\displaystyle f_{z}(x)}بتربيع كثير الحدود الحالي وأخذ الباقي moduloوz(x){\displaystyle f_{z}(x)}،
  3. باستخدام عملية الأسس بالتربيع وكثيرات الحدود المحسوبة في الخطوات السابقة، احسب باقي قسمةx(ص-1)/2{\textstyle x^{(p-1)/2}}moduloوz(x){\textstyle f_{z}(x)}،
  4. لوx(ص-1)/2±1(تعديلوz(x)){\textstyle x^{(p-1)/2}\not \equiv \pm 1{\pmod {f_{z}(x)}}}ثمالقاسم المشترك الأكبر{\displaystyle \gcd }تُقدّم الأمثلة المذكورة أدناه تحليلًا غير تافه لـوz(x){\displaystyle f_{z}(x)}،
  5. وإلا فإن جميع جذوروz(x){\displaystyle f_{z}(x)}إما أن تكون بقايا أو غير بقايا في آن واحد، ويتعين على المرء اختيار الآخر.z{\displaystyle z}.

لوو(x){\displaystyle f(x)}يقبل القسمة على كثير حدود أولي غير خطيز(x){\displaystyle g(x)}زيادةFص{\displaystyle \mathbb {F} _{p}}ثم عند الحسابالقاسم المشترك الأكبر{\displaystyle \gcd }معز0(x){\displaystyle g_{0}(x)}وز1(x){\displaystyle g_{1}(x)}سيحصل المرء على تحليل غير تافه لـوz(x)/زz(x){\displaystyle f_{z}(x)/g_{z}(x)}وبالتالي، تسمح هذه الخوارزمية بإيجاد جميع جذور كثيرات الحدود العشوائية علىFص{\displaystyle \mathbb {F} _{p}}.

الجذر التربيعي المعياري

ضع في اعتبارك المعادلةx2أ(تعديلص){\textstyle x^{2}\equiv a{\pmod {p}}}يحتوي على عناصرβ{\displaystyle \beta }و-β{\displaystyle -\beta }كجذورها. حل هذه المعادلة يكافئ تحليل كثير الحدودو(x)=x2-أ=(x-β)(x+β){\textstyle f(x)=x^{2}-a=(x-\beta )(x+\beta )}زيادةFص{\displaystyle \mathbb {F} _{p}}في هذه الحالة المحددة، يكفي حساب فقطالقاسم المشترك الأكبر(وz(x)؛ز0(x)){\displaystyle \gcd(f_{z}(x);g_{0}(x))}بالنسبة لهذه المعادلة متعددة الحدود، ستتحقق إحدى الخصائص التالية فقط:

  1. القاسم المشترك الأكبر يساوي1{\displaystyle 1}وهذا يعني أنz+β{\displaystyle z+\beta }وz-β{\displaystyle z-\beta }كلاهما عبارة عن بواقي غير تربيعية،
  2. القاسم المشترك الأكبر يساويوz(x){\displaystyle f_{z}(x)}وهذا يعني أن كلا العددين عبارة عن بواقي تربيعية،
  3. القاسم المشترك الأكبر يساوي(x-ت){\displaystyle (xt)}وهذا يعني أن أحد هذه الأرقام بالضبط هو الباقي التربيعي.

في الحالة الثالثة، يكون القاسم المشترك الأكبر مساوياً لأحد القيمتين التاليتين:(x-z-β){\displaystyle (xz-\beta )}أو(x-z+β){\displaystyle (x-z+\beta )}يسمح ذلك بكتابة الحل على النحو التالي:β=(ت-z)(تعديلص){\textstyle \beta =(tz){\pmod {p}}}[ 1 ]

مثال

لنفترض أننا بحاجة إلى حل المعادلةx25(تعديل11){\textstyle x^{2}\equiv 5{\pmod {11}}}لذلك نحتاج إلى تحليل المقدار إلى عوامله الأولية.و(x)=x2-5=(x-β)(x+β){\displaystyle f(x)=x^{2}-5=(x-\beta )(x+\beta )}ضع في اعتبارك بعض القيم المحتملة لـz{\displaystyle z}:

  1. يتركz=3{\displaystyle z=3}. ثموz(x)=(x-3)2-5=x2-6x+4{\displaystyle f_{z}(x)=(x-3)^{2}-5=x^{2}-6x+4}، هكذاالقاسم المشترك الأكبر(x2-6x+4؛x5-1)=1{\displaystyle \gcd(x^{2}-6x+4;x^{5}-1)=1}كلا الرقمين3±β{\displaystyle 3\pm \beta }هي بواقي غير تربيعية، لذلك نحتاج إلى أخذ شيء آخرz{\displaystyle z}.
  1. يتركz=2{\displaystyle z=2}. ثموz(x)=(x-2)2-5=x2-4x-1{\displaystyle f_{z}(x)=(x-2)^{2}-5=x^{2}-4x-1}، هكذاالقاسم المشترك الأكبر(x2-4x-1؛x5-1)x-9(تعديل11){\textstyle \gcd(x^{2}-4x-1;x^{5}-1)\equiv x-9{\pmod {11}}}ويترتب على ذلكx-9=x-2-β{\textstyle x-9=x-2-\beta }، لذاβ7(تعديل11){\displaystyle \beta \equiv 7{\pmod {11}}}و-β-74(تعديل11){\textstyle -\beta \equiv -7\equiv 4{\pmod {11}}}.

أظهر فحص يدوي أن ذلك بالفعل72495(تعديل11){\textstyle 7^{2}\equiv 49\equiv 5{\pmod {11}}}و42165(تعديل11){\textstyle 4^{2}\equiv 16\equiv 5{\pmod {11}}}.

إثبات صحة النتائج

تجد الخوارزمية تحليل العوامل لـوz(x){\displaystyle f_{z}(x)}في جميع الحالات باستثناء الحالات التي تكون فيها جميع الأرقامz+λ1،z+λ2،...،z+λن{\displaystyle z+\lambda _{1},z+\lambda _{2},\ldots ,z+\lambda _{n}}تكون البقايا التربيعية أو غير البقايا في آن واحد. وفقًا لنظرية القطع الدائري ، [ 8 ] فإن احتمال حدوث مثل هذا الحدث في الحالة التيλ1،...،λن{\displaystyle \lambda _{1},\ldots ,\lambda _{n}}هل جميعها بقايا أو غير بقايا في آن واحد (أي عندماz=0{\displaystyle z=0}(سيفشل) يمكن تقديره على النحو التالي2-ك{\displaystyle 2^{-k}}أينك{\displaystyle k} يمثل عدد القيم المميزة فيλ1،...،λن{\displaystyle \lambda _{1},\ldots ,\lambda _{n}}[ 1 ] وبهذه الطريقة حتى في أسوأ الحالاتك=1{\displaystyle k=1}وو(x)=(x-λ)ن{\displaystyle f(x)=(x-\lambda )^{n}}يمكن تقدير احتمال الخطأ على النحو التالي:1/2{\displaystyle 1/2}أما بالنسبة لحالة الجذر التربيعي المعياري، فإن احتمال الخطأ يكون على الأكثر1/4{\displaystyle 1/4}.

تعقيد

ليكن كثير الحدود من الدرجةن{\displaystyle n}نستنتج تعقيد الخوارزمية على النحو التالي:

  1. بسبب نظرية ذات الحدين(x-z)ك=أنا=0ك(كأنا)(-z)ك-أناxأنا{\textstyle (xz)^{k}=\sum \limits _{i=0}^{k}{\binom {k}{i}}(-z)^{ki}x^{i}}، قد ننتقل منو(x){\displaystyle f(x)}لو(x-z){\displaystyle f(xz)}فييا(ن2){\displaystyle O(n^{2})}وقت.
  2. يمكن إجراء عملية ضرب كثيرات الحدود وأخذ باقي قسمة كثيرة حدود على أخرى فييا(ن2){\textstyle O(n^{2})}وبالتالي حسابx2كتعديلوz(x){\textstyle x^{2^{k}}{\bmod {f}}_{z}(x)}يتم ذلك فييا(ن2سجلص){\textstyle O(n^{2}\log p)}.
  3. يعمل الأس الثنائي فييا(ن2سجلص){\displaystyle O(n^{2}\log p)}.
  4. أخذالقاسم المشترك الأكبر{\displaystyle \gcd }يعمل حل كثيرتي حدود باستخدام خوارزمية إقليدية فييا(ن2){\displaystyle O(n^{2})}.

وبالتالي يمكن إتمام الإجراء بأكمله فييا(ن2سجلص){\displaystyle O(n^{2}\log p)}باستخدام تحويل فورييه السريع وخوارزمية نصف القاسم المشترك الأكبر، [ 9 ] يمكن تحسين تعقيد الخوارزمية إلىيا(نسجلنسجلصن){\displaystyle O(n\log n\log pn)}في حالة الجذر التربيعي المعياري، تكون الدرجة هين=2{\displaystyle n=2}وبالتالي، فإن التعقيد الكلي للخوارزمية في هذه الحالة يكون محدودًا بـيا(سجلص){\displaystyle O(\log p)}لكل تكرار. [ 7 ]

مراجع

  1. 1 2 3 4 5 6 7 بيرلكامب، إي آر (1970). "تحليل كثيرات الحدود على الحقول المنتهية الكبيرة" . رياضيات الحساب . 24 (111): 713-735 . doi : 10.1090/S0025-5718-1970-0276200-X . ISSN 0025-5718 . 
  2. 1 2 3 4 م. رابين (1980). “الخوارزميات الاحتمالية في المجالات المحدودة”. مجلة SIAM للحوسبة . 9 (2): 273-280 . سيتيسيركس 10.1.1.17.5653 . دوى : 10.1137/0209024 . ISSN 0097-5397 .  
  3. دونالد إي كنوث (1998). فن برمجة الحاسوب. المجلد 2. أديسون-ويسلي. ISBN 978-0201896848. OCLC 900627019 . 
  4. تسز-وو سزي (2011). "حول إيجاد الجذور التربيعية بدون بواقي غير تربيعية على الحقول المنتهية". رياضيات الحساب . 80 (275): 1797-1811 . arXiv : 0812.2591 . doi : 10.1090 /s0025-5718-2011-02419-1 . ISSN 0025-5718 . S2CID 10249895 .  
  5. ر. بيرالتا (نوفمبر 1986). "خوارزمية احتمالية بسيطة وسريعة لحساب الجذور التربيعية بتردد عدد أولي (مراسلات)". معاملات IEEE في نظرية المعلومات . 32 (6): 846-847 . doi : 10.1109/TIT.1986.1057236 . ISSN 0018-9448 . 
  6. سي بادرو، جي سايز (أغسطس 2002). "إيجاد الجذور التكعيبية في Zm". رسائل الرياضيات التطبيقية . 15 (6): 703-708 . doi : 10.1016/s0893-9659(02)00031-9 . ISSN 0893-9659 . 
  7. 1 2 3 4 ألفريد ج. مينيز، إيان ف. بليك، شوهونغ غاو، رونالد س. مولين، سكوت أ. فانستون (1993). تطبيقات الحقول المنتهية . سلسلة سبرينغر الدولية في الهندسة وعلوم الحاسوب. سبرينغر الولايات المتحدة. ISBN 9780792392828.{{cite book}}: صيانة CS1: أسماء متعددة: قائمة المؤلفين ( رابط )
  8. مارشال هول (1998). نظرية التوافيق . جون وايلي وأولاده. ISBN 9780471315186.
  9. أهو، ألفريد ف. (1974). تصميم وتحليل خوارزميات الحاسوب . شركة أديسون-ويسلي للنشر. ISBN 0201000296.