Prime Numbers on kenji.blog

فرضية ريمان وتوزيع الأعداد الأولية: علاقة عميقة مع تقنيات التشفير الحديثة

Fri, 11 Sep 2026 16:00:00 +0900

1. مقدمة: غموض الأعداد الأولية في الكون وفرضية ريمان

“الأعداد الأولية (Prime Numbers)” هي الأعداد الطبيعية التي لا تقبل القسمة إلا على 1 وعلى نفسها، وتُعرف أيضاً بـ “ذرات” عالم الرياضيات. تبدو هذه السلسلة من الأعداد 2, 3, 5, 7, 11, 13… عشوائية وبدون نظام واضح للوهلة الأولى. منذ أن أثبت عالم الرياضيات اليوناني القديم إقليدس أن “الأعداد الأولية لا نهائية”، حاول عدد لا يحصى من علماء الرياضيات كشف القواعد الخفية وراء تسلسل هذه الأعداد.

أقرب محاولة لحل لغز الأعداد الأولية كانت “فرضية ريمان (Riemann Hypothesis)” التي اقترحها عالم الرياضيات الألماني برنارد ريمان (Bernhard Riemann) في عام 1859. تعتبر فرضية ريمان واحدة من أهم المسائل المستعصية في الرياضيات الحديثة، وهي إحدى مسائل جائزة الألفية التي حددها معهد كلاي للرياضيات بجائزة قدرها مليون دولار.

للوهلة الأولى، قد يبدو لغز التوزيع للأعداد الأولية في الرياضيات البحتة غير ذي صلة بحياتنا اليومية. ومع ذلك، فإن أمن الإنترنت الذي يدعم البنية التحتية للمجتمع الحديث، وخاصة تقنيات التشفير الحديثة مثل تشفير RSA وتشفير المنحنيات الإهليلجية (ECC)، يعتمد بشدة على خصائص الأعداد الأولية الضخمة.

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

2. مبرهنة الأعداد الأولية وتوزيع الأعداد الأولية: اكتشاف غاوس

لفهم كيفية توزيع الأعداد الأولية، فكر علماء الرياضيات في دالة عد الأعداد الأولية (Prime-counting function) $\pi(x)$ التي تعبر عن “عدد الأعداد الأولية الموجودة أقل من أو تساوي عدداً معيناً $x$”.

على سبيل المثال:

$\pi(10) = 4$ (2, 3, 5, 7)
$\pi(100) = 25$
$\pi(1000) = 168$

اكتشف عالم الرياضيات العبقري كارل فريدريش غاوس (Carl Friedrich Gauss) في سن الخامسة عشرة، بعد حساب جداول ضخمة للأعداد الأولية، أن تكرار ظهور الأعداد الأولية يتناقص بشكل يتناسب عكسياً مع اللوغاريتم الطبيعي $\ln x$. بمعنى آخر، توقع أن احتمال العثور على عدد أولي بالقرب من عدد معين $x$ هو حوالي $\frac{1}{\ln x}$.

تم التعبير عن هذا باستخدام التكامل في التكامل اللوغاريتمي (Logarithmic integral) $\text{Li}(x)$:

$$ \text{Li}(x) = \int_{2}^{x} \frac{dt}{\ln t} $$

تم إثبات توقع غاوس لاحقاً وبشكل مستقل من قبل جاك هادامار وشارل جان دو لا فالي بوسان في عام 1896، وتأسست كـ مبرهنة الأعداد الأولية (Prime Number Theorem, PNT):

$$ \lim_{x \to \infty} \frac{\pi(x)}{\text{Li}(x)} = 1 $$

أو يمكن التعبير عنها تقريبياً على النحو التالي:

$$ \pi(x) \sim \frac{x}{\ln x} $$

بفضل هذه المبرهنة، أصبح من الواضح أن الأعداد الأولية، عند النظر إليها على نطاق واسع (ماكروسكوبي)، تمتلك توزيعاً سلساً للغاية ويمكن التنبؤ به. ومع ذلك، عند النظر إليها على نطاق ضيق (ميكروسكوبي)، هناك دائماً “خطأ” أو “تذبذب” بين $\pi(x)$ و $\text{Li}(x)$. الهوية الحقيقية لهذا التذبذب هي اللغز الأكبر الذي تحاول فرضية ريمان حله.

3. دالة زيتا لريمان وحاصل ضرب أويلر

أقوى سلاح في تحليل توزيع الأعداد الأولية هو دالة زيتا لريمان (Riemann Zeta Function). في الأصل، كانت متسلسلة لانهائية تم تعريفها من قبل ليونهارت أويلر (Leonhard Euler) للأعداد الحقيقية $s > 1$:

$$ \zeta(s) = \sum_{n=1}^\infty \frac{1}{n^s} = 1 + \frac{1}{2^s} + \frac{1}{3^s} + \frac{1}{4^s} + \dots $$

أحد أعظم إنجازات أويلر هو إثباته أن هذه المتسلسلة اللانهائية يمكن التعبير عنها كحاصل ضرب لانهائي يشمل جميع الأعداد الأولية $p$. هذا هو حاصل ضرب أويلر (Euler Product Formula):

$$ \zeta(s) = \prod_{p \text{ prime}} \frac{1}{1 - p^{-s}} = \left( \frac{1}{1 - 2^{-s}} \right) \left( \frac{1}{1 - 3^{-s}} \right) \left( \frac{1}{1 - 5^{-s}} \right) \dots $$

الفهم الحدسي لهذا الإثبات هو أنه إذا قمنا بنشر كل حد في الجانب الأيمن كمتسلسلة هندسية وضربناها معاً، فبفضل المبرهنة الأساسية في الحساب (كل عدد طبيعي يمكن التعبير عنه بشكل فريد كحاصل ضرب أعداد أولية)، يتم إعادة بناء مجموع مقلوبات الأعداد الطبيعية في الجانب الأيسر بالكامل.

هذه المعادلة الرياضية الواحدة أصبحت جسراً يربط بين التحليل الرياضي (المتسلسلات اللانهائية والدوال المستمرة) ونظرية الأعداد (الأعداد الأولية والأعداد المنفصلة). دراسة دالة زيتا مرادفة لدراسة توزيع الأعداد الأولية.

4. الامتداد التحليلي والتوسيع إلى المستوى المركب

تكمن عبقرية ريمان في أنه وسع المتغير $s$ في الدالة $\zeta(s)$، الذي كان أويلر يعتبره عدداً حقيقياً فقط، ليصبح عدداً مركباً $s = \sigma + it$ ($\sigma$ هو الجزء الحقيقي و $t$ هو الجزء التخيلي).

المتسلسلة اللانهائية الأصلية تتقارب فقط عندما يكون $\sigma > 1$، ولكن باستخدام تقنية “الامتداد التحليلي (Analytic Continuation)"، وسع ريمان التعريف بحيث أصبح للدالة $\zeta(s)$ معنى في كل المستوى المركب باستثناء القطب عند $s = 1$.

وعلاوة على ذلك، استنتج معادلة دالية (Functional equation) جميلة تحققها دالة زيتا:

$$ \zeta(s) = 2^s \pi^{s-1} \sin\left(\frac{\pi s}{2}\right) \Gamma(1-s) \zeta(1-s) $$

حيث $\Gamma(x)$ هي دالة غاما. من خلال هذه المعادلة، يمكننا معرفة خصائص النصف الأيسر من المستوى بناءً على خصائص النصف الأيمن.

الأصفار (Zeros of the Zeta Function)

تُسمى الأعداد المركبة $s$ التي تجعل قيمة دالة زيتا تساوي 0 بـ “الأصفار”. من المعادلة الدالية، عندما يكون $s$ عدداً زوجياً سالباً ($-2, -4, -6, \dots$)، تصبح $\sin(\pi s / 2)$ صفراً، وبالتالي تكون $\zeta(s) = 0$. وتُعرف هذه بـ الأصفار البديهية (Trivial zeros).

ومع ذلك، ما يهم في توزيع الأعداد الأولية هي الأصفار الأخرى، وتحديداً الأصفار غير البديهية (Non-trivial zeros) التي تقع في “الشريط الحرج (Critical strip)” حيث $0 \le \sigma \le 1$.

5. جوهر فرضية ريمان والصيغة الصريحة

حسب ريمان عدداً قليلاً من الأصفار وطرح فرضية مذهلة. هذه هي فرضية ريمان.

فرضية ريمان (Riemann Hypothesis) جميع الأصفار غير البديهية لدالة زيتا لريمان $\zeta(s)$ تقع على خط مستقيم جزؤه الحقيقي يساوي $1/2$ (أي $\text{Re}(s) = 1/2$).

يسمى هذا الخط الذي يكون فيه الجزء الحقيقي 1/2 بـ “الخط الحرج (Critical line)”.

graph TD A["دالة زيتا لريمان ζ(s)"] --> B["التوسيع إلى المستوى المركب بواسطة الامتداد التحليلي"] B --> C["الأصفار البديهية (s = -2, -4, -6 ...)"] B --> D["الأصفار غير البديهية (0 <= Re(s) <= 1)"] D --> E["فرضية ريمان"] E --> F["جميع الأصفار غير البديهية تقع على Re(s) = 1/2"] F --> G["تؤدي إلى إثبات نهاية مصطلح الخطأ في توزيع الأعداد الأولية"]

لماذا تعتبر فرضية ريمان على هذا القدر من الأهمية؟ لأن أصفار دالة زيتا تحدد بالكامل توزيع الأعداد الأولية.

