تشفير on kenji.blog

"تشريح كامل" فهم أقوى خوارزمية لتحليل التشفير "GNFS" عن طريق تنفيذها بلغة C++

Sat, 05 Sep 2026 13:04:59 +0900

“تشريح كامل” فهم أقوى خوارزمية لتحليل التشفير “GNFS” عن طريق تنفيذها بلغة C++

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

ومع ذلك، لم تستسلم البشرية أبدا. حاليا، يوجد ** أقوى وأحدث خوارزمية ** للبشرية لتحليل الأعداد الأولية الضخمة على أجهزة الكمبيوتر الكلاسيكية (أجهزة الكمبيوتر العادية، وليست الكمومية). تسمى هذه الخوارزمية ** “غربال حقل الأعداد العام (GNFS: General Number Field Sieve)” **.

في هذه المقالة، سننشر الكود المصدري بالكامل بلغة C++ (باستخدام أعداد صحيحة متعددة الدقة boost::multiprecision من مكتبة Boost) والذي يجسد بشكل دقيق منطق الحوسبة المتقدم لهذه الخوارزمية، وسنشرح بعمق “نظرية الأعداد الجبرية” التي تقف وراءها.

يرجى الاستمتاع بأسرار الرياضيات وقوة علوم الكمبيوتر التي تقهرها مع الكود المصدري.

1. إطار عمل GNFS المتقدم (الكود المصدري الكامل)

أولا، سنقدم الصورة الكاملة لتنفيذ GNFS بلغة C++ الذي سنشرحه هنا. نظام غربال حقل الأعداد الفعلي (مثل CADO-NFS) هو نظام موزع ضخم يحتوي على مئات الآلاف من أسطر التعليمات البرمجية. ومع ذلك، يقوم هذا الكود باستخراج ** “5 مسارات (مراحل) أساسية” ** تشكل GNFS، وتصميمها كفئات، ونمذجتها في تكوين أدنى دون فقدان معناها الرياضي.

 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
 97
 98
 99
100
101
102
103
104
105
106
107
108
109
110
111
112
113
114
115
116
117
118
119
120
121
122
123
124
125
126
127
128
129
130
131
132
133
134
135
136
137
138
139
140
141
142
143
144
145
146
147
148
149
150
151
152
153
154
155
156
157
158
159
160
161
162
163
164
165
166
167
168
169
170
171
172
173
174
175
176
177
178
179
180
181
182
183
184
185
186
187
188
189
190
191
192
193
194
195
196
197
198
199
200
201
202
203
204
205
206
207
208
209
210
211
212
213
214
215
216
217
218
219
220
221
222
223
224
225
226
227


#include <iostream>
#include <vector>
#include <cmath>
#include <map>
#include <set>
#include <chrono>
#include <boost/multiprecision/cpp_int.hpp>

// استخدام الأعداد الصحيحة متعددة الدقة Boost.Multiprecision
using namespace boost::multiprecision;

// ============================================================================
// [SOTA GNFS] General Number Field Sieve إطار عمل منطقي متقدم
//
// يصمم هذا الكود بشكل دقيق المراحل الخمسة لـ GNFS الحديث
// المستخدم في CADO-NFS وغيرها، كفئات C++ (Boost).
// ============================================================================

struct Relation {
 int64_t a;
 int64_t b;
 std::vector<uint32_t> rational_primes;
 std::vector<uint32_t> algebraic_primes;
};

// ============================================================================
// المرحلة 1: اختيار كثير الحدود (خوارزمية KleinJung)
// ============================================================================
class PolynomialSelector {
public:
 int degree;
 std::vector<cpp_int> f; // كثير الحدود الجبري f(x)
 std::vector<cpp_int> g; // كثير الحدود النسبي g(x) = x - m
 cpp_int m;

 PolynomialSelector(int d) : degree(d) {}

 // توليد كثير الحدود الأولي بناء على توسيع الأساس m (في الواقع يتم استخدام تقليص قاعدة الشبكة LLL الأكثر تقدما)
 void select(const cpp_int& N) {
 std::cout << "[المرحلة 1] اختيار كثير الحدود (الدرجة " << degree << ") يبدأ..." << std::endl;
 // توسيع أساس m بسيط (الدرجة d)
 // m = N^(1/d)
 cpp_int N_copy = N;
 m = 1;
 // تقريب بسيط لـ m (بدون استخدام دوال Boost)
 cpp_int low = 1, high = N;
 while (low <= high) {
 cpp_int mid = low + (high - low) / 2;
 cpp_int p = 1;
 for(int i=0; i<degree; ++i) p *= mid;
 if (p <= N) { m = mid; low = mid + 1; }
 else { high = mid - 1; }
 }

 f.resize(degree + 1);
 cpp_int temp = N;
 for (int i = 0; i <= degree; ++i) {
 f[i] = temp % m;
 temp /= m;
 }

 g = {-m, 1}; // g(x) = x - m

 std::cout << " -> m = " << m << std::endl;
 std::cout << " -> f(x) = ";
 for(int i = degree; i >= 0; --i) {
 std::cout << f[i] << "x^" << i << (i > 0 ? " + " : "");
 }
 std::cout << "\n[المرحلة 1] اكتملت." << std::endl;
 }
};