استنتج ريمان ولاحقاً عالم الرياضيات فون مانجولدت “الصيغة الصريحة (Explicit formula)” التي تصف توزيع الأعداد الأولية بدقة. باستخدام دالة تشيبيشيف $\psi(x)$، يمكن التعبير عنها كما يلي:

$$ \psi(x) = x - \sum_{\rho} \frac{x^\rho}{\rho} - \ln(2\pi) - \frac{1}{2}\ln(1 - x^{-2}) $$

حيث $\rho$ يمثل المجموع عبر جميع الأصفار غير البديهية لدالة زيتا. المصطلح الرئيسي هو $x$ (والذي يتوافق مع مبرهنة الأعداد الأولية)، ومن خلال إضافة أو طرح المصطلحات الشبيهة بالتموجات التي تعتمد على الأصفار $\rho$، يتم إعادة بناء التوزيع الدقيق المتدرج للأعداد الأولية. يمكن القول إن الأصفار غير البديهية تمثل “ترددات (موجات)” توزيع الأعداد الأولية.

إذا كانت فرضية ريمان صحيحة، وكان الجزء الحقيقي لجميع الأصفار غير البديهية $\rho$ هو بالضبط $1/2$، فإن مصطلح الخطأ في مبرهنة الأعداد الأولية سيندرج ضمن أصغر نطاق يمكن تصوره نظرياً.

$$ |\pi(x) - \text{Li}(x)| \le \frac{1}{8\pi} \sqrt{x} \ln x \quad \text{for} \quad x \ge 2657 $$

بعبارة أخرى، إذا كانت فرضية ريمان صحيحة، فسيثبت ذلك أن الأعداد الأولية موزعة بأكثر الطرق “انتظاماً وجمالاً” التي يمكننا تخيلها.

6. العلاقة التي لا تنفصل بين تقنيات التشفير الحديثة والأعداد الأولية

حتى هذه النقطة كنا في عالم الرياضيات البحتة العميق، ولكن هذه الخصائص للأعداد الأولية تدعم من الأساس مجتمعنا الرقمي الحديث. والمثال الأبرز على ذلك هو تشفير المفتاح العام مثل تشفير RSA.

تعتمد أمان كافة الاتصالات مثل مدفوعات بطاقات الائتمان على الإنترنت، وإرسال كلمات المرور، والتوقيعات الإلكترونية لتقنية البلوك تشين على “الأعداد الأولية”.

كيفية عمل تشفير RSA

يعتمد أمان تشفير RSA على الحقيقة الرياضية القائلة بأن “تحليل رقم مركب كبير جداً إلى عوامله الأولية يمثل صعوبة بالغة” (مسألة التحليل إلى العوامل الأولية).

توليد المفاتيح: يتم اختيار عددين أوليين ضخمين بشكل عشوائي $p$ و $q$ (مثلاً 2048 بت لكل منهما). يتم ضربهما معاً لحساب $N = p \times q$. هذا الـ $N$ يصبح جزءاً من المفتاح العام. باستخدام دالة مؤشر أويلر $\phi(N) = (p-1)(q-1)$، يتم توليد المفتاح الخاص $d$.
$$ e \times d \equiv 1 \pmod{\phi(N)} $$
التشفير وفك التشفير: يتم تحويل النص الواضح $M$ إلى نص مشفر $C$ باستخدام المفتاح العام $e, N$.
$$ C \equiv M^e \pmod{N} $$
فقط من يمتلك المفتاح الخاص $d$ يمكنه فك التشفير.
$$ M \equiv C^d \pmod{N} $$

graph LR A["النص الواضح (Plaintext)"] --> B["التشفير بالمفتاح العام (e, N)"] B --> C["النص المشفر (Ciphertext)"] C --> D["فك التشفير بالمفتاح الخاص (d)"] D --> E["النص الواضح الأصلي"] F["المهاجم (Attacker)"] -- "محاولة تحليل N إلى عوامله الأولية" --> C F -.-> G["بدون معرفة p و q، لا يمكن حساب d"]

لكسر تشفير RSA، من الضروري إيجاد الأعداد الأولية الأصلية $p$ و $q$ (التحليل إلى عوامل) انطلاقاً من العدد الضخم $N$. حتى باستخدام الخوارزميات السائدة حالياً (مثل منخل حقل الأعداد العام: GNFS)، فإن تحليل عدد يتكون من مئات الأرقام قد يستغرق وقتاً يتجاوز بكثير عمر الكون، حتى باستخدام أجهزة الكمبيوتر الفائقة.

7. تأثير فرضية ريمان على تقنيات التشفير

إذن، كيف تتقاطع “فرضية ريمان” الموجودة في قمة الرياضيات البحتة مع “تقنيات التشفير”؟

7.1. خوارزميات توليد الأعداد الأولية (اختبار الأولية) وفرضية ريمان المعممة (GRH)

لتشغيل تشفير RSA، يجب أولاً توليد أعداد أولية ضخمة $p$ و $q$. ومع ذلك، فإن التحديد السريع والمؤكد لـ “ما إذا كان عدد ما أولياً أم لا” ليس أمراً سهلاً.

حالياً، تستخدم الخوارزمية الاحتمالية المعروفة بـ اختبار ميلر-رابين لأولية عدد (Miller-Rabin primality test) بشكل عملي. هذه الخوارزمية سريعة جداً، ولكنها تحمل خطراً ضئيلاً جداً لحدوث “أعداد شبه أولية”، حيث يتم تحديد عدد مركب بشكل خاطئ على أنه عدد أولي.

ومع ذلك، إذا افترضنا أن “فرضية ريمان المعممة (Generalized Riemann Hypothesis, GRH)”، التي توسع فرضية ريمان إلى دوال L لدركليه، صحيحة، فإن القصة تتغير بشكل جذري. إذا كانت GRH صحيحة، فإن الحد الأقصى لعدد الاختبارات في خوارزمية ميلر-رابين يصبح مضموناً رياضياً، وترتقي من خوارزمية احتمالية إلى “خوارزمية حتمية ذات وقت كثير الحدود” (كانت هذه حقيقة مهمة معروفة قبل اكتشاف اختبار AKS للأولية).

بمعنى آخر، تلعب فرضية ريمان (وتعميماتها) دوراً في إعطاء ضمانات مباشرة وقاطعة لتوليد الأعداد الأولية الضخمة بسرعة وموثوقية مطلقة، وهو ما يمثل الأساس لتوليد التشفير.

7.2. العلاقة بخوارزميات التحليل إلى العوامل

تعتبر معرفة توزيع الأعداد الأولية أمراً حتمياً عند تقييم التعقيد الحسابي للخوارزميات المستخدمة لفك التشفير (مثل منخل حقل الأعداد العام). تعتمد العديد من خوارزميات التحليل إلى عوامل على توزيع “الأعداد السلسة (Smooth numbers: الأعداد التي تحتوي فقط على عوامل أولية صغيرة)”.

لتقييم مدى تكرار ظهور الأعداد السلسة بدقة، يلزم وجود فهم عميق لتوزيع الأعداد الأولية، وهنا أيضاً يتم استخدام تقنيات نظرية الأعداد التحليلية المرتبطة مباشرة بدالة زيتا وفرضية ريمان. إذا تم إثبات فرضية ريمان وتم تحديد خطأ توزيع الأعداد الأولية بالكامل، فسيصبح من الممكن تحديد حدود الأداء لخوارزميات التحليل إلى عوامل بشكل أكثر دقة.

8. إذا تم إثبات فرضية ريمان، هل سيتم كسر التشفير؟

غالباً ما يُقال كنوع من الأساطير الحضرية: “إذا حُلت فرضية ريمان، سينهار تشفير RSA في لحظة”، ولكن هذا غير دقيق رياضياً.

إثبات فرضية ريمان بحد ذاته لا يخلق فوراً خوارزمية سحرية تسرع بشكل كبير من التحليل إلى عوامل أولية. فرضية ريمان هي مجرد نظرية حول “الانتظام الماكروسكوبي لتوزيع” الأعداد الأولية، ولا تخبرنا مباشرة (كخاصية محلية) عن الأعداد الأولية التي تقسم عدداً معيناً $N$.

ومع ذلك، فإن التأثير ليس صفراً. لأنه في عملية إثبات فرضية ريمان، هناك احتمال كبير جداً لاكتشاف “أدوات رياضية جديدة” أو “أساليب تحليل غير معروفة”. بالنظر إلى التاريخ، عندما تم إثبات مبرهنة فيرما الأخيرة أو حدسية بوانكاريه، فإن النظريات الجديدة التي طُورت خلال تلك العملية دفعت الرياضيات بأكملها إلى قفزات هائلة.

إذا تم تأسيس طرق هندسية جبرية غير معروفة أو أساليب هندسية غير تبادلية قادرة على التحكم الكامل في خصائص أصفار دالة زيتا لريمان، فلا يمكن إنكار احتمال أن يؤدي ذلك كنتيجة إلى اكتشاف خوارزميات ثورية للتحليل إلى العوامل (على سبيل المثال، خوارزمية كلاسيكية تقلل وقت الحساب إلى وقت كثير الحدود). بهذا المعنى، لا يمكن لعلماء التشفير أن يغضوا الطرف أبداً عن تطورات فرضية ريمان.

الحواسيب الكمومية وخوارزمية شور

التهديد الأكثر مباشرة وواقعية لتقنيات التشفير ليس إثبات فرضية ريمان، بل الحواسيب الكمومية. “خوارزمية شور” التي نشرها بيتر شور (Peter Shor) في عام 1994، أثبتت أنه بوجود حاسوب كمي ذو أداء كافٍ، يمكن حل التحليل إلى العوامل في وقت كثير الحدود. هذا يعني أن تشفير RSA وتشفير المنحنيات الإهليلجية سيتم كسرهما بشكل جذري.

حالياً، يجري التحول عالمياً نحو “تشفير ما بعد الكوانتم (Post-Quantum Cryptography, PQC)” الذي لا يمكن فك تشفيره حتى بواسطة الحواسيب الكمومية (مثل التشفير القائم على الشبكات). ربما يقترب العصر الذهبي لتقنيات التشفير المعتمدة على الأعداد الأولية من نهايته بطريقة ما، لكن القيمة الرياضية للأعداد الأولية بحد ذاتها لن تضيع أبداً.

9. الخاتمة: تقاطع التجريد الرياضي مع المجتمع الواقعي

graph TD A["استكشاف الرياضيات البحتة"] --> B["توضيح فرضية ريمان"] B --> C["فهم كامل لتوزيع الأعداد الأولية"] C --> D["تطور هائل في نظرية الأعداد والهندسة الجبرية"] D -.-> E["إمكانية خوارزميات جديدة للتحليل إلى العوامل"] E -.-> F["تحديث تقييم أمان تقنيات التشفير"] A --> G["الرياضيات التطبيقية وعلوم الكمبيوتر"] G --> H["تحسين كفاءة اختبار الأولية وتوليد التشفير"] H --> F

الاستكشاف اللانهائي للأعداد الأولية، والذي استمر منذ اليونان القديمة، ارتقى بفضل العبقري ريمان ليصبح سيمفونية جميلة على المستوى المركب (أصفار دالة زيتا). ومن المثير للدهشة أن هذه البلورة الرياضية النقية، بعد مرور قرون، تُستخدم كأقوى درع لضمان أمن مجتمع الإنترنت.

فرضية ريمان هي رمز يجسد في وقت واحد “الجمال التجريدي” للرياضيات و “القدرة المذهلة على التطبيق في العالم المادي والمجتمع الواقعي”.

عندما يتم قهر هذا الجبل الرياضي الضخم الذي لم يصل إلى قمته أحد بعد، سنفهم تماماً الحقيقة الكونية المتمثلة في الأعداد الأولية، وسنمتلك منظوراً جديداً لبنية مجتمع المعلومات. إن دراسة تقنيات التشفير هي، في حد ذاتها، رحلة لتتبع تاريخ حكمة البشرية.

خوارزمية سريعة لاختبار الأولية مطبقة بلغة C++ (مثل طريقة Miller-Rabin)

Fri, 11 Sep 2026 14:00:00 +0900

مقدمة: لماذا نحتاج إلى خوارزمية سريعة لاختبار الأولية؟

في عالم علوم الحاسوب، ونظرية التشفير، أو البرمجة التنافسية، يعتبر التحديد السريع والدقيق لـ “ما إذا كان عدد ما أولياً أم لا” تحدياً أساسياً وهاماً جداً. على سبيل المثال، يعتمد أمان أنظمة تشفير المفتاح العام مثل تشفير RSA - والذي يمثل الأساس الأمني لمجتمع الإنترنت الحديث - على توليد أعداد أولية ضخمة وصعوبة ضربها (صعوبة تحليل العوامل الأولية). لذلك، ليس من المبالغة القول إن تقنية التمييز الفوري لما إذا كان عدد ضخم أولياً أم لا هي تقنية تدعم أساس المجتمع الرقمي.

بالإضافة إلى ذلك، في البرمجة التنافسية (مثل AtCoder و Codeforces)، يعتبر اختبار الأولية موضوعاً متكرراً. بالنسبة للمدخلات الضخمة التي قد تصل إلى $N \le 10^{18}$، في الحالات التي يتم فيها إجراء اختبارات الأولية عشرات الآلاف من المرات في غضون ثانية واحدة، لا يمكن للخوارزميات التقليدية البسيطة أن تلبي متطلبات وقت الحساب وتتسبب في (Time Limit Exceeded: TLE).

في هذا المقال، سنبدأ بخوارزمية اختبار الأولية البسيطة، ثم “اختبار فيرما” وهو اختبار احتمالي، وصولاً إلى الخوارزمية الأسرع والأكثر عملية “طريقة ميلر-رابين لاختبار الأولية (Miller-Rabin Primality Test)” التي تغلبت على نقاط الضعف في الخوارزميات السابقة. سنشرح ذلك بالتفصيل من الخلفية الرياضية إلى التنفيذ عالي التحسين بلغة C++. بشكل خاص، بالنسبة للأعداد الصحيحة بحجم 64 بت ($N < 2^{64}$)، سنشرح بالتفصيل طريقة “تحديد الأولية بنسبة 100% (الاختبار الحتمي)” بدلاً من مجرد التحديد الاحتمالي، وسنوفر كود مصدري بلغة C++ جاهز للاستخدام العملي.

1. أساسيات اختبار الأولية وطريقة القسمة التجريبية (Trial Division)

العدد الأولي (Prime number) هو عدد طبيعي أكبر من أو يساوي 2، وليس له أي قواسم موجبة سوى 1 والعدد نفسه. باتباع تعريف العدد الأولي بدقة، لتحديد ما إذا كان العدد الصحيح $N$ أولياً، يمكننا محاولة قسمة $N$ على جميع الأعداد الصحيحة من $2$ إلى $N-1$. إذا لم يقبل القسمة على أي منها فهو عدد أولي، وإذا قَبِل القسمة ولو لمرة واحدة فهو عدد مؤلف (ليس أولياً).

ولكن، تعقيد الوقت لهذه الطريقة هو $O(N)$، وإذا كان $N$ عدداً ضخماً مثل $10^{18}$، فسيستغرق الأمر وقتاً هائلاً لحسابه حتى مع أجهزة الكمبيوتر الحديثة.

تحسين طريقة القسمة التجريبية: البحث حتى $\sqrt{N}$

عندما يتم التعبير عن العدد المؤلف $N$ كـ $a \times b = N$ (حيث $a \le b$)، فمن المؤكد أن $a \le \sqrt{N}$. لذلك، لا توجد حاجة لتكرار حلقة اختبار الأولية حتى $N-1$، ويكفي التحقق حتى $\sqrt{N}$.

 1
 2
 3
 4
 5
 6
 7
 8
 9
10
11
12
13
14


#include <iostream>

// اختبار الأولية باستخدام طريقة القسمة التجريبية (O(sqrt(N)))
bool is_prime_trial_division(long long n) {
 if (n <= 1) return false;
 if (n == 2 || n == 3) return true;
 if (n % 2 == 0) return false;

 // فحص الأعداد الفردية فقط من 3 فصاعداً
 for (long long i = 3; i * i <= n; i += 2) {
 if (n % i == 0) return false;
 }
 return true;
}

التعقيد الزمني لهذه الخوارزمية هو $O(\sqrt{N})$. إذا كان $N \le 10^{12}$ تقريباً، فيمكن حسابه في لحظة، ولكن في حالة $N \approx 10^{18}$، يصبح عدد التكرارات حوالي $10^9$ مرة، وحتى مع مترجمات C++، سيستغرق ذلك من مئات الميلي ثانية إلى عدة ثوانٍ، مما يجعله غير مناسب للاختبارات المتعددة.

2. اختبار فيرما: بداية اختبارات الأولية الاحتمالية

لتجاوز حدود طريقة القسمة التجريبية، تم ابتكار “الخوارزميات الاحتمالية (Probabilistic Algorithm)” باستخدام نظريات من نظرية الأعداد. والمثال النموذجي لذلك هو “اختبار فيرما للأولية (Fermat Primality Test)” الذي يستخدم مبرهنة فيرما الصغرى.

مبرهنة فيرما الصغرى (Fermat’s Little Theorem)

اكتشف بيير دي فيرما (Pierre de Fermat) هذه النظرية، والتي تنص على ما يلي:

لأي عدد أولي $p$، ولأي عدد صحيح $a$ أولي نسبياً مع $p$ (أي ليس من مضاعفات $p$)، تتحقق معادلة التطابق التالية:
$$ a^{p-1} \equiv 1 \pmod p $$

إذا أخذنا المعاكس الإيجابي لهذه النظرية، فيمكننا القول: “بالنسبة لعدد صحيح $N$، وعدد صحيح $a$ أولي نسبياً مع $N$، إذا كان $a^{N-1} \not\equiv 1 \pmod N$، فمن المؤكد أن $N$ عدد مؤلف”. باستخدام هذه الخاصية، يتم اختيار أساس (base) عشوائي $a$ للعدد المراد اختباره $N$، وحساب $a^{N-1} \pmod N$ للتحقق مما إذا كان يساوي $1$. هذا هو اختبار فيرما.

الأسس النمطية السريعة (التربيع المتكرر)

لإجراء اختبار فيرما، يجب حساب الأس الضخم $a^{N-1} \pmod N$ بسرعة. لهذا، نستخدم “التربيع المتكرر (Modular Exponentiation / Binary Exponentiation)”. التعقيد الزمني هو $O(\log N)$، وهو سريع جداً.

 1
 2
 3
 4
 5
 6
 7
 8
 9