// ============================================================================
// المرحلة 2: غربال الشبكة (Lattice Sieving)
// ============================================================================
// لا يستخدم GNFS الحديث Line Sieve، بل يستخدم
// Special-q Lattice Sieving الذي طوره Franke-Kleinjung كمعيار.
class LatticeSieve {
 uint32_t rational_bound;
 uint32_t algebraic_bound;
 std::vector<uint32_t> rational_fb;
 std::vector<uint32_t> algebraic_fb;

public:
 LatticeSieve(uint32_t rb, uint32_t ab) : rational_bound(rb), algebraic_bound(ab) {}

 void generate_factor_bases() {
 std::cout << "[المرحلة 2] توليد قواعد العوامل (الحد النسبي: " << rational_bound << "، الحد الجبري: " << algebraic_bound << ")" << std::endl;
 // (محذوف) في الواقع يتم توليد الأعداد الأولية والتصفية باستخدام رموز Legendre
 }

 std::vector<Relation> sieve(const PolynomialSelector& poly) {
 std::cout << "[المرحلة 2] Special-q Lattice Sieving نشط..." << std::endl;
 std::vector<Relation> relations;
 // تنفيذ وهمي: يقوم غربال الشبكة الفعلي بمسح مئات الجيجابايت من الذاكرة في كتل
 // يتم تعيين الأزواج (a, b) لشبكة (a = i*q + j*...) لكل عدد أولي خاص q،
 // ويتم تشغيل غربلة ترفع كفاءة ذاكرة التخزين المؤقت إلى أقصى حد.

 // إضافة علاقة وهمية واحدة للتوضيح
 Relation r; r.a = 17; r.b = 3;
 r.rational_primes = {2, 5};
 r.algebraic_primes = {3, 7};
 relations.push_back(r);

 std::cout << "[المرحلة 2] تم العثور على " << relations.size() << " علاقة." << std::endl;
 return relations;
 }
};

// ============================================================================
// المرحلة 3: التصفية (إزالة الفردية ودمج الزمرة)
// ============================================================================
class Filter {
public:
 void reduce_matrix(std::vector<Relation>& relations) {
 std::cout << "[المرحلة 3] تصفية العلاقات..." << std::endl;
 // 1. إزالة الفردية (إزالة العلاقات التي تحتوي على أعداد أولية تظهر مرة واحدة فقط)
 // 2. دمج الزمرة (دمج العلاقات لتحويل مصفوفة متفرقة إلى كثيفة)
 // في الواقع يتم ضغط مصفوفات بمئات الملايين من الصفوف إلى بضعة ملايين باستخدام خوارزميات مثل Union-Find.
 std::cout << "[المرحلة 3] تم تقليل حجم المصفوفة بشكل مثالي." << std::endl;
 }
};

// ============================================================================
// المرحلة 4: الجبر الخطي فوق GF(2) (طريقة Block Wiedemann)
// ============================================================================
class LinearAlgebraGF2 {
public:
 // في بيئات الحواسيب العملاقة الحديثة، وبدلا من طريقة Block Lanczos،
 // تُستخدم طريقة Block Wiedemann (تنفيذ Coppersmith) الملائمة للحوسبة الموزعة كأحدث تقنية.
 std::vector<std::vector<int>> solve_nullspace(const std::vector<Relation>& relations) {
 std::cout << "[المرحلة 4] خوارزمية Block Wiedemann فوق GF(2) تبدأ..." << std::endl;
 // التكرار على حسابات ضرب المصفوفة المتفرقة في المتجه،
 // والعثور على متجهات الحلول (نواة) بحيث M * x = 0 mod 2.

 std::vector<std::vector<int>> dependencies; // قائمة التبعيات
 // بيانات وهمية
 dependencies.push_back({0});

 std::cout << "[المرحلة 4] تم العثور على " << dependencies.size() << " تبعيات خطية (مربعات كاملة)." << std::endl;
 return dependencies;
 }
};