10
11


// حساب a^b mod m باستخدام التربيع المتكرر
long long mod_pow(long long a, long long b, long long m) {
 long long res = 1;
 a %= m;
 while (b > 0) {
 if (b & 1) res = (__int128_t)res * a % m;
 a = (__int128_t)a * a % m;
 b >>= 1;
 }
 return res;
}

※ هنا، لمنع حدوث فيضان الأعداد (Overflow)، نستخدم امتداد GCC/Clang المسمى __int128_t (عدد صحيح بحجم 128 بت) للاحتفاظ بالضرب الوسيط.

الأعداد شبه الأولية وأعداد كارمايكل (Carmichael Numbers)

على الرغم من أن اختبار فيرما قوي للغاية، إلا أن لديه نقطة ضعف قاتلة. وهي وجود أعداد $N$ مؤلفة، ولكنها تحقق $a^{N-1} \equiv 1 \pmod N$ لجميع قيم $a$ (الأولية نسبياً مع $N$).

تسمى هذه الأعداد بـ “الأعداد شبه الأولية المطلقة” أو “أعداد كارمايكل”. أصغر عدد كارمايكل هو $561 = 3 \times 11 \times 17$. نظراً لوجود أعداد كارمايكل، لا يمكن إجراء اختبار حتمي بـ “احتمال 100%” باستخدام اختبار فيرما وحده. بغض النظر عن قيم $a$ المختلفة التي نختبرها، فإن عدداً مثل $561$ سيتظاهر دائماً بأنه عدد أولي (يخدعنا).

3. طريقة ميلر-رابين لاختبار الأولية (Miller-Rabin Primality Test)

لقد تغلب كل من غاري ميلر (Gary L. Miller) ومايكل رابين (Michael O. Rabin) ببراعة على نقطة ضعف اختبار فيرما (وجود أعداد كارمايكل) من خلال ابتكار “طريقة ميلر-رابين لاختبار الأولية”. حالياً، كخوارزمية عملية وسريعة لاختبار الأولية، تُستخدم على نطاق واسع في المكتبات الداخلية لمختلف لغات البرمجة وتوليد المفاتيح في أنظمة التشفير.

المبادئ الرياضية

بالإضافة إلى مبرهنة فيرما الصغرى، تستخدم خوارزمية ميلر-رابين الخاصية القائلة بأنه “في حقل البواقي للعدد الأولي ($\mathbb{Z}/p\mathbb{Z}$)، تقتصر حلول $x^2 \equiv 1 \pmod p$ على $x \equiv 1$ أو $x \equiv -1$” (في حالة كان المعامل عدداً مؤلفاً، قد توجد جذور تربيعية غير بديهية أخرى).

العدد الفردي $N$ المراد اختباره مطروحاً منه $1$، أي $N-1$، سيكون بالتأكيد عدداً زوجياً. لذلك، نقوم بقسمة $N-1$ على $2$ قدر الإمكان، ونعبر عنه بالشكل التالي:

$$ N-1 = d \cdot 2^s $$

(حيث $d$ عدد فردي، و $s \ge 1$)

بالنسبة لأي أساس $a$ ($1 < a < N-1$)، نتحقق مما إذا كان $a^{N-1} \equiv 1 \pmod N$ وفقاً لمبرهنة فيرما الصغرى، ولكننا نقوم بهذا الحساب تدريجياً. على وجه التحديد، نقوم بتكرار التربيع بالترتيب التالي: $a^d, a^{d \cdot 2}, a^{d \cdot 4}, \ldots, a^{d \cdot 2^s}$.

لكي تقرر طريقة ميلر-رابين أن $N$ “عدد أولي (أو عدد أولي باحتمال كبير)"، يجب أن يتحقق أي من الشروط التالية:

$a^d \equiv 1 \pmod N$
يوجد $r$ ($0 \le r < s$) بحيث يتحقق $a^{d \cdot 2^r} \equiv -1 \pmod N$. ※ في عملية المعامل (Modulo) بلغة C++، يصبح $-1 \pmod N$ مساوياً لـ $N-1$.

إذا كان $N$ عدداً أولياً، فإن هذا الشرط سيتحقق بالتأكيد لأي $a$. على العكس، إذا كان $N$ عدداً مؤلفاً، فقد تم إثبات رياضياً أن احتمال تحقق هذا الشرط (احتمال الانخداع) عند اختيار $a$ عشوائياً هو أقل من أو يساوي $\frac{1}{4}$. إذا قمنا بإجراء $k$ اختبارات مستقلة، فإن احتمال الخطأ في التحديد يصبح أقل من أو يساوي $\left(\frac{1}{4}\right)^k$، ويمكن اعتباره عملياً صفراً. لا توجد أعداد “يمكنها الخداع تماماً” مثل أعداد كارمايكل.

سير خوارزمية ميلر-رابين (مخطط Mermaid)

يوضح الشكل التالي التدفق المنطقي لاختبار واحد في طريقة ميلر-رابين (اختبار لأساس واحد $a$).

graph TD Start["بدء الاختبار (الإدخال: N, a)"] --> CalcDS["حساب d (فردي) و s التي تحقق N-1 = d * 2^s"] CalcDS --> CalcX["حساب x = a^d mod N"] CalcX --> CheckX1{"هل x == 1 أو x == N-1 ؟"} CheckX1 -- "نعم" --> ReturnTrue["يحتمل أن يكون أولياً (Probably Prime)"] CheckX1 -- "لا" --> LoopStart["بدء حلقة r من 1 إلى s-1"] LoopStart --> LoopCondition{"هل r < s ؟"} LoopCondition -- "لا" --> ReturnFalse["مؤلف بالتأكيد (Composite)"] LoopCondition -- "نعم" --> SquareX["حساب x = (x * x) mod N"] SquareX --> CheckXMinus1{"هل x == N - 1 ؟"} CheckXMinus1 -- "نعم" --> ReturnTrue CheckXMinus1 -- "لا" --> CheckXOne{"هل x == 1 ؟"} CheckXOne -- "نعم" --> ReturnFalse CheckXOne -- "لا" --> LoopNext["زيادة r بمقدار 1 والانتقال للتالي"] LoopNext --> LoopCondition

4. الاختبار الحتمي للأعداد الصحيحة بحجم 64 بت

على الرغم من أن طريقة ميلر-رابين هي أساساً خوارزمية “احتمالية”، إلا أنه عندما يكون الحد الأقصى لـ $N$ ثابتاً، يمكننا إجراء اختبار أولية “مؤكد بنسبة 100%” من خلال تجربة مجموعة محددة من قيم $a$ (الأساس). يسمى هذا اختبار ميلر-رابين الحتمي (Deterministic Miller-Rabin Test).

من خلال بحث Jim Sinclair وآخرين، تبين أنه بالنسبة لجميع الأعداد الصحيحة $N < 2^{64}$ (حوالي $1.8 \times 10^{19}$)، يكفي اختيار الأعداد الأولية السبعة التالية كأساس $a$ واختبارها للحصول على تقييم حتمي وكامل بما فيه الكفاية.

قائمة الأساس $a$ المراد اختبارها: {2, 325, 9375, 28178, 450775, 9780504, 1795265022}

أو، كمجموعة أخرى معروفة، يمكن التحديد بشكل مثالي لـ $N < 2^{64}$ باستخدام الأعداد الأولية الـ $12$ التالية: {2, 3, 5, 7, 11, 13, 17, 19, 23, 29, 31, 37}

هذه المرة، من أجل تعزيز بساطة الخوارزمية وموثوقيتها، سنعتمد طريقة تستخدم الأعداد الأولية الـ $12$ الأخيرة كقاعدة (أو القواعد السبعة المحسنة). في تنفيذ C++، نقوم بتحسين يقلل من عدد الاختبارات إلى الحد الأدنى عن طريق تقسيم النطاقات باستخدام التفرع الشرطي.

5. تنفيذ متقدم بلغة C++ (Highly Optimized C++ Implementation)

الآن، سنقوم بتجميع النظرية الرياضية وتصميم الخوارزمية التي تناولناها حتى الآن، وسنقدم الكود المصدري لتنفيذ دالة اختبار ميلر-رابين للأولية بأعلى مستوى في لغة C++ الحديثة.

نقاط التنفيذ

تجنب فيضان الضرب للأعداد الصحيحة 64 بت: في حالة $N \approx 10^{18}$، يمكن أن يصل الحد الأقصى لـ $x \times x$ في الضرب النمطي إلى $10^{36}$، مما يتجاوز بسهولة الحد الأقصى للأنواع الصحيحة ذات 64 بت (uint64_t أو long long) وهو $1.8 \times 10^{19}$. لحل هذه المشكلة، نستخدم النوع الممتد لـ GCC و Clang وهو __int128_t (أو unsigned __int128)، ونحسب المعامل (Modulo) بعد الحساب بدقة 128 بت. هذا يسمح بالضرب النمطي السريع دون استخدام خوارزميات معقدة.
اختيار الأساس الحتمي: عندما تكون قيمة $N$ صغيرة، نقوم بالتحسين بحيث لا نحتاج إلا لاختبار عدد قليل من الأسس (bases).

كود C++ المصدري الكامل

فيما يلي نعرض الكود المصدري للنسخة النهائية الجاهزة للاستخدام العملي. يمكن نسخ هذا الكود واستخدامه كما هو في بيئات البرمجة التنافسية وغيرها.

 1
 2
 3
 4
 5
 6
 7
 8
 9