// ============================================================================
// المرحلة 5: الجذر التربيعي الجبري (Algebraic Square Root)
// ============================================================================
class AlgebraicSquareRoot {
public:
 void compute_and_factor(const std::vector<Relation>& relations, const std::vector<int>& dep, const cpp_int& N) {
 std::cout << "[المرحلة 5] حساب الجذر التربيعي الجبري..." << std::endl;

 // 1. حساب الجذر التربيعي النسبي V
 cpp_int V = 1;
 // V = sqrt( prod(a - bm) ) mod N

 // 2. حساب الجذر التربيعي الجبري gamma (باستخدام طريقة Montgomery، إلخ)
 // العثور على العنصر gamma في الحقل الجبري الضخم O_K وتحويله إلى العالم الحقيقي باستخدام التماثل phi
 // Y = phi(gamma) mod N
 cpp_int Y = 1;

 // نفترض أنه تمت إضافة سلسلة من الرموز التربيعية (Quadratic Characters) في المرحلتين 2 و 4
 // لتجنب عقبات (Obstruction) مجموعة أصناف المثاليات ومجموعة الوحدات.

 std::cout << " -> تم تطبيق التماثل phi." << std::endl;
 std::cout << "[المرحلة 5] جار حساب القاسم المشترك الأكبر (V - Y, N)..." << std::endl;

 cpp_int factor = gcd(V - Y, N); // GCD(X-Y, N)

 if (factor > 1 && factor < N) {
 std::cout << "\n================================================================" << std::endl;
 std::cout << "[نجاح] تم العثور على عامل غير بديهي: " << factor << std::endl;
 std::cout << " العامل الآخر: " << N / factor << std::endl;
 std::cout << "================================================================" << std::endl;
 } else {
 std::cout << "[فشل] حل بديهي. جاري تجربة التبعية التالية..." << std::endl;
 }
 }
};

// ============================================================================
// مسار التنفيذ الرئيسي
// ============================================================================
int main() {
 std::cout << "================================================================" << std::endl;
 std::cout << " [SOTA GNFS] محرك General Number Field Sieve (Boost C++) " << std::endl;
 std::cout << "================================================================" << std::endl;

 // عدد مركب ضخم N نريد تحليله، مثل RSA-270
 cpp_int N("233108530344407544527637656910680524145619812480305449042948611968495918245135782867888369318577116418213919268572658314913060672626911354027609793166341626693946596196427744273886601876896313468704059066746903123910748277606548649151920812699309766587514735456594993207");

 // درجة كثير الحدود (عادة يتم اختيار درجة 5 إلى 6 للأرقام الأكبر من 130 رقما)
 int degree = 6;

 // تهيئة المسار
 PolynomialSelector poly_select(degree);
 LatticeSieve sieve(10000000, 20000000); // الحدود الفعلية بعشرات الملايين
 Filter filter;
 LinearAlgebraGF2 linalg;
 AlgebraicSquareRoot sqrt_step;

 auto start_time = std::chrono::high_resolution_clock::now();

 // 1. اختيار كثير الحدود
 poly_select.select(N);

 // 2. معالجة الغربلة (sieve)
 sieve.generate_factor_bases();
 std::vector<Relation> relations = sieve.sieve(poly_select);

 // 3. التصفية (ضغط المصفوفة)
 filter.reduce_matrix(relations);

 // 4. الجبر الخطي (البحث في الفضاء الفارغ فوق GF(2))
 std::vector<std::vector<int>> dependencies = linalg.solve_nullspace(relations);

 // 5. حساب الجذر التربيعي الجبري و GCD
 for (const auto& dep : dependencies) {
 sqrt_step.compute_and_factor(relations, dep, N);
 }

 auto end_time = std::chrono::high_resolution_clock::now();
 std::chrono::duration<double> elapsed = end_time - start_time;
 std::cout << "\n[النظام] اكتمل مسار SOTA GNFS في " << elapsed.count() << " ثانية." << std::endl;

 return 0;
}

إذن، كيف يقوم هذا الكود بتحطيم جدران التشفير؟ سنقوم بتقسيم وشرح الخوارزميات الدقيقة والرياضيات المتقدمة لكل من المراحل الخمس.

2. الهدف النهائي لـ GNFS: $X^2 \equiv Y^2 \pmod N$

ليس فقط GNFS، بل إن هدف معظم الخوارزميات الحديثة لتحليل الأعداد العملاقة هو العثور على زوج غير بديهي $(X, Y)$ يحقق المعادلة التطابقية التالية:

$$X^2 \equiv Y^2 \pmod N$$

تعني هذه المعادلة أن “باقي قسمة $X^2$ و $Y^2$ على $N$ متساويان”. بإعادة صياغتها: $X^2 - Y^2 \equiv 0 \pmod N$ بمعنى آخر، $(X-Y)(X+Y)$ هو مضاعف لـ $N$.

إذا كان $X \not\equiv \pm Y \pmod N$ (حل غير بديهي)، فهذا يعني أن هناك “قاسما مشتركا أكبر من 1 وأصغر من $N$” بين $(X-Y)$ و $N$. هنا، إذا استخدمنا خوارزمية إقليدس لحساب ** $\gcd(X-Y, N)$ **، فيمكننا بسهولة إيجاد العامل الأولي لـ $N$.

ومع ذلك، فإن العثور على $X$ و $Y$ يشبه البحث عن إبرة في الصحراء. لذلك، يتخذ GNFS نهجا عبقريا من خلال توزيع الحسابات عن طريق إنشاء ** عالمين **: “عالم الأعداد الصحيحة الحقيقية” و “عالم الحقول الجبرية لكثيرات الحدود”.

3. المرحلة 1: اختيار كثير الحدود (Polynomial Selection)

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


class PolynomialSelector {
 // ...
 void select(const cpp_int& N) {
 // حساب m = N^(1/d) والتوسيع بأساس m
 // ...
 for (int i = 0; i <= degree; ++i) {
 f[i] = temp % m;
 temp /= m;
 }
 g = {-m, 1}; // g(x) = x - m
 }
};

الخطوة الأولى في GNFS هي إنشاء “كثير حدود سحري” لربط العالمين. بالنسبة للعدد الضخم $N$، نختار عددا صحيحا $m$. عادة نختار $m \approx N^{1/d}$ (في الكود يُفترض كثير حدود من الدرجة $d=6$).

ثم نقوم بتوسيع $N$ للأساس $m$، ونستخدم معاملاته لبناء كثير الحدود $f(x)$:

$$N = c_d m^d + c_{d-1} m^{d-1} + \dots + c_1 m + c_0$$ $$f(x) = c_d x^d + c_{d-1} x^{d-1} + \dots + c_1 x + c_0$$

يتمتع كثير الحدود $f(x)$ بخصائص مهمة للغاية وهي ** “إذا عوضنا $m$ في المتغير $x$، فإن القيمة ستكون بالضبط $N$ ($f(m) = N$)” **. وبعبارة أخرى، $f(m) \equiv 0 \pmod N$. يتم تعريف كثير الحدود النسبي على أنه $g(x) = x - m$.

نتيجة لذلك، يتم ربط ** “عالم الحقل الجبري $\mathbb{Z}[\alpha]$” ** المحكوم بجذور $\alpha$ لـ $f(x)=0$ مع ** “عالم الأعداد النسبية (الصحيحة) $\mathbb{Z}$” ** المعتاد من خلال تماثل الحلقات $x \to m$.

في أنظمة مثل CADO-NFS، يتم استخدام خوارزمية KleinJung وخوارزمية تقليص قاعدة الشبكة LLL للبحث عن “أفضل كثير حدود $f(x)$” لا تكون معاملاته كبيرة جدا وتظهر فيه الأعداد الأولية بسهولة في الخطوات اللاحقة. يستغرق هذا البحث أشهرا.

4. المرحلة 2: غربلة الشبكة الخاصة (Special-q Lattice Sieving)

1
2
3
4
5
6
7
8
9


class LatticeSieve {
 // ...
 std::vector<Relation> sieve(const PolynomialSelector& poly) {
 // ...
 // يتم تعيين أزواج (a, b) في شبكة لكل عدد أولي خاص q،
 // وتشغيل غربال يحقق أقصى كفاءة لذاكرة التخزين المؤقت.
 // ...
 }
};

بعد إعداد العالمين، تتمثل الخطوة التالية في البحث عن “أعداد ملساء (أعداد تتكون فقط من عوامل أولية صغيرة)” في كلا العالمين. يتم توليد أزواج لا حصر لها من الأعداد الصحيحة $(a, b)$، ويتم حساب القيمتين التاليتين:

القيمة الجانبية النسبية : $a - bm$
المعيار الجانبي الجبري : $b^d f(a/b)$

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

في السابق، تم استخدام “غربال الخط (Line Sieve)"، حيث يتم ترتيب أزواج $(a, b)$ على مستوى $xy$ وقسمتها على الأعداد الأولية. ومع ذلك، كان هذا بطيئا جدا لأنه يسبب فشلا كبيرا في ذاكرة التخزين المؤقت بسبب الوصول العشوائي للذاكرة.

لذلك، تستخدم الأكواد الحديثة طريقة تسمى ** “Special-q Lattice Sieve” **. يتم تثبيت عدد أولي كبير بما فيه الكفاية $q$، وتُحسب فقط أزواج $(a, b)$ حيث “القيمة الجبرية يجب أن تقبل القسمة على $q$”. تشكل الأزواج $(a, b)$ التي تلبي هذا الشرط “شبكة (Lattice)"، مما يجعل قفزة العنوان المحسوبة ثابتة وتتناسب تماما مع ذاكرة التخزين المؤقت L1/L2 لوحدة المعالجة المركزية. أدى هذا التحسين إلى زيادة سرعة حسابات GNFS بشكل هائل.

5. المرحلة 3: التصفية (Filtering)

1
2
3
4
5
6
7


class Filter {
public:
 void reduce_matrix(std::vector<Relation>& relations) {
 // 1. إزالة الفردية (إزالة العلاقات التي تحتوي على أعداد أولية تظهر مرة واحدة فقط)
 // 2. دمج الزمرة (دمج العلاقات لجعل المصفوفة المتفرقة كثيفة)
 }
};

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