10
11
12
13
14
15
16
17
18
19
20
21
22
23
24
25
26
27
28
29
30
31
32
33
34
35
36
37
38
39
40
41
42
43
44
45
46
47
48
49
50
51
52
53
54
55
56
57
58
59
60
61
62
63
64
65
66
67
68
69
70
71
72
73
74
75
76
77
78
79
80
81
82
83
84
85
86
87
88
89
90
91
92
93
94
95
96


#include <iostream>
#include <vector>
#include <cstdint>
#include <initializer_list>

using namespace std;

// (a * b) mod m السريع باستخدام أعداد صحيحة 128 بت
inline uint64_t mod_mul(uint64_t a, uint64_t b, uint64_t m) {
 return (uint64_t)((unsigned __int128)a * b % m);
}

// حساب (base^exp) mod m باستخدام التربيع المتكرر
uint64_t mod_pow(uint64_t base, uint64_t exp, uint64_t m) {
 uint64_t res = 1;
 base %= m;
 while (exp > 0) {
 if (exp & 1) res = mod_mul(res, base, m);
 base = mod_mul(base, base, m);
 exp >>= 1;
 }
 return res;
}

// الاختبار الحتمي للأعداد الصحيحة 64-بت باستخدام طريقة ميلر-رابين
bool is_prime_miller_rabin(uint64_t n) {
 // التحقق المسبق من القيم الحدية والأعداد الأولية الصغيرة
 if (n < 2) return false;
 if (n == 2 || n == 3 || n == 5 || n == 7) return true;
 if (n % 2 == 0 || n % 3 == 0 || n % 5 == 0 || n % 7 == 0) return false;

 // تحليل إلى الشكل n-1 = d * 2^s
 uint64_t d = n - 1;
 int s = 0;
 while ((d & 1) == 0) {
 d >>= 1;
 s++;
 }

 // قائمة الأسس (bases) المستخدمة في الاختبار
 // تحسين يقلل من عدد الأسس المختبرة بناءً على حجم N
 vector<uint64_t> bases;
 if (n < 4759123141ULL) {
 bases = {2, 7, 61};
 } else if (n < 1122004669633ULL) {
 bases = {2, 13, 23, 1662803};
 } else {
 // 7 أسس حتمية لجميع الأعداد N < 2^64
 bases = {2, 325, 9375, 28178, 450775, 9780504, 1795265022};
 }

 // تنفيذ الاختبار لكل أساس
 for (uint64_t a : bases) {
 a %= n;
 if (a == 0) continue; // إذا كان a من مضاعفات n فلا يمكن التحديد ولكنه ليس أولياً

 uint64_t x = mod_pow(a, d, n);
 if (x == 1 || x == n - 1) continue; // اجتياز الشرط الأول، الانتقال للأساس التالي

 bool composite = true;
 // حلقة التكرار s-1 مرة (x = x^2 mod n)
 for (int r = 1; r < s; r++) {
 x = mod_mul(x, x, n);
 if (x == n - 1) {
 composite = false; // اجتياز الشرط الثاني، يحتمل أن يكون أولياً
 break;
 }
 }

 // إذا لم يستوفِ أي شرط، فهو بالتأكيد عدد مؤلف
 if (composite) return false;
 }

 // إذا اجتاز الشروط لجميع الأسس، فهو أولي بالتأكيد
 return true;
}

int main() {
 // أمثلة للاختبار
 vector<uint64_t> test_cases = {
 1000000007, // عدد أولي مشهور
 998244353, // عدد أولي مشهور
 1000000000000000003, // عدد أولي بالقرب من 10^18
 1000000000000000007, // عدد مؤلف (10^18 + 7)
 561, // عدد كارمايكل (مؤلف)
 18446744073709551557ULL // أحد أكبر الأعداد الأولية بالقرب من 2^64
 };

 for (uint64_t n : test_cases) {
 cout << n << " is "
 << (is_prime_miller_rabin(n) ? "Prime" : "Composite")
 << endl;
 }

 return 0;
}

6. تقييم أداء وتعقيد الخوارزمية

سنقوم بدراسة أداء الخوارزمية التي تم تنفيذها.

التعقيد الزمني (Time Complexity)

طريقة القسمة التجريبية: $O(\sqrt{N})$
اختبار فيرما: حساب الأس $O(\log N) \times k$ ($k$ هو عدد المحاولات)
طريقة ميلر-رابين: حساب الأس والتكرار $O(\log N) \times k$

في بيئة 64 بت ($N \le 2^{64}$)، تتحقق طريقة ميلر-رابين الحتمية المذكورة أعلاه من $7$ أسس كحد أقصى. لذلك، يمكن اعتبار $k \le 7$ كثابت، ويكون التعقيد الزمني الإجمالي هو $O(\log N)$ بدقة. حتى في الحالة القصوى ($N \approx 10^{19}$)، لا يتجاوز عدد خطوات التنفيذ $7 \times 64 = 448$ خطوة من العمليات الأساسية، ويكون وقت التنفيذ بضعة ميكروثانية ($10^{-6}$ ثانية) أو أقل. مقارنة بـ $O(\sqrt{N})$ لطريقة القسمة التجريبية (عدد التكرارات $\approx 4 \times 10^9$ مرة)، يتم تحقيق تسريع بملايين المرات.

المزيد من التحسين: ضرب مونتغمري (Montgomery Multiplication)

في تنفيذ هذا المقال، نستخدم النوع الممتد 128 بت __int128_t لإجراء القسمة (عملية المعامل %). حتى في وحدات المعالجة المركزية الحديثة، تُعتبر قسمة الأعداد الصحيحة (تعليمات DIV) أمرًا مكلفًا يتطلب عشرات الدورات مقارنة بالجمع والضرب.

قد يستخدم صانعو المكتبات والمبرمجون التنافسيون الذين يسعون إلى التحسين الأقصى تقنية تُعرف باسم ضرب مونتغمري (Montgomery Multiplication). ضرب مونتغمري هو خوارزمية مذهلة تستبدل عملية المعامل (القسمة) المكلفة بـ “إزاحة البتات (Bit Shifts) والضرب فقط” عن طريق تعيين الأرقام إلى “مساحة مونتغمري” الخاصة. بدمج هذا في الضرب النمطي لاختبار ميلر-رابين، يمكن زيادة سرعة التنفيذ بمرتين إلى ثلاث مرات. نظراً لأن هذا موضوع عميق جداً، أود أن أشرحه بالتفصيل في مقال آخر.

7. الخلاصة

في هذا المقال، قمنا بشرح شامل من أساسيات اختبار الأولية إلى المواضيع المتقدمة. دعونا نراجع النقاط الرئيسية.

طريقة القسمة التجريبية موثوقة، ولكن نظرًا لتعقيدها الزمني $O(\sqrt{N})$، تفتقر إلى الفائدة العملية عندما يتجاوز $N$ $10^{12}$.
اختبار فيرما سريع جدًا بـ $O(\log N)$، لكن له نقطة ضعف قاتلة تتمثل في انخداعه بالأعداد شبه الأولية المطلقة مثل أعداد كارمايكل.
طريقة ميلر-رابين لاختبار الأولية هي الخوارزمية الأكثر عملية وقوة والتي تقضي على نقاط ضعف اختبار فيرما.
في تنفيذ C++، يمكن التعامل بأمان مع فيضان ضرب الأعداد الصحيحة 64 بت من خلال الاستفادة من __int128_t.
ضمن نطاق الأعداد الصحيحة 64 بت ($N < 2^{64}$)، من خلال اختيار $7$ أو $12$ عدداً أولياً محدداً كأساس، يمكن اختبار الأولية بشكل حتمي (دقيق 100%) وليس باحتمال.

يعد اختبار الأولية السريع تقنية لا مفر منها في الحسابات التي تتعامل مع أعداد ضخمة. كود المصدر لطريقة ميلر-رابين بـ C++ المقدم في هذا المقال قوي ويمكن استخدامه عملياً كما هو. يرجى الاستفادة منه في مشاريعك الخاصة أو مسابقات الخوارزميات.

graph LR TrialDivision["طريقة القسمة التجريبية (O(√N))"] --> Fermat["اختبار فيرما (O(log N), توجد نقاط ضعف)"] Fermat --> MillerRabin["طريقة ميلر-رابين (O(log N), يمكن جعلها حتمية)"] MillerRabin --> Montgomery["+ ضرب مونتغمري (تسريع المضاعف الثابت)"] style MillerRabin fill:#f9f,stroke:#333,stroke-width:2px

عالم الخوارزميات حيث تتقاطع البرمجة والرياضيات هو عالم جميل وعميق جداً. نأمل أن يكون هذا المقال مفيداً لتعلمك المستقبلي.

Reference:

Pomerance, C., Selfridge, J. L., & Wagstaff, S. S. (1980). The pseudoprimes to 25.10^9. Mathematics of Computation.
Sinclair, J. (2011). Deterministic Miller-Rabin primality testing.

【شرح رياضيات】شرح آلية عمل تشفير RSA بطريقة يفهمها حتى طلاب الثانوية

Fri, 11 Sep 2026 13:00:00 +0900

إحدى التقنيات التي تدعم أمن مجتمع الإنترنت من جذوره هي “تشفير RSA”. العديد من الاتصالات التي نستخدمها يوميًا دون تفكير، مثل مدفوعات بطاقات الائتمان في التسوق عبر الإنترنت، والتواصل مع الأصدقاء على شبكات التواصل الاجتماعي، وإرسال واستقبال المعلومات السرية للشركات، محمية بواسطة تشفير RSA أو التقنيات التي خلفته.