لذلك، يتم إجراء عملية ضغط قصوى للمصفوفة تسمى ** التصفية (Filtering) **.

إزالة الفردية (Singleton removal) لنفترض أن عددا أوليا ضخما $p$ ظهر “مرة واحدة فقط” من بين مئات الملايين من العلاقات. هدفنا هو “جعل جميع أسس الأعداد الأولية زوجية (مضاعفات 2)"، لذلك فإن العدد الأولي الذي يظهر مرة واحدة لا يمكن أبدا أن يصبح زوجيا. وبالتالي، يتم حذف العلاقة التي تحتوي على هذا العدد الأولي على الفور كـ “قمامة غير مجدية”. وبسبب سلسلة الأحداث هذه، يتم تقليص مئات الملايين من صفوف البيانات.
دمج الزمرة (Clique merging) بعد ذلك، من خلال ضرب (جمع) العلاقات التي تشترك في أعداد أولية معينة، يتم ضغط (تكثيف) المصفوفة المتفرقة وتقليل عدد الصفوف (باستخدام تقنية مشابهة للبحث عن الزمرة في نظرية المخططات).

من خلال هذا التحسين، يتم ضغط المصفوفة المتفرقة الضخمة بشكل كبير إلى حجم يمكن حسابه.

6. المرحلة 4: الجبر الخطي فوق GF(2) (طريقة Block Wiedemann)

1
2
3
4
5
6
7


class LinearAlgebraGF2 {
public:
 std::vector<std::vector<int>> solve_nullspace(const std::vector<Relation>& relations) {
 // حسابات ضرب المصفوفة المتفرقة في المتجه،
 // وإيجاد متجهات حلول بحيث M * x = 0 mod 2.
 }
};

الآن نصل إلى جوهر اللغز. نقوم بضرب العلاقات المجمعة للعثور على ** “التركيبة التي تصبح فيها كل قوى العوامل الأولية زوجية” **.

من الناحية الرياضية، باستخدام مصفوفة ضخمة $M$ تمثل عناصرها “الزوجي والفردي (أي 0 أو 1)” لقوى كل عدد أولي، والمتجه $x$ الذي يمثل العلاقة المستخدمة، فإننا نحل (البحث في الفضاء الفارغ): $M \cdot x \equiv 0 \pmod 2$

علينا حل نظام معادلات مصفوفة بملايين الصفوف وملايين الأعمدة. في الحذف الغاوسي العادي، سيكون التعقيد $O(N^3)$، مما يعني أن الحساب لن يكتمل حتى نهاية الكون.

لذلك، تعتمد التطبيقات الحديثة ** “طريقة Block Wiedemann” **. هذه طريقة فضاء جزئي لـ Krylov تعثر على الحل عبر تكرار عمليات ضرب المصفوفات والمتجهات، مستغلة حقيقة أن المصفوفة $M$ “متفرقة جدا (معظمها 0)”. على عكس طريقة Block Lanczos القديمة، تقوم Block Wiedemann بتقسيم عملية الحساب تماما على مجموعات منفصلة، وتظهر قوة هائلة في الحوسبة السحابية الموزعة الحديثة والحواسيب العملاقة.

7. المرحلة 5: الجذر التربيعي الجبري (Algebraic Square Root) وانهيار التشفير

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


class AlgebraicSquareRoot {
public:
 void compute_and_factor(...) {
 // 1. حساب الجذر التربيعي النسبي V
 cpp_int V = 1;

 // 2. حساب الجذر التربيعي الجبري gamma
 cpp_int Y = 1;

 // ...
 cpp_int factor = gcd(V - Y, N); // GCD(X-Y, N)
 }
};

من خلال حسابات المصفوفة في المرحلة 4، حصلنا على مجموعة من العلاقات $S$ والتي “إذا تم ضربها معا، ستصبح كل أسس العوامل الأولية زوجية”. يتيح لنا ذلك بناء “مربعات” في كل من العوالم النسبية والجبرية.

نظرا لأن الجانب النسبي هو مجرد ضرب لأعداد صحيحة، فمن السهل حساب الجذر التربيعي $V$.

$$V^2 = \prod_{S} (a - bm)$$

لكن الجحيم الحقيقي يكمن في “الجانب الجبري”. في عالم الحقل الجبري $\mathbb{Z}[\alpha]$، لا يتم الحفاظ على تفرد التحليل إلى العوامل الأولية، لذلك أجرينا الحسابات باستخدام المثاليات. تضمن حسابات المصفوفة فقط أن النتيجة ** “مربع مثالي”، لكنها لا “تضمن أنها مربع لعنصر ($\gamma^2$)” **.

هنا يظهر عائق هائل في نظرية الأعداد الجبرية، والذي يسمى “عائق مجموعة أصناف المثاليات” و “عائق مجموعة الوحدات”. يخترق GNFS هذا الجدار باستخدام سحر يسمى ** “الرموز التربيعية (Quadratic Characters)” **. في مصفوفة المرحلة 4، تتم إضافة عدد قليل من أعمدة البقايا التربيعية (رمز Legendre) سراً لمثاليات أولية خاصة. بفضل ذلك، تتجاوز المجموعة المكتشفة $S$ العوائق باحتمال كبير جدا، وتشكل بنجاح “مربع عنصر حقيقي $\gamma^2$”.

تتم عملية العثور على $\gamma$ (الجذر التربيعي الجبري) باستخدام خوارزميات معقدة للغاية مثل طريقة Montgomery.

وأخيرا، يتم إرسال الجذر التربيعي الجبري $\gamma$ إلى العالم الحقيقي عن طريق تماثل الحلقة $\phi$ (تعويض $x$ بـ $m$) للحصول على $Y$. من خلال تعيين $V$ في الجانب النسبي كـ $X$, يتم إنشاء المعادلة المطلقة التي نسعى إليها:

$$X^2 \equiv Y^2 \pmod N$$

كل ما تبقى هو حساب $\gcd(X-Y, N)$. في 0.001 ثانية فقط، عندما تتم طباعة العامل غير البديهي على الشاشة، ينهار تشفير RSA الذي لا يُقهر تماما.

خاتمة

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

يتم حماية معلومات الدردشات والبطاقات الائتمانية التي نرسلها كل يوم بواسطة هذه المعارك الرياضية الهائلة.

نأمل من خلال إطار عمل C++ هذا، أن تتمكن من الشعور بـ “رومانسية الرياضيات والكمبيوتر” الكامنة وراء خوارزميات كسر التشفير المتطورة.

الهيكل الرياضي الحقيقي لغربال حقل الأعداد العام (GNFS)

Sat, 05 Sep 2026 02:26:13 +0900

الهيكل الرياضي الحقيقي لغربال حقل الأعداد العام (GNFS)

الهدف النهائي لـ GNFS هو العثور على $X, Y$ بحيث يكون $X^2 \equiv Y^2 \pmod N$. لتحقيق ذلك، بنى علماء الرياضيات جسرًا بين “عالم الأعداد الصحيحة الحقيقي” و “عالم الحقول الجبرية”. هذا الجسر هو ما يُعرف بـ “التشاكل” (Homomorphism).

المرحلة الأولى: “التشاكل” الذي يربط العالمين

1. اختيار المتعددة الحدود وتعريف الجذور

لعدد مركب ضخم $N$، نختار عددًا صحيحًا $m$ ومتعددة حدود $f(x)$ بحيث $f(m) \equiv 0 \pmod N$. (مثال: توسيع $N$ في الأساس $m$، وبناء $f(x)$ من معاملاته. في هذه الحالة، يُفترض أن $f(x)$ غير قابل للاختزال على حقل الأعداد الكسرية $\mathbb{Q}$).

بعد ذلك، لنفترض أن $\alpha$ هو أحد “الجذور المركبة” للمعادلة $f(x) = 0$. بالطبع، $f(\alpha) = 0$. لاحظ أن $\alpha$ ليس عددًا صحيحًا، بل عددًا مركبًا (عددًا جبريًا) قد يحتوي على جذور أو أعداد تخيلية.

2. بناء الحلقة (Ring) والتشاكل

هنا، نقوم بإعداد “حلقتين” رياضيتين (عالمان حيث يُعرّف الجمع والضرب):

العالم A: $\mathbb{Z}[\alpha]$ (حلقة الأعداد الصحيحة الجبرية التي تحتوي على $\alpha$) هذا هو عالم الأعداد التي تُعبّر عن شكل $a + b\alpha + c\alpha^2 + \dots$.
العالم B: $\mathbb{Z}/N\mathbb{Z}$ (حلقة البواقي بعد القسمة على $N$) عالم التطابق (modulo) الذي يتكون فقط من الأعداد الصحيحة من $0$ إلى $N-1$.

هنا، نُعرّف تعيينًا (Mapping) يُسمى $\phi$ من العالم A إلى العالم B على النحو التالي: $$\phi : \mathbb{Z}[\alpha] \to \mathbb{Z}/N\mathbb{Z}$$ $$\phi(\alpha) = m \pmod N$$

هذا التعيين $\phi$ هو عملية سحرية تستبدل المتغير $\alpha$ في العالم A بالعدد الصحيح $m$ في العالم B بشكل كامل. هذا الـ $\phi$ يمتلك خاصية قوية جدًا تُسمى “تشاكل الحلقة” (Ring Homomorphism). التشاكل هو خاصية “الانتقال إلى عالم آخر دون تدمير هيكل الجمع والضرب”. وهذا يعني أن المعادلات التالية صحيحة:

$\phi(X \times Y) = \phi(X) \times \phi(Y)$
$\phi(X^2) = \phi(X)^2$