ولكن عندما تسمع كلمة “تشفير”، قد تتخيل آلات تشفير معقدة كتلك التي تظهر في أفلام الجواسيس، أو رياضيات متقدمة للغاية لا يفهمها سوى بعض العباقرة. صحيح أن نظرية التشفير الحديثة تعتمد على رياضيات متقدمة، إلا أن الآلية الأساسية لتشفير RSA يمكن فهمها تمامًا إذا كان لديك معرفة برياضيات المرحلة الثانوية (خصائص الأعداد الصحيحة، الأعداد الأولية، التطابق وغيرها).

في هذا المقال، ومن منطلق المعرفة برياضيات المرحلة الثانوية، سنشرح خطوة بخطوة وبشكل شامل المبادئ الرياضية التي يعمل بها تشفير RSA، ولماذا يصعب فك تشفيره. سنشرح بعناية مع استخدام أمثلة ملموسة حتى يتمكن حتى من يجدون صعوبة قليلة في الرياضيات من الفهم.

1. التشفير بالمفتاح المتماثل والتشفير بالمفتاح المعلن

قبل الخوض في الآلية الرياضية لتشفير RSA، دعونا أولاً ننظم الفكرة الأساسية للتشفير. تنقسم طرق التشفير بشكل عام إلى نوعين: “التشفير بالمفتاح المتماثل (المشترك)” و"التشفير بالمفتاح المعلن (العام)".

1.1 حدود طريقة التشفير بالمفتاح المتماثل

العديد من الشفرات المستخدمة منذ القدم تسمى “طريقة التشفير بالمفتاح المتماثل”. وهي طريقة تستخدم نفس المفتاح في “التشفير (تحويل الرسالة إلى نص مشفر سري)” و"فك التشفير (إعادة النص المشفر إلى الرسالة الأصلية)".

على سبيل المثال، لنفترض أن أليس تريد إرسال رسالة سرية إلى بوب. تضع أليس الرسالة في صندوق وتغلقه باستخدام قفل (مفتاح متماثل). لكي يفتح بوب ذلك الصندوق، يحتاج إلى امتلاك نفس المفتاح الذي استخدمته أليس.

هذه الطريقة بها مشكلة كبيرة، وهي “مشكلة توزيع المفاتيح”. عندما تتواصل أليس وبوب، اللذان تفصل بينهما مسافة كبيرة، لأول مرة، كيف يمكنهما مشاركة المفتاح دون أن يتم التنصت عليهما؟ إذا سُرق المفتاح من قبل طرف ثالث أثناء إرساله بالبريد، فإن جميع الاتصالات المشفرة اللاحقة ستكون مكشوفة تمامًا.

1.2 الاختراع الثوري “طريقة التشفير بالمفتاح المعلن”

تم ابتكار “طريقة التشفير بالمفتاح المعلن” لحل مشكلة توزيع المفاتيح هذه. وتشفير RSA هو نوع من هذا التشفير.

في طريقة التشفير بالمفتاح المعلن، يتم استخدام مفتاحين مختلفين: “مفتاح للتشفير (المفتاح المعلن/العام)” و"مفتاح لفك التشفير (المفتاح السري/الخاص)".

يقوم المتلقي بوب بإنشاء زوج من “المفتاح المعلن” و"المفتاح السري".
ينشر بوب “المفتاح المعلن” للعالم أجمع (لا يهم من يحصل عليه).
تستخدم المرسلة أليس “المفتاح المعلن” الخاص ببوب لتشفير رسالتها وإرسالها.
لا يمكن فك تشفير الرسالة المشفرة إلا بواسطة “المفتاح السري” الذي يمتلكه بوب فقط.

إذا شبهنا ذلك بالقفل، فإن بوب يصنع الكثير من “الأقفال المفتوحة (المفاتيح المعلنة)” وينثرها في جميع أنحاء العالم. تضع أليس رسالتها الموجهة لبوب في صندوق وتغلقه باستخدام أحد أقفال بوب التي وجدتها. بمجرد إغلاق القفل، لا يمكن فتحه إلا باستخدام “المفتاح الرئيسي (المفتاح السري)” الذي يمتلكه بوب. حتى لو سرق شخص ما الصندوق في الطريق، فلا يمكنه فتحه لأنه لا يمتلك المفتاح الرئيسي.

graph TD A["أليس (المرسلة)"] --> B["نص صريح (رسالة)"] B --> C["عملية التشفير"] D["مفتاح بوب المعلن (متاح للجميع)"] --> C C --> E["إرسال عبر الإنترنت: نص مشفر"] E --> F["عملية فك التشفير"] G["مفتاح بوب السري (يمتلكه بوب فقط)"] --> F F --> H["النص الصريح المسترجع (الرسالة)"] H --> I["بوب (المتلقي)"]

لتحقيق هذا النظام الثوري، نحتاج إلى نوع من “الدوال أحادية الاتجاه (لغز رياضي في اتجاه واحد)”، حيث “يمكن التشفير بسهولة باستخدام المفتاح المعلن، ولكن من المستحيل مطلقًا فك التشفير بدون المفتاح السري”. القطعة التي تم التركيز عليها لحل هذا اللغز هي “الأعداد الأولية” التي نعرفها جيدًا.

2. الأساس الرياضي الذي يدعم تشفير RSA 1: الأعداد الأولية والتحليل إلى العوامل الأولية

يعتمد أمان تشفير RSA على حقيقة رياضية مفادها أن “تحليل الأعداد الضخمة إلى عواملها الأولية أمر صعب للغاية”.

2.1 ما هو العدد الأولي

العدد الأولي هو “عدد طبيعي أكبر من 1، لا يقبل القسمة إلا على 1 وعلى نفسه”. مثال: $2, 3, 5, 7, 11, 13, 17, 19, 23...$

الأعداد الأولية تشبه “الذرات” لجميع الأعداد الصحيحة. يمكن تحليل أي عدد طبيعي إلى شكل ضرب أعداد أولية. وهذا ما يسمى التحليل إلى العوامل الأولية. على سبيل المثال، إمكانية التحليل إلى عوامل أولية بطريقة واحدة فقط (إذا تجاهلنا الترتيب) كما في $60 = 2^2 \times 3 \times 5$ تُعرف باسم “المبرهنة الأساسية في الحساب”.

2.2 صعوبة التحليل إلى العوامل الأولية (الدالة أحادية الاتجاه)

المهم هنا هو عدم التماثل حيث أن “عملية الضرب سهلة، لكن التحليل إلى العوامل الأولية صعب”.

على سبيل المثال، حاول حساب حاصل ضرب العددين الأوليين التاليين ذهنيًا: $11 \times 13 = ?$ هذا سهل، أليس كذلك؟ الجواب هو $143$.

إذًا، ماذا عن العدد التالي؟ حلل العدد $323$ إلى عوامله الأولية. ما رأيك؟ سيستغرق الأمر بعض الوقت. (الجواب هو $17 \times 19$).

إذا كانت الأعداد صغيرة، فيمكن للإنسان حسابها بطريقة ما، ولكن عندما تصبح الأعداد كبيرة، يصبح الحساب صعبًا بشكل هائل حتى باستخدام أجهزة الكمبيوتر. في تشفير RSA السائد حاليًا، يتم استخدام عدد $N = p \times q$ ناتج عن ضرب عددين أوليين ضخمين بشكل لا يصدق $p$ و $q$ بحجم 2048 بت (حوالي 600 رقم بالنظام العشري).

عند إعطاء عددين أوليين ضخمين $p$ و $q$، فإن حساب $N$ يستغرق جزءًا من الثانية بالنسبة للكمبيوتر (أقل من مللي ثانية). ولكن على العكس من ذلك، إذا تم إعطاء $N$ فقط لإيجاد $p$ و $q$ الأصليين، فسيستغرق الأمر وقتًا طويلاً لا يمكن لأسرع حاسوب عملاق حالي حله حتى لو استمر في العمل لترليونات السنين.

هذا “عدم التماثل في الحساب (طريق الذهاب سهل، وطريق العودة صعب)” هو الأساس الذي يخلق العلاقة بين المفتاح المعلن والمفتاح السري.

3. الأساس الرياضي الذي يدعم تشفير RSA 2: التطابق (حساب المقياس)

لا تتم حسابات تشفير RSA بالجمع والضرب حيث تكبر الأعداد بشكل لا نهائي كما نستخدمها عادة، بل تتم في عالم “البواقي” بعد القسمة على عدد معين. وهذا ما يسمى التطابق (حساب المقياس أو الحساب المعياري).

3.1 رياضيات الساعة

غالبًا ما يُشبه حساب المقياس بـ “رياضيات الساعة”. إذا كانت الساعة الآن 10، فكم ستكون الساعة بعد 5 ساعات؟ إنها $10 + 5 = 15$، ولكن في نظام الساعة المعتاد المكون من 12 ساعة، نجيب “الساعة 3”. هذا لأن باقي قسمة 15 على 12 هو 3.

في عالم الرياضيات، يُكتب هذا على النحو التالي:

$$ 15 \equiv 3 \pmod{12} $$

ويُقرأ: “15 و 3 متطابقان بمقياس 12 (باقي قسمتهما على 12 متساوٍ)”.

3.2 الخصائص الأساسية للتطابق