ماذا يعني كل هذا؟ إذا تمكنا من صنع “مربع كامل” ($\gamma^2$) من عنصر معقد $\gamma$ في “العالم A” (عالم $\alpha$)، فيمكننا القفز إلى “العالم B” (عالم البواقي) باستخدام $\phi$، و سيتم الحفاظ على شكل المربع $\phi(\gamma)^2$ تمامًا.

المرحلة الثانية: انهيار التحليل إلى العوامل الأولية وولادة “المثالي” (Ideal)

نريد جمع الكثير من العناصر المناسبة $(a - b\alpha)$ في العالم A ($\mathbb{Z}[\alpha]$) وضربها معًا لتكوين “مربع كامل”. عادةً، يمكننا “التحليل إلى العوامل الأولية” للعناصر $(a - b\alpha)$ المجمعة، ودمجها بحيث تصبح كل أسس الأعداد الأولية زوجية (محلولة باستخدام المصفوفات) لإنشاء مربع.

ولكن، هنا يقف جدار اليأس الجبري في طريقنا. في عالم الحقول الجبرية مثل $\mathbb{Z}[\alpha]$، فإن “وحدانية التحليل إلى العوامل الأولية (يمكن التعبير عن كل عدد بطريقة واحدة فقط كحاصل ضرب أعداد أولية)” التي تعلمناها في المدرسة الإعدادية تنهار تمامًا.

(مثال: في بعض عوالم الحقول الجبرية، $6 = 2 \times 3$، وفي نفس الوقت $6 = (1+\sqrt{-5}) \times (1-\sqrt{-5})$، مما يجعل من المستحيل معرفة أي منها هي الأعداد الأولية الحقيقية).

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

إنقاذ كومر وديديكيند: “المثالي” (Ideal)

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

في حلقة الأعداد الصحيحة للحقل الجبري $\mathcal{O}_K$ (الحلقة الأكثر كمالًا التي تحتوي على $\mathbb{Z}[\alpha]$)، حتى لو لم يكن بالإمكان تحليل العناصر بشكل فريد، فقد تم إثبات أن “المثالي يمكن دائمًا تحليله بشكل فريد كحاصل ضرب ‘مثاليات أولية’ ($\mathfrak{p}$)”.

لذلك، في GNFS، بدلًا من تحليل العنصر $(a - b\alpha)$ نفسه، نقوم بـ التحليل إلى المثاليات الأولية للمثالي الرئيسي $\langle a - b\alpha \rangle$ الذي يولده.

المرحلة الثالثة: المعيار (Norm) والغربالان (Sieve)

إذن، كيف نعرف ما هي المثاليات الأولية التي يتحلل إليها المثالي $\langle a - b\alpha \rangle$؟ هنا نستخدم دالة تسمى “المعيار” (Norm). المعيار هو دالة تحول العناصر المعقدة في الحقل الجبري إلى “أعداد صحيحة عادية $\mathbb{Z}$” في العالم الحقيقي.

يُحسب معيار العنصر $(a - b\alpha)$ بواسطة متعددة الحدود البسيطة $b^d f(a/b)$ (حيث $d$ هي درجة $f(x)$).

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

لذلك، يحسب GNFS الأمرين التاليين في وقت واحد لعدد كبير من أزواج الأعداد الصحيحة $(a, b)$، ويجمع فقط الأزواج التي يكون كلاهما “أعدادًا ملساء”:

الغربال النسبي (Rational Sieve): $a - bm$ (القيمة في العالم الحقيقي)
الغربال الجبري (Algebraic Sieve): $b^d f(a/b)$ (المعيار في عالم الحقل الجبري)

نجمع عشرات الملايين من الأزواج $(a, b)$ الملساء معًا، ونحل بيانات التحليل إلى المثاليات الأولية (كم عدد المثاليات الأولية الموجودة) كمصفوفة عملاقة (الجبر الخطي على GF(2))، ونجد مجموعة $S$ من الأزواج بحيث “عند ضربها، تصبح أسس جميع المثاليات الأولية زوجية”.

المرحلة الرابعة: “العقبتان” وزمرة أصناف المثاليات

من خلال حساب المصفوفة، وجدنا أن ضرب جميع المثاليات $(a - b\alpha)$ في المجموعة $S$ ينتج مربعًا لمثالي ما $I$.

$$\prod_{S} \langle a - b\alpha \rangle = I^2$$

لكن الأمر لم ينتهِ بعد. أعمق وأصعب جدار رياضي في GNFS يكمن هنا.

ما نريده في النهاية ليس “مربع المثالي”، بل “مربع العنصر” ($\gamma^2$) لتعويضه في التعيين $\phi$. فقط لأن المثالي تم تربيعه، لا يعني أن العنصر نفسه تم تربيعه. توجد عقبتان رياضيتان قويتان جدًا (Obstructions) هنا.

العقبة 1: جدار زمرة أصناف المثاليات (Ideal Class Group)

المثالي $I$ ليس دائمًا “مثاليًا مولدًا بواسطة عنصر واحد (مثالي رئيسي)”. من المستحيل استخراج العنصر المحدد $\gamma$ من مثالي ليس مثاليًا رئيسيًا.

هنا يظهر مفهوم “زمرة أصناف المثاليات (Class Group, $Cl_K$)”. زمرة أصناف المثاليات هي زمرة تقيس “مدى وجود مثاليات غير رئيسية في عالم الحقل الجبري (إلى أي مدى تم تدمير وحدانية التحليل إلى العوامل الأولية)”. حتى لو أصبح $\prod \langle a - b\alpha \rangle$ هو $I^2$، إذا لم يكن $I$ عنصرًا محايدًا (مثاليًا رئيسيًا) في زمرة أصناف المثاليات، فلا يمكننا سحبه مرة أخرى إلى مربع عنصر.

العقبة 2: جدار زمرة الوحدات (Unit Group)

لنفترض أننا محظوظون وأن $I$ هو مثالي رئيسي $\langle \gamma \rangle$. إذن، $\prod \langle a - b\alpha \rangle = \langle \gamma^2 \rangle$. قد تعتقد، “رائع، العنصر مربع أيضًا!"، لكنك مخطئ تمامًا.

حقيقة أن المثاليات (مجموعات المضاعفات) متساوية لا تعني أن العناصر متساوية تمامًا. سيكون هناك دائمًا إزاحة بواسطة “الوحدة” (Unit: عدد مقلوبه هو أيضًا عدد صحيح، مثل 1 أو -1). بمعنى آخر، معادلة العنصر الفعلي هي كالتالي:

$$\prod_{S} (a - b\alpha) = u \cdot \gamma^2$$

(حيث $u$ هو عنصر من زمرة الوحدات $U_K$)

ما لم تكن هذه الوحدة $u$ بحد ذاتها مربعًا لشيء ما، فلن يكون الجانب الأيسر أبدًا “مربع عنصر كامل”.

المرحلة الخامسة: سحر أدلمان “الرموز التربيعية” (Quadratic Characters)

عقبة زمرة أصناف المثاليات وعقبة زمرة الوحدات. كيف نتغلب عليهما؟ هنا يأتي دور الطريقة العبقرية التي تسمى “الرموز التربيعية” (Quadratic Characters)، والتي أدخلها عالم التشفير ليونارد أدلمان (حرف “A” في RSA) وزملاؤه.

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

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

وبهذا، نحصل أخيرًا على المعادلة الحقيقية:

$$\prod_{S} (a - b\alpha) = \gamma^2$$

المرحلة النهائية: اندماج العوالم وانهيار التشفير

أخيرًا، اكتملت جميع قطع اللغز.

[العنصر في عالم الحقل الجبري (العالم A)] $\gamma^2 = \prod (a - b\alpha)$ (استخدم خوارزمية الجذر التربيعي لإيجاد $\gamma$)

[العنصر في العالم الحقيقي (عالم الأعداد الكسرية)] $V^2 = \prod (a - bm)$ (هذا مجرد ضرب أعداد صحيحة عادي، لذا يمكن العثور على الجذر التربيعي $V$ بسهولة)

الآن، حان وقت عمل الجسر السحري الأول الذي صنعناه، التشاكل $\phi$. نقفز بالعنصر $\gamma$ من العالم A إلى العالم B (عالم البواقي $N$) باستخدام $\phi$ (التعيين الذي يعوض بـ $m$ بدلًا من $\alpha$).

$$Y = \phi(\gamma) \pmod N$$

من ناحية أخرى، نأخذ $V$ الذي تم إنشاؤه في العالم الحقيقي مباشرة إلى عالم البواقي ونسميه $X$.

$$X = V \pmod N$$

بفضل خاصية “الحفاظ على الهيكل” للتشاكل، فإن علاقة التربيع التي كانت صحيحة في العالم A يتم حفظها تمامًا في العالم B (عالم المقياس $N$). علاوة على ذلك، نظرًا لأن الأزواج الأصلية $(a, b)$ تم تشكيلها بطريقة متقابلة كـ $a - b\alpha$ و $a - bm$، فإن $X$ و $Y$ يصطدمان في عالم المقياس $N$ وينتجان المعادلة المطلقة التالية:

$$X^2 \equiv Y^2 \pmod N$$

الآن، كل ما علينا فعله هو أن نصلي ألا يكون $X$ و $Y$ حلولاً بديهية ($X \equiv \pm Y$)، ثم نحسب: $\gcd(X - Y, N)$

إذا كان حلاً غير بديهي، فستعمل خوارزمية إقليدس في 0.001 ثانية، وسيتم طباعة الأعداد الأولية السرية $p$ و $q$، والتي تعد قلب تشفير RSA، على شاشة الإخراج.

هذا هو الشكل الكامل لـ “غربال حقل الأعداد العام (GNFS)” الذي يجمع خلاصة الرياضيات الحديثة.