للتطابق خصائص مفيدة جدًا تشبه إلى حد كبير خصائص المعادلة ($=$). لنفترض أن المقياس (القاسم) هو $N$. عندما يكون $a \equiv b \pmod N$ و $c \equiv d \pmod N$، يتحقق ما يلي:

الجمع: $a + c \equiv b + d \pmod N$
الطرح: $a - c \equiv b - d \pmod N$
الضرب: $a \times c \equiv b \times d \pmod N$
الرفع لأس: $a^k \equiv b^k \pmod N$ (حيث $k$ عدد طبيعي)

الخاصية الأكثر أهمية هي “الرفع لأس”. هذا يعني أن “أس الباقي يساوي باقي الأس”. على سبيل المثال، لنفترض أننا نريد إيجاد باقي قسمة $7^{100}$ على $5$. سيكون من الصعب جدًا ضرب $7$ في نفسها 100 مرة ثم القسمة على $5$، ولكن باستخدام خصائص التطابق، نظرًا لأن $7 \equiv 2 \pmod 5$، فإن $7^{100} \equiv 2^{100} \pmod 5$، مما يبسط الحساب بشكل هائل. هذه الخاصية لا غنى عنها في عالم التشفير لأننا نتعامل مع أسس لأعداد كبيرة جدًا.

4. الأساس الرياضي الذي يدعم تشفير RSA 3: مؤشر أويلر ومبرهنة أويلر

من هنا تبدأ الرياضيات السحرية التي تشكل جوهر تشفير RSA. ستظهر “مبرهنة أويلر”، وهي تعميم لـ “مبرهنة فيرما الصغرى”.

4.1 دالة مؤشر أويلر (دالة توتنت) $\phi(N)$

دالة مؤشر أويلر (دالة $\phi$) هي دالة تعيد، بالنسبة لعدد طبيعي معين $N$، “عدد الأعداد الطبيعية من 1 إلى $N$ التي تكون أولية نسبيًا مع $N$ (أي القاسم المشترك الأكبر بينهما هو 1)”.

دعونا نلقي نظرة على بعض الأمثلة:

$\phi(5)$: من بين الأعداد 1, 2, 3, 4, 5، الأعداد الأولية نسبيًا مع 5 هي 1, 2, 3, 4 (عددها 4). إذن $\phi(5) = 4$.
$\phi(6)$: من بين الأعداد 1, 2, 3, 4, 5, 6، الأعداد الأولية نسبيًا مع 6 هي 1, 5 (عددها 2). إذن $\phi(6) = 2$.

【خاصية مميزة في حالة العدد الأولي】 إذا كان $p$ عددًا أوليًا، فإن جميع الأعداد من 1 إلى $p-1$ تكون أولية نسبيًا مع $p$. وبالتالي:

$$ \phi(p) = p - 1 $$

【خاصية مميزة في حالة حاصل ضرب عددين أوليين】 بالنسبة لعددين أوليين مختلفين $p$ و $q$، إذا كان $N = p \times q$، فيمكن حساب $\phi(N)$ بسهولة كالتالي:

$$ \phi(N) = \phi(p) \times \phi(q) = (p - 1)(q - 1) $$

هذه الخاصية تعمل كـ “باب خلفي سري (Trapdoor)” لتشفير RSA. فالشخص الذي يعرف $p$ و $q$ (منشئ المفتاح) يمكنه حساب $\phi(N)$ في لحظة، لكن الطرف الثالث الذي يعرف $N$ فقط لا يمكنه حساب $\phi(N)$ ما لم يقم بتحليل $N$ إلى عوامله الأولية.

4.2 مبرهنة أويلر

أثبت ليونهارد أويلر هذه المبرهنة الجميلة باستخدام $\phi(N)$:

مبرهنة أويلر: عندما يكون العددان الصحيحان $a$ و $N$ أوليين نسبيًا، يتحقق التطابق التالي:

$$ a^{\phi(N)} \equiv 1 \pmod N $$

هذه خاصية مذهلة تعني أن “حاصل ضرب أي عدد $a$ في نفسه $\phi(N)$ مرة ثم قسمته على $N$، سيعطي باقيًا يساوي $1$ دائمًا”. (إذا كان $N$ عددًا أوليًا $p$، فإنها تصبح $a^{p-1} \equiv 1 \pmod p$، وتُعرف باسم مبرهنة فيرما الصغرى).

دعونا نعدل مبرهنة أويلر هذه. سنضرب كلا الطرفين في $a$ مرة أخرى:

$$ a^{\phi(N) + 1} \equiv a \pmod N $$

علاوة على ذلك، بالنسبة لأي عدد صحيح $k$، فإن $a^{k \cdot \phi(N)}$ سيكون أيضًا $1^k = 1$، وبالتالي تتحقق المعادلة التالية:

$$ a^{k \cdot \phi(N) + 1} \equiv a \pmod N $$

هذه المعادلة بالذات هي المبدأ الأساسي الذي يجعل سحر تشفير RSA “العودة إلى الحالة الأصلية بعد التشفير وفك التشفير” ممكنًا.

5. خوارزمية تشفير RSA: خطوات إنشاء المفتاح، التشفير، وفك التشفير

بعد أن اكتملت المعرفة الأساسية، دعونا أخيرًا نلقي نظرة على الخطوات المحددة لتشفير RSA. ينقسم تشفير RSA بشكل رئيسي إلى ثلاث مراحل: “1. إنشاء المفتاح”، “2. التشفير”، و"3. فك التشفير".

flowchart TD A1["1. اختيار عددين أوليين p, q"] --> A2["حساب N = p × q"] A1 --> A3["حساب φ(N) = (p-1)(q-1)"] A3 --> A4["اختيار e بحيث يكون أوليًا نسبيًا مع φ(N)"] A3 --> A5["حساب d بحيث يتحقق التطابق e × d ≡ 1 (mod φ(N))"] A2 --> A6["المفتاح المعلن (N, e)"] A4 --> A6 A5 --> A7["المفتاح السري d"] B1["2. النص الصريح (رسالة) M"] --> B2["حساب C ≡ M^e (mod N)"] A6 -.-> B2 B2 --> B3["إرسال النص المشفر C"] B3 --> C1["3. النص المشفر المُستلَم C"] C1 --> C2["حساب M ≡ C^d (mod N)"] A7 -.-> C2 C2 --> C3["الحصول على النص الصريح الأصلي M"]

5.1 إنشاء المفتاح (Key Generation)

يقوم المتلقي بوب بإنشاء “مفتاح معلن” و"مفتاح سري" خاص به.

اختيار الأعداد الأولية: يختار عشوائيًا عددين أوليين كبيرين $p$ و $q$.
حساب المقياس $N$: يحسب $N = p \times q$. هذا الـ $N$ سيتم إعلانه للجميع.
حساب $\phi(N)$: يحسب دالة أويلر $\phi(N) = (p - 1)(q - 1)$. هذا هو رقم بوب السري.
اختيار المفتاح المعلن $e$: يختار عددًا صحيحًا $e$ بحيث $1 < e < \phi(N)$ ويكون أوليًا نسبيًا مع $\phi(N)$.
حساب المفتاح السري $d$: يجد عددًا صحيحًا $d$ يحقق الشرط التالي: $$ e \times d \equiv 1 \pmod{\phi(N)} $$ وهذا يعني، “العدد $d$ الذي يجعل باقي قسمة $e \times d$ على $\phi(N)$ يساوي $1$”.

وبهذا تكتمل تجهيزات المفاتيح.

المفتاح المعلن: الزوج $(N, e)$. يتم نشره للعالم أجمع.
المفتاح السري: $d$. لا يخبر به أحداً على الإطلاق.

5.2 التشفير (Encryption)

لنفترض أن أليس تريد إرسال رسالة سرية $M$ إلى بوب. ($M$ هي أحرف تم تحويلها إلى أرقام، بحيث يكون $0 \le M < N$). تستخدم أليس المفتاح المعلن لبوب $(N, e)$ للحساب كما يلي:

$$ C \equiv M^e \pmod N $$

تحسب “باقي قسمة رسالة $M$ مرفوعة للأس $e$ على $N$، وهو $C$”. هذا الـ $C$ هو النص المشفر.

5.3 فك التشفير (Decryption)

يتلقى بوب النص المشفر $C$. يستخدم بوب المفتاح السري $d$ للحساب كما يلي:

$$ M \equiv C^d \pmod N $$

عندما يحسب “باقي قسمة النص المشفر $C$ مرفوع للأس $d$ على $N$"، فإنه بشكل مذهل يسترجع الرسالة الأصلية $M$!

6. لماذا يعود إلى حالته الأصلية عند فك التشفير؟ (إثبات رياضي)

قد تتساءل: “كيف يمكن مجرد رفع $C$ للأس $d$ أن يعيده إلى $M$ الأصلي؟”. هنا تأتي قوة “مبرهنة أويلر” التي ذكرناها سابقًا.

دعونا نعوض معادلة التشفير $C = M^e$ في معادلة فك التشفير $C^d \pmod N$:

$$ C^d \equiv (M^e)^d \equiv M^{ed} \pmod N $$

هنا، تذكر الخطوة 5 من عملية إنشاء المفتاح. اختار بوب $d$ بحيث يكون $e \times d \equiv 1 \pmod{\phi(N)}$. هذا يعني أن “$ed$ هو من مضاعفات $\phi(N)$ مضافًا إليه $1$”. باستخدام عدد صحيح $k$، يمكن كتابتها على النحو التالي:

$$ ed = k \cdot \phi(N) + 1 $$

نعوض هذا في جزء الأس ونفككه باستخدام قوانين الأسس:

$$ M^{ed} = M^{k \cdot \phi(N) + 1} = M^{k \cdot \phi(N)} \times M^1 = (M^{\phi(N)})^k \times M $$

هنا، بافتراض أن الرسالة $M$ و $N$ أوليان نسبيًا، فإنه وفقًا لـ مبرهنة أويلر، يكون $M^{\phi(N)} \equiv 1 \pmod N$.

$$ (M^{\phi(N)})^k \times M \equiv 1^k \times M \equiv M \pmod N $$

وبالتالي، تتحقق المعادلة التالية بنجاح تام:

$$ C^d \equiv M \pmod N $$

أليس لا تعرف $d$، والمتنصت أيضًا لا يعرف $d$، لذا فإن الوحيد الذي يمكنه استخراج $M$ من $C$ هو بوب الذي يمتلك $d$.

7. مثال تطبيقي: دعونا نجرب تشفير RSA بالحساب اليدوي باستخدام أعداد أولية صغيرة

دعونا نقوم بتجربة تواصل مشفر من أليس إلى بوب باستخدام أرقام صغيرة (أعداد أولية) في الواقع.

【مرحلة إنشاء المفتاح لبوب】

يختار عددين أوليين $p=11$, $q=13$.
يحسب $N = 11 \times 13 = 143$.
يحسب $\phi(N) = (11 - 1) \times (13 - 1) = 10 \times 12 = 120$.
يختار مفتاحًا معلنًا $e$ يكون أوليًا نسبيًا مع $\phi(N)=120$. لنختر هنا $e=7$.
يجد المفتاح السري $d$. يبحث عن $d$ الذي يحقق $7 \times d \equiv 1 \pmod{120}$. في المعادلة $7d = 120k + 1$، عندما يكون $k=6$ يصبح الناتج $721$، و $721 \div 7 = 103$. وبالتالي، $d = 103$.

المفتاح المعلن: $(N=143, e=7)$
المفتاح السري: $d=103$

【مرحلة التشفير لأليس】 لنفترض أنها تريد إرسال الرسالة $M = 9$. المعادلة: $C \equiv 9^7 \pmod{143}$ $9^7 = 4,782,969$. بقسمة هذا العدد على 143 نحصل على $33447$ والباقي $48$. أصبح النص المشفر $C = 48$.

【مرحلة فك التشفير لبوب】 يتلقى بوب النص المشفر $C = 48$، ويستخدم المفتاح السري $d = 103$ لفك التشفير. المعادلة: $M \equiv 48^{103} \pmod{143}$ عند تنفيذ (48 ** 103) % 143 على الآلة الحاسبة، ستجد أن النتيجة بشكل رائع هي “9”! لقد تمكن من استلام الرسالة الأصلية بنجاح.

8. كيفية إيجاد المفتاح السري $d$: خوارزمية إقليدس الممتدة

في مثال الحساب اليدوي، اعتمدنا على التخمين للبحث عن $k$ وإيجاد $d=103$، لكن هذه الطريقة مستحيلة عندما تتكون الأعداد من مئات الأرقام. في البرامج الحقيقية، تُستخدم خوارزمية تسمى “خوارزمية إقليدس الممتدة”.

حل التطابق $7d \equiv 1 \pmod{120}$ هو نفسه إيجاد الأعداد الصحيحة $d, y$ التي تحقق المعادلة $7d + 120y = 1$. من خلال الحساب العكسي لخوارزمية إقليدس، يمكن إيجاد ذلك ميكانيكيًا.

$120 \div 7 = 17$ والباقي $1$
بتعديل هذه المعادلة نحصل على: $1 = 120 - 17 \times 7$
بمعنى آخر، $-17 \times 7 \equiv 1 \pmod{120}$

الرقم $-17$ يحمل نفس المعنى مثل $120 - 17 = 103$ في عالم المقياس $120$. وبالتالي، يمكن إيجاد $d = 103$ في لحظة. تتيح هذه الطريقة إجراء حسابات سريعة جدًا بغض النظر عن مدى ضخامة الأعداد.

9. الوجه الآخر لتشفير RSA: التوقيع الرقمي

من الميزات الرائعة لتشفير RSA أنه يمكن استخدامه أيضًا كـ “توقيع رقمي” عن طريق عكس أدوار المفتاح المعلن والمفتاح السري.

في التشفير كان الإجراء “التشفير بالمفتاح المعلن $\Rightarrow$ فك التشفير بالمفتاح السري”، أما في التوقيع الرقمي، نتبع خطوات “التشفير بالمفتاح السري $\Rightarrow$ فك التشفير بالمفتاح المعلن”.

flowchart TD A1["1. أليس تنشئ توقيعًا باستخدام مفتاحها السري"] --> A2["S ≡ M^d (mod N)"] A2 --> A3["إرسال الرسالة M والتوقيع S"] A3 --> B1["2. بوب يتحقق من التوقيع باستخدام المفتاح المعلن"] B1 --> B2["حساب M' ≡ S^e (mod N)"] B2 --> B3["التحقق مما إذا كان M' يتطابق مع M"]

تستخدم أليس مفتاحها السري $d$ لتحويل الرسالة (وهذا هو التوقيع $S$)، وترسله إلى بوب. يستخدم بوب المفتاح المعلن لأليس $e$ لإجراء عملية التحقق. إذا كانت نتيجة الحساب تتطابق مع الرسالة الأصلية، فإن ذلك يثبت في نفس الوقت أن “هذه البيانات لا يمكن إنشاؤها إلا بواسطة المفتاح السري لأليس”، و"لم يتم التلاعب بالرسالة في الطريق”.

10. تجربة تشفير RSA من خلال البرمجة

حسابات الأسس التي تعتبر صعبة في الحساب اليدوي، يمكن تنفيذها بسهولة بالغة باستخدام بايثون (Python). وفيما يلي كود بايثون يمكنك من خلاله تجربة المنطق الأساسي لتشفير RSA.

 1
 2
 3
 4
 5
 6
 7
 8
 9
10
11
12
13
14
15
16
17
18
19
20
21
22
23
24
25
26
27


def gcd(a, b):
 """إيجاد القاسم المشترك الأكبر"""
 while b != 0:
 a, b = b, a % b
 return a

def mod_inverse(e, phi):
 """إيجاد المفتاح السري d (باستخدام الدالة المدمجة في بايثون 3.8 وما بعده)"""
 return pow(e, -1, phi)

# 1. إنشاء المفتاح
p, q = 11, 13
N = p * q
phi = (p - 1) * (q - 1)
e = 7
d = mod_inverse(e, phi)

print(f"المفتاح المعلن: (N={N}, e={e}), المفتاح السري: d={d}")

# 2. التشفير
message = 9
ciphertext = pow(message, e, N)
print(f"النص المشفر: {ciphertext}")

# 3. فك التشفير
decrypted_message = pow(ciphertext, d, N)
print(f"الرسالة بعد فك التشفير: {decrypted_message}")

تستخدم دالة pow(base, exp, mod) في بايثون خوارزمية سريعة تسمى “التربيع المتكرر (Exponentiation by squaring)” داخليًا، لذلك حتى بالنسبة للأعداد المكونة من مئات الأرقام، ينتهي الحساب في لحظة.

11. الخلاصة وتقنيات التشفير المستقبلية

لقد قمنا بتوضيح آلية عمل تشفير RSA بناءً على المعرفة برياضيات المرحلة الثانوية.

صعوبة التحليل إلى العوامل الأولية: حساب $p \times q = N$ سهل، لكن إيجاد $p, q$ من $N$ صعب للغاية.
التطابق ومبرهنة أويلر: بفضل القاعدة $a^{\phi(N)} \equiv 1 \pmod N$، يكتمل “الباب الخلفي” السحري الذي يسمح بـ “العودة إلى الحالة الأصلية عند الرفع لأس عدد معين”.
المفتاح المعلن والمفتاح السري: يمكن لأي شخص التشفير، ولكن المستلم الشرعي فقط هو من يمكنه فك التشفير.

الـ $N$ المستخدم في تشفير RSA الحالي يتكون من أكثر من 600 رقم، وحتى لو تم حشد جميع الحواسيب العملاقة في العالم، فإن التحليل إلى العوامل الأولية سيستغرق وقتًا أطول من عمر الكون. ومع ذلك، عندما تصبح “الحواسيب الكمومية”، التي تتقدم الأبحاث بشأنها في السنوات الأخيرة، قابلة للتطبيق العملي في المستقبل، فهناك احتمال أن يتم حل هذا التحليل إلى العوامل الأولية في لحظة بواسطة “خوارزمية شور”. لذلك، يجري حاليًا تطوير “تشفير مقاوم للحواسيب الكمومية”، والذي لا يمكن حتى للحواسيب الكمومية فك تشفيره، بوتيرة سريعة في جميع أنحاء العالم.

الرياضيات المتقدمة التي غالبًا ما يُعتقد أنها “غير مفيدة”، تحمي في الواقع حياتنا اليومية من جذورها. يعتبر تشفير RSA من أفضل المواد التعليمية التي توضح لنا عمق وجمال هذه الرياضيات. نأمل أن تكون قد شعرت ولو قليلاً بمتعة التشفير والرياضيات من خلال هذا المقال.