Tài liệu Kỹ thuật lập trình - Chương 2: Giới thiệu lý thuyết số

Chia sẻ: | Ngày: | Loại File: DOC | Số trang:20

Thêm vào BST

Báo xấu

98
lượt xem 15
download

Download Vui lòng tải xuống để xem tài liệu đầy đủ

Trong chương này sẽ trình bày ngắn gọn những kiến thức cơ bản về lý thuyết số, nó là công cụ hữu ích để hiểu sâu một thuật toán mật mã nào đó.

Chủ đề:

Bình luận(0) Đăng nhập để gửi bình luận!

Lưu

Nội dung Text: Tài liệu Kỹ thuật lập trình - Chương 2: Giới thiệu lý thuyết số

CHƯƠNG 2 GIỚI THIỆU VỀ LÝ THUYẾT SỐ Hầu hết các thuật toán mật mã khóa công khai được xây dựng dựa trên các khái ni ệm c ủa lý thuyết số. Trong chương này sẽ trình bày ngắn gọn những kiến thức cơ b ản v ề lý thuy ết s ố, nó là công cụ hữu ích để hiểu sâu một thuật toán mật mã nào đó. 2.1 Các số nguyên tố và các số nguyên tố cùng nhau 2.1.1 Các ước số Nói rằng, b (một số khác 0) là ước số của a nếu a = mb, với một giá trị m nào đó, ở đây a, b, m là các số nguyên. Như vậy, b là ước số của a, nếu như chia a cho b không còn lại số dư. Để kí hiệu b là ước số của a thường sử dụng cách viết ba. Có các quan hệ sau: o Nếu a1, thì a = ± 1. o Nếu ba và ab, thì a = ± b. o Số b bất kỳ khác 0 là ước số của 0. o Nếu bg và bh, thì b(mg + nh) đối với bất kì số nguyên m, n. Để chứng minh khẳng định cuối cùng, cần chú ý như sau: nếu bg thì g có dạng g = b * g1, đối với số nguyên g1 nào đó, nếu bh thì h có dạng h = b * h1, đối với giá trị nguyên h1 nào đó, Vì vậy: mg + nh = mbg1 + nbh1 = b * (mg1+ nh1) Cuối cùng, b là ước số của mg + nh. Ví dụ 2.1: Các số 1, 2, 3, 4, 6, 8, 12 và 24 là các ước số của 24. Ví dụ 2.2: b = 7, g = 14, h = 63, m = 3, n = 2 714 và 763. Yêu cầu chứng minh rằng 7(3 * 14 + 2 * 63). Chúng ta có: ( 3 * 14 + 2 * 63 ) = 7(3 * 2 + 2 * 9) Như vậy, rõ ràng rằng 77(3 * 2 + 2 * 9). 2.1.2 Các số nguyên tố Một số nguyên p > 1 được gọi là số nguyên tố, nếu chỉ có ± 1 và ± p là ước số của nó. Một số nguyên bất kỳ a > 1 có thể phân tích thành các thừa số và được trình bày dưới dạng: a = p1 1 p2 2 ... ptαt , α α ở đây: p1 < p2 < … < pt là các số nguyên tố, còn các giá trị αi > 0. Ví dụ 2.3: 91 = 7 * 13; 11011 = 7 * 112 * 13. Nếu P là kí hiệu tập hợp tất cả các số nguyên tố, thì đối với m ột s ố nguyên d ương b ất kì đ ược viết duy nhất dưới dạng: a a = ∏ p p, ở đây tất cả các ap ≥ 0 P 1
Trong công thức này, biểu thức ở vế phải sau dấu bằng, ký hiệu tích theo t ất cả khả năng của các số nguyên tố p. Đối với mỗi giá trị cụ thể của a thì giá trị lớn nhất của chỉ số ap sẽ bằng 0. Ví dụ 2.4: 3600 = 24 * 32 * 52. Các giá trị của một số nguyên dương bất kỳ có thể liệt kê dưới một dạng đơn giản của t ất cả các chỉ số khác không theo công thức ở trên. Ví dụ 2.5: Số nguyên 12 có thể trình bày như {a2 = 2, a3 = 1}. Số nguyên 18 có thể trình bày như {a2 = 1, a3 = 2}. Phép nhân hai số nguyên tương đương với phép cộng các giá trị các chỉ số phù hợp: k = m * n → kp = mp + np, đối với tất cả các p. Ví dụ 2.6: k = 12 * 18 = 216, k2 = 2 + 1 = 3, k3 = 1 + 2 = 3, 216 = 23 * 33 Bổ đề: Một số nguyên dương bất kì dạng pk chỉ có thể chia hết cho một số nguyên, khi số bị chia có bậc của số nguyên tố p với chỉ số không vượt hơn k, nghĩa là số pj với j ≤ k. Như vậy: ab → ap ≤ bp, đối với tất cả p. Ví dụ 2.7: a = 12 , b = 36 , 1236, 12 = 22 * 3, 36 = 22 * 32 a2 = 2 = b2, a3 = 1 ≤ 2 = b3. 2.1.3 Các số nguyên tố cùng nhau Chúng ta sẽ sử dụng ký hiệu gcd( a, b) để chỉ ước số chung lớn nhất (UCLN) của số a và b. Nói rằng, một số nguyên dương c là UCLN của a và b, nếu: o c là ước số của a và b. o Ước số bất kỳ của a và b đều là ước số của c. Có thể định nghĩa tương đương như sau: gcd(a, b) = max [k, khi ka và kb]. Bởi vì, đòi hỏi rằng UCLN là một số dương, chúng ta có gcd( a, b) = gcd(a, –b ) = gcd(–a, b) = gcd(–a, –b ). Nói chung, gcd(a, b) = gcd(|a|, |b|). Ví dụ 2.8: gcd(60, 24 ) = gcd(60, –24 ) = 12. Bởi vì tất cả các số nguyên khác không đều là ước số của số 0, chúng ta luôn luôn có: gcd( a, 0 ) = | a|. Dễ dàng xác định được UCLN của hai số nguyên dương, nếu các số này được trình bày dưới dạng tích của các thừa số nguyên tố. Ví dụ 2.9: 300 = 22 * 31 * 52, 18 = 21 * 32. gcd(18, 300) = 21 × 31 × 50 = 6. 2
Trong trường hợp chung: k = gcd(a, b) → kp = min (ap, bp), đối với tất cả các p. Việc xác định các thừa số nguyên tố của các số lớn là bài toán không đ ơn gi ản, b ởi vì r ằng các trình bày ở trên không cho một khả năng thực tiễn để tính UCLN của hai số. Các số nguyên a và b là nguyên tố cùng nhau, nếu chúng không có ước số nguyên t ố chung, hay ước số chung duy nhất của chúng là 1. Nói một cách khác a và b là hai số nguyên tố cùng nhau nếu gcd(a, b) = 1. Ví dụ 2.10: Số 8 và số 15 là các số nguyên tố cùng nhau, bởi vì ước s ố của 8 là 1, 2, 4 và 8, còn các ước số của 15 là 1, 3, 5 và 15. Như vậy, 1 là ước số chung duy nhất của hai số này. 2.1.4 Số học trong lớp số dư Đối với bất kỳ một số nguyên dương n và một số nguyên a bất kỳ, khi chia a cho n, ta nhận được một số nguyên q nào đó và số dư r, phù hợp với quan hệ sau: a = qn + r, 0 ≤ r < n, q = a / n. ở đây x ký hiệu số nguyên lớn nhất, không lớn hơn x. Đối với một số dương a và một số nguyên dương n, luôn luôn có thể tìm được q và r, phù hợp với quan hệ tính toán trên. Ví dụ 2.11: a = 11; n = 7; 11 = 1 * 7 + 4; r = 4. Nếu a là một số nguyên, còn n là một số nguyên dương, thì a mod n, được xác định như phần dư của phép chia a cho n. Như vậy, đối với một số nguyên a bất kỳ có thể viết : a = a / n * n + (a mod n) Ví dụ 2.12: 11 mod 7 = 4 ; –11 mod 7 = 3. 2.2 Lý thuyết về đồng dư Định nghĩa 2.1 (Đồng dư) Cho a, b ∈ Z , n ∈ N . Ta nói a đồng dư với b theo modulo n, khi a và b chia cho n có cùng số dư, và được viết dưới dạng sau: a ≡ b mod n . Chứng minh: Giả sử chia a và b cho n và thu được các thương nguyên và phần dư. Các phần dư nằm giữa 0 và n – 1, nghĩa là a = q1n + r1 và b = q2 n + r2 . Trong đó 0 ≤ r1 ≤ n − 1 và 0 ≤ r2 ≤ n − 1 . Khi đó có thể dễ dàng thấy rằng a ≡ b mod n khi và chỉ khi r1 = r2 . Như vậy: a ≡ b mod n khi và chỉ khi a mod n = b mod n . Chúng ta đi tìm hiểu một số tính chất quan trọng của đồng dư. Tính chất 2.1: Nếu a1 ≡ b1 mod n và a2 ≡ b2 mod n thì a1 + a 2 ≡ b1 + b2 mod n Chứng minh: Từ giả thuyết ta có a1 = nt1 + b1 và a2 = nt 2 + b2 , suy ra a1 + a2 =b1 +b2 + (t1 + t2 )n , điều này có nghĩa là a1 + a 2 ≡ b1 + b2 mod n . Tính chất 2.2: Nếu a1 ≡ b1 mod n và a2 ≡ b2 mod n thì a1 *a 2 ≡ b1 * b2 mod n 3
Chứng minh: Từ giả thuyết ta có a1 = nt1 + b1 và a2 = nt 2 + b2 , suy ra a1 * a2 =b1 *b2 + (b1t 2 + b2 t1 + nt1t 2 ) n , điều này có nghĩa là a1 *a 2 ≡ b1 * b2 mod n Tính chất 2.3: Nếu như gcd(d , n) = 1, a ≡ b mod n, a = a1d , b = b1d , thì a1 ≡ b1 mod n Chứng minh: Từ giả thuyết ta có (a1 − b1 ) d | n , mà gcd(d , n) = 1 , nên suy ra (a1 − b1 ) | n , hay a1 ≡ b1 mod n . Ví dụ 2.13: Chứng minh rằng 37n+2+16n+1+23n chia hết cho 7. Giải: Rõ ràng chúng ta thấy rằng 37 ≡ 2 mod 7 , 16 ≡ 2 mod 7 và 23 ≡ 7 mod 7 , từ tính chất 2.2 có: 37 n+ 2 ≡ 2 n+ 2 mod 7 , 16 n+1 ≡ 2 n+1 mod 7 , 23n ≡ 7 n mod 7 . Từ tính chất 2.1 có: 37n+2+16n+1+23n ≡ 7.2 n mod 7 , và từ tính chất 2.3 suy ra điều phải chứng minh. 2.3 Các số nguyên modulo n Các số nguyên modulo n (ký hiệu Z n ) là tập hợp các số nguyên { 0, 1, ..., n − 1} bao gồm 2 phép toán cộng và nhân. Việc cộng và nhân trong Z n được thực hiện giống như cộng và nhân các s ố nguyên, ngoại trừ một điểm là các kết quả sẽ được rút gọn theo modulo n. Ví dụ 2.14: Tính 11 * 13 trong Z16 . Tương tự như với các số nguyên ta có 11 * 13 = 143 . Để rút gọn 143 theo modulo 16, ta thực hiện phép chia bình thường: 143 = 8 × 16 + 15 , bởi vậy 143 mod16 = 15 . Do đó 11 × 13 = 15 trong Z16 . Các định nghĩa trên phép cộng và phép nhân trong Z n thoả mãn hầu hết các quy tắc quen thuộc trong số học. Sau đây ta sẽ liệt kê mà không chứng minh các tính chất này: o Phép cộng là đóng, tức với bất kì a, b ∈ Z n , a + b ∈ Z n . o Phép cộng là giao hoán, tức là với bất kì a, b ∈ Z n thì: a + b = b + a. o Phép cộng là kết hợp, tức với bất kì a, b, c ∈ Z n : (a + b) + c = a + (b + c). o 0 là phần tử đơn vị của phép cộng, có nghĩa là với bất kì a ∈ Z n : a + 0 = 0 + a = a. o Phần tử nghịch đảo của phép cộng của một phần tử bất kì a ∈ Z n là m – a, nghĩa là a + (m – a) = (m – a) + a = 0 với bất kì a ∈ Z n . o Phép nhân là đóng, tức với bất kì a, b ∈ Z n , a * b ∈ Z n . o Phép nhân là giao hoán, nghĩa là với bất kì a, b ∈ Z n , a * b = b * a . o Phép nhân là kết hợp, nghĩa là với a, b, c ∈ Z n , (a * b) * c = a * (b * c) . o 1 là phần tử đơn vị của phép nhân, tức là với bất kì a ∈ Z n : a *1=1* a = a . 4
o Phép nhân có tính chất phân phối đối với phép cộng, tức là đối với a, b, c ∈ Z n , (a + b) * c = a * c + b * c và a * (b + c) = a * b + a * c . 2.4 Hàm Euler, định lý Euler và định lý Fermat Định nghĩa 2.2 (Hàm Euler) Cho n là số nguyên dương, đặt φ (n) là số các phần tử của tập hợp, mà tập này là các số nguyên trong khoảng [1, n] và nguyên tố cùng nhau với n, thì φ (n) gọi là hàm Euler. Ta công nhận một số tính chất quan trọng của hàm Euler: 1. φ (1) = 1 2. Nếu p là số nguyên tố thì φ ( p) = p − 1 3. Nếu như p và q là hai số nguyên tố cùng nhau thì φ ( pq) = φ ( p ) * φ ( q) 4. φ ( p s ) = p s −1 ( p − 1) e e e 5. Nếu như n = p11 p22 ... pk k , ở đây p1 , p2 , …, pk là số nguyên tố, thì  1  1   1  φ (n) = n1 − 1 − ...1 −    p1    p2   pk    Trên bảng 2.1 trình bày 30 giá trị đầu tiên của φ (n) . Giá trị φ(1) là không xác định, nhưng coi rằng nó bằng 1. Bảng 2.1: Một số giá trị của hàm Euler φ (n) n φ(n) n φ(n) n φ(n) 1 1 11 10 21 12 2 1 12 4 22 10 3 2 13 12 23 22 4 2 14 6 24 8 5 4 15 8 25 20 6 2 16 8 26 12 7 6 17 16 27 18 8 4 18 6 28 12 9 6 19 18 29 28 10 4 20 8 30 8 Định lý Euler: Cho n > 1, gcd(a, n) = 1 , và φ (n) là hàm Euler. Khi đó ta có: a φ ( n ) ≡ 1 mod n Ví dụ 2.15: a = 3, n = 10, φ(10) = 4, 34 = 81 ≡ 1 mod 10, 5
a = 2, n = 11, φ(11) = 10, 210 = 1024 ≡ 1 mod 11. Chứng minh: Giả sử x1 , ..., xφ ( n ) – là các số tự nhiên khác nhau, nhỏ hơn n và nguyên tố cùng nhau với n. Hãy xét tất cả các khả năng của tích xi a , với i ∈ 1, φ ( n) . Bởi vì a nguyên tố cùng nhau với n và xi nguyên tố cùng nhau với n, nên tích xi a cũng nguyên tố cùng nhau với n, do đó có xi a ≡ x j mod n . Chú ý rằng các phần dư của phép chia xi a cho n là khác nhau. Nếu điều này không đúng, có nghĩa là tồn tại i1 ≠ i2 , sao cho xi1 a ≡ xi 2 a mod n Cho nên ( xi1 − xi2 ) a ≡ 0 mod n . Bởi vì a nguyên tố cùng nhau với n, nên biểu thức cuối cùng tương đương với xi1 − xi2 ≡ 0 mod n Điều này là mâu thuẫn, bởi các số x1 , ..., xφ ( n ) là các cặp khác nhau theo modulo n. Hãy nhân tất cả đẳng thức xi a ≡ x j mod n , thì thu được: x1... xφ ( n ) a φ ( n ) ≡ x1... xφ ( n ) mod n Hay x1...xφ ( n ) ( aφ ( n ) − 1) ≡ 0 mod n Bởi vì số x1... xφ ( n ) mod n nguyên tố cùng nhau với n nên đẳng thức cuối cùng tương đương với a φ ( n ) − 1 ≡ 0 mod n hay a φ ( n ) ≡ 1 mod n φ (n ) + 1 Ta có một công thức quan trọng sau: a ≡ a mod n. Định lý Fermat nhỏ: Cho p là số nguyên tố, a là số nguyên dương không chia hết cho p. Khi đó ta có a p −1 ≡ 1 mod p . Chứng minh: Ta có φ ( p) = p − 1 , áp dụng định lý Euler ta có điều phải chứng minh. Từ định lý Fermat chúng ta có các hệ quả quan trọng sau: 1. Cho a ∈ Z , p là số nguyên tố, thì ta có: a p ≡ a mod p 2. Cho a, b ∈ Z , p là số nguyên tố (a + b) n ≡ a n + b n mod p Ví dụ 2.16: Chứng mình rằng 118 +218 +318 +418 +518 +618 ≡ –1 mod 7 Giải: Các số 1, 2, 3, 4, 5, 6 là số nguyên tố cùng nhau v ới 7. Nên theo đ ịnh lý Fermat chúng ta có các phương trình sau: 16 ≡ 1 mod 7 ; 2 6 ≡ 1mod 7 ; 36 ≡ 1mod 7 ; 4 6 ≡ 1 mod 7 ; 56 ≡ 1 mod 7 ; 66 ≡ 1 mod 7 Lấy bậc ba hai vế từng phương trình và cộng lại ta được: 118 +218 +318 +418 +518 +618 ≡ 6 mod 7 = −1 mod 7 6
2.5 Thuật toán Euclide và thuật toán Euclide mở rộng 2.5.1 Thuật toán Euclide Định lý Euclide: Cho a, b ∈ Z , b ≠ 0 , tồn tại duy nhất cặp giá trị ( q, r) với q ∈ Z và r ∈ N thỏa mãn: a = bq + r , 0 ≤ r < b . ở đây r gọi là số dư. Có một số ký hiệu cho số dư như sau: r = Rb (a ) , r = a mod b . Một số tính chất đơn giản của về số dư: 1. R−b ( a) = Rb ( a) , bởi vì a = qb + r  a = (− q )(−b) + r ⇔ 0 ≤ r < b  0 ≤ r < b 2. R b ( a + ib) = Rb ( a), ∀i ∈ Z , bởi vì a + ib = ( q + i )b + r . Nếu như r = 0 thì ta nói a chia hết cho b, ký hiệu là a . b Định lý 2.1: Đối với bất kỳ các số a, b, i ∈ Z : gcd(a, b) = gcd(a + ib, b) . Chứng minh: Nếu như ad và bd , thì a = q1d , b = q2 d . Lúc này a + ib = d ( q1 + iq2 ) ⇒ (a + ib) . Có nghĩa là: d gcd(a, b) ≤ gcd(a + ib, b) . (*) Tương tự giả sử rằng (a + ib) d và bd ⇒ a + ib = q1d , b = q2 d . Lúc này a = d (q1 − iq2 ) ⇒ ad . Có nghĩa là: gcd(a, b) ≥ gcd(a + ib, b) . (**) Từ (*) và (**) chúng ta có đẳng thức gcd(a,b)=gcd(a+ib,b). Định lý 2.2: Đối với bất kỳ a, b ∈ Z , b ≠ 0 , ta có: gcd(a, b) = gcd(b, Rb (a )) Chứng minh. Bởi vì a = bq + Rb (a) , gcd(a, b) = gcd(a − qb, b) = gcd( Rb ( a), b) . Từ định lý 2.2 ta có thuật toán Euclide. Đây là thuật toán giúp tìm UCLN của hai s ố nguyên dương a0 và a1 với a0 > a1. Thuật toán được miêu tả bằng dãy các phép chia như sau: a0 = a1q1 + a 2 , 0 < a2 < a1 a1 = a 2 q 2 + a3 , 0 < a3 < a2 7
… ak −2 = ak −1qk −1 + ak , 0 < ak < |ak-1| ak −1 = ak qk + 0 Dựa vào định lý 2.1, nhận được gcd(a0, a1) = gcd(a1, a2) = … = gcd(ak, 0) = ak. Ví dụ 2.17 (Thuật toán Euclide): Tìm gcd(814, 187). Giải: Khai triển dãy phép tính theo thuật toán Euclide: 814 = 4 * 187 + 66 187 = 2 * 66 + 55 66 = 1 * 55 + 11 55 = 5 * 11 + 0 Vậy gcd(814, 187) = 11. Thuật toán có thể viết như sau: Vào: Hai số nguyên dương a và b với a > b Ra: UCLN của a và b (1) While b # 0 do r ← a mod b, a ← b, b ← r (2) Return (a). 2.5.2 Thuật toán Euclide mở rộng Định nghĩa 2.3 (Phần tử nghịch đảo) Cho a ∈ Z n . Phần tử nghịch đảo của a là phần tử b ∈ Z n sao cho a * b ≡ b * a ≡ 1 mod n . Kí hiệu phần tử nghịch đảo của a là a-1. Định lý 2.2: Cho số nguyên a > 0 nguyên tố cùng nhau với n, thì luôn tồn tại phần tử nghịch đảo của a theo modulo n. Chứng minh: Hãy xét tập hợp {1, 2, ..., n − 1} . Nhân từng phần tử của tập hợp với a theo modulo n, nhận được tập hợp { ( a mod n), ( 2a mod n), ..., ((n − 1) a mod n)} , tập này sẽ bao gồm các số 1, 2,…, n – 1, có nghĩa đối với một số giá trị i nào đó sẽ thỏa mãn điều kiện ia mod n = 1 . Điều này dẫn đến một mâu thuẫn nếu như tồn tại hai giá trị h và k thỏa mãn điều kiện trên, nghĩa là ha mod n = ka mod n . Điều này dẫn đến h ≡ k mod n , vì gcd(a, n) = 1, suy ra h = k. Vậy ta tìm được i là phần tử nghịch đảo của a và i là duy nhất. Hệ quả 2.1: Nếu như p là số nguyên tố, thì bất kỳ số a, 0 < a < p, luôn tồn tại phần tử nghịch đảo theo modulo p. Chúng ta sẽ tìm hiểu về cách tìm phần tử nghịch đảo thông qua thuật toán Euclide mở rộng. Từ dãy chia của thuật toán Euclide a0 = a1q1 + a2 , 0 < a2 < a1 a1 = a2 q2 + a3 , 0 < a3 < a2 … 8
ak − 2 = ak −1qk −1 + ak , 0 < ak < ak-1 Ta dễ dàng rút ra dãy sau: ak = ak − 2 − qk −1ak −1 mà ak −1 = ak − 3 − qk − 2 ak − 2 nên suy ra ak = ak − 2 − qk −1 ( ak − 3 − qk − 2 ak − 2 ) , tương tự như thế, chúng ta rút ra ak − 2 , đến cuối cùng chúng ta có được biểu thức dạng sau: ak = a0 x + a1 y (2.1) Nếu như gcd(a0 , a1 ) = 1 , chúng ta có x là phần tử nghịch đảo của a0 theo modulo a1. Ví dụ 2.18 (Thuật toán Euclide mở rộng): 1. Tìm x và y trong biểu thức (2.1) với a0 = 814 và a1 = 187. Giải: Theo ví dụ 2.17 ta thu được dãy: 814 = 4 * 187 + 66 187 = 2 * 66 + 55 66 = 1 * 55 + 11 55 = 5 * 11 + 0 Suy ra: 11 = 66 – 1 * 55 = 66 – 1 * (187 – 2 * 66) = 3 * 66 – 1 * 187 = 3 * (814 – 4 * 187) –187 = 3 * 814 – 13 * 187. Vậy x = 3 và y = –13. 2. Tìm phần tử nghịch đảo của 17 theo modulo 74. Giải: Theo ví dụ trên ta có được 3 * 74 – 13 * 17 = 1. Nên ph ần t ử nghich đ ảo c ủa 17 là – 13, nhưng chỉ lấy số dương, nên phần tử nghịch đảo của 17 là 74 – 13 = 61. 2.6 Giải phương trình đồng dư trong nhóm thương Z n Có nhiều phương pháp giải phương trình đồng dư, nhưng chúng tôi mu ốn trình bày v ới bạn đọc một phương pháp khá hay dựa vào định lý sau: Định lý 2.3: Cho n > 1, gcd(a, n) = 1. Khi đó phương trình đồng dư ax ≡ b mod n có nghiệm duy nhất, và nghiệm đó là: x = baφ ( n ) − 1 mod n Chứng minh: Theo định lý Euler, ta có aφ ( n ) ≡ 1 mod n , suy ra a.baφ ( n ) −1 ≡ b mod n . Vậy nghiệm x = baφ ( n ) −1 mod n. Ví dụ 2.19: Giải phương trình 7x ≡ 3 mod 10. Chúng ta tính như sau: φ (10) = 4; x ≡ 3 * 74-1 mod 10 ≡ 1029 mod 10 ≡ 9. Thế nhưng cách này sẽ khó thực hiện nếu bậc của a đủ lớn. Bởi vậy chúng ta xem cách sau. Chúng ta đã biết định lý về phần tử nghịch đảo, nếu như gcd( a, n) = 1, thì luôn tồn tại phần tử nghịch đảo a-1. Nên từ phương trình ax ≡ b mod n , suy ra x ≡ ba −1 mod n. Ví dụ 2.20: Giải phương trình 7 x ≡ 3 mod 10 . Chúng ta thấy 7 −1 mod 10 = 3 . Nên nghiệm của phương trình x = 3 * 3 mod 10 = 9. 9
Định lý 2.4: Cho n >1, để phương trình ax ≡ b mod n có nghiệm thì điều kiện cần và đủ là gcd( a, n)|b. Chứng minh: Phương trình ax ≡ b mod n có thể viết dưới dạng phương trình tuyến tính sau ax + kn = b, (2.2), k là số nguyên. Ta chứng minh điều kiện cần. Giả sử hoàn thành điều kiện (2.2). B ời vì s ố gcd( a, n) là ước của vế trái, nên nó cũng phải là ước của vế phải. Chứng minh điều kiện đủ. Ứng dụng thuật toán Euclide đối với hai số a và n, có thể tính aλ + µn = gcd(a, n) (2.3) Bởi vì b/gcd(a, n) là số nguyên, nên nhân hai vế (2.3) cho b/gcd(a, n) và nhận được công thức (2.2), λb và ở đây x = mod n là một nghiệm. gcd(a, n) Dễ dàng kiểm tra phương trình đã có nghiệm dạng ni x+ mod n , i = 0, 1, 2, .., gcd(a, n) – 1. gcd(a, n) Tức là phương trình đã cho có gcd(a, n) nghiệm khác nhau, nghiệm này nhỏ hơn n. Ví dụ 2.21: 1. Giải phương trình 2 x ≡ 5 mod 10 Vì gcd(2, 10) = 2, vì 2 không là ước của 5 nên phương trình trên vô nghi ệm, n ếu không d ựa vào định lý cũng dễ dàng thấy rằng 2 x + 10y = 5 cũng vô nghiệm, vì vế trái là số chẳn, còn vế phải là s ố lẻ. Điều này là không thể. 2. Giải phương trình 6 x ≡ 18 mod 36 . Vì gcd(6, 36)|18. Nên phương trình này có 6 nghiệm: 3, 9, 15, 21, 27 và 33. 2.7 Định lý phần dư Trung Hoa. Định lý này nhằm giúp chúng ta giải hệ phương trình đồng dư thức, định lý phát biểu như sau: Giả sử cho các số n1, n2, …, nk là các số nguyên dương nguyên tố cùng nhau từng đôi m ột và c1, c2, …, ck là các số nguyên khi đó hệ phương trình đồng dư: x ≡ c1 mod n1 x ≡ c2 mod n2 … x ≡ ck mod nk Có nghiệm duy nhất theo modulo n và nghiệm đó là: k x = ∑ ci N i ( N i−1 mod ni ) mod n (2.4) i =1 M ở đây n = n1n2 ...nk và N i = . ni Chứng minh: 10
Nếu ta chứng minh hệ trên tương đương với một phương trình sau: x ≡ x0 mod n k −1 ở đây x0 = ∑ ci N i ( N i mod ni ) thì coi như chúng ta đã chứng minh được định lý trên. i =1 Ta có Ni chia hết cho ns, khi s ≠ i. Điều này dẫn đến x0 ≡ N i ( N i−1 mod ni )ci (mod ni ) , từ đó x0 ≡ ci (mod ni ) . Điều này có nghĩa hệ trên tương đương với hệ sau: x ≡ x0 mod n1 x ≡ x0 mod n2 ... x ≡ x0 mod nk Điều này tương đương với một phương trình sau x ≡ x0 mod n . Mà ta đã biết phương trình này có nghiệm duy nhất là x0 mod n. Ví dụ 2.22: Giải hệ phương trình đồng dư sau: x ≡ 5 mod 7 x ≡ 3 mod 11 x ≡ 10 mod 13 Giải: Ta có n = 7 * 11 * 13 = 1001; N1 = n/5 = 143; N2= n/11 = 91; N3 = n/13 = 77 Ta đi tính phần tử nghịch đảo của Ni theo modulo ni 143-1 mod 7 = 5; 91-1 mod 11 = 4 và 77-1 mod 13 = 12. Áp dụng công thức (2.4) ta có được nghiệm của hệ phương trình là: x = (5.143.5 + 3.91.4 + 10.77.12) mod 1001 = 894. 2.8 Bậc và căn nguyên thủy Định nghĩa 2.4 (Bậc) Cho gcd(a, n) = 1. Bậc a mod n là số nguyên dương nhỏ nhất γ thỏa mãn phương trình đồng dư thức: a γ ≡ 1 mod n , Kí hiệu bậc là Ordn(a). Ví dụ 2.23: Tìm bậc của 2 mod 5 Giải: 21 mod 5 = 2; 22 mod 5 = 4; 23 mod 5 = 3; 24 mod 5 = 1; vậy Ord5(2) = 4. Định nghĩa 2.5 (Căn nguyên thủy) Nếu như bậc của a modulo n bằng giá trị của hàm Euler của n. Thì a gọi là căn nguyên thủy của n. Ví dụ 2.24: Số 2 là căn nguyên thủy của 5. Vì Ord5(2) = φ (5) = 4. 2.9 Thặng dư bậc hai 11
Thặng dư bậc hai đóng vai trò rất quan trọng trong lý thuyết s ố. Ví d ụ, thu ật toán phân tích s ố nguyên ra thừa số. Ngoài ra thặng dư bậc hai cũng ứng dụng lớn trong mật mã cũng như trong các giao thức mã hóa. Định nghĩa 2.6 (Thặng dư bậc hai) Cho n là số nguyên và n > 1 . Số từ nhóm Z n gọi là thặng dư bậc hai theo modulo n, nếu trong * nhóm Z n tồn tại số x, thỏa mãn điều kiện x 2 ≡ a mod n , trường hợp ngược lại gọi là không thặng dư * bậc hai. Tập hợp các số thặng dư bậc hai theo modulo n ký hiệu QRn, tập hợp không thặng dư bậc hai ký hiệu là QNRn. Ví dụ 2.25: Tính QR11 là tập tất thặng dư bậc hai theo modulo 11. { } Chúng ta có QR11 = 12 , 2 2 , 32 , 4 2 , 52 , 6 2 , 7 2 , 82 , 9 2 ,10 2 mod 11 = {1, 3, 4, 5, 9} , trong ví dụ, chúng ta * tính QR11 bằng cách lựa chọn các phần tử của nhóm Z11 . { } Độc giả có thể kiểm tra lại rằng QR11 = 12 , 2 2 , 32 , 4 2 , 52 mod 11 = {1, 3, 4, 5, 9} . Điều này sẽ tìm hiểu qua định lý sau: Định lý 2.5: Cho p là số nguyên tố. Khi đó những khẳng định sau đây là đúng. 1. { QR p = x 2 (mod p) | 0 < x ≤ ( p − 1) / 2 .} * 2. Nhóm Z p chia ra hai tập con, QRp và QNRp có số phần tử bằng nhau. Chứng minh: Chúng ta chứng minh điều khẳng định thứ nhất. Rõ ràng rằng { } S = x 2 (mod p) | 0 < x ≤ ( p − 1) / 2 ⊆ QR p . Để chứng minh QR p = S , điều này sẽ đúng nếu như chúng ta chứng minh được rằng QR p ⊆ S . Chọn một số a bất kỳ a ∈ QR p . Tồn tại số x < p, thỏa mãn điều kiện x 2 ≡ a (mod p ) . Nếu như x ≤ ( p − 1) / 2 , thì a ∈ S . Giả sử rằng x > (p – 1)/2. Như thế y = p − x ≤ ( p − 1) / 2 và y 2 ≡ ( p − x) 2 ≡ p 2 − 2 px + x 2 ≡ x 2 ≡ a(mod p ) . Chứng tỏ, QR p ⊆ S . Ta đi chứng minh điều khẳng định thứ hai. Để chứng minh # QR p = ( p − 1) / 2 ,thì nó sẽ tương đương với việc chứng minh rằng đối với tất cả số x, thỏa mãn bất đẳng thức 0 < x < y ≤ ( p − 1) / 2 , thì thỏa mãn điều kiện x 2 ≠ y 2 (mod p ) . Giả sử điều này sai: x 2 − y 2 ≡ ( x + y )( x − y ) ≡ 0(mod p). Dẫn đến, p|(x + y) hoặc p|(x – y). Bởi vì x + y < p và p là số nguyên tố. Nên chỉ có thể là x = y. Điều này trái ngược với giả thuyết. Vậy # QR p = ( p − 1) / 2 . Mà # Z p = p − 1 , nên điều khẳng định thứ hai hoàn toàn * đúng. Thường xuất hiện bài toán, xác định xem số n có phải là thặng dư bậc hai theo modulo đã cho hay không. Vấn đề này gọi là bài toán về thặng dư bậc hai. Bài toán này đ ược giải quyết thông qua đ ịnh lý sau: 12
* Định lý 2.6 (Tiêu chuẩn Euler): Cho p là số nguyên tố. Đối với bất kỳ số x ∈ Z p , x là thặng dư bậc hai khi và chi khi thỏa mãn điều kiện sau: x ( p −1) / 2 ≡ 1 mod p . * Chứng minh: Đối với x ∈ QR p thì tồn tại số y ∈ Z p sao cho y 2 ≡ x(mod p ) . Theo định lý Fermat chúng ta có, x ( p −1) / 2 ≡ y p −1 ≡ 1 mod p . Ngược lại, nếu như x ( p −1) / 2 ≡ 1 mod p , như vậy x là nghiệm của đa thức y ( p −1) / 2 − 1 ≡ 0 mod p . Chúng ta đã biết Z p là một trường. Mà ta dễ dàng thấy rằng mỗi phần tử của trường là nghiệm của đa thức y p − y ≡ 0 mod p . Hay nói cách khác, mỗi phần tử khác không của trường là nghi ệm của đa thức: y p −1 − 1 ≡ ( y ( p −1) / 2 − 1)( y ( p −1) / 2 + 1) ≡ 0 mod p. Mà tất cả các nghiệm này phải khác nhau, bởi vì đa thức có bậc p – 1. Điều này dẫn đến (p – 1)/2 nghiệm của đa thức y ( p −1) / 2 − 1 ≡ 0 mod p phải khác nhau. Mà theo định lý về số phần tử của t ập thặng dư bậc hai ta có tập QRp có (p – 1)/2 phần tử, mỗi phần tử là nghiệm của phương trình * y ( p −1) / 2 − 1 ≡ 0 mod p . Và bất kỳ phần tử nào khác của nhóm Z p cũng phải thỏa mãn phương trình y ( p −1) / 2 + 1 ≡ 0 mod p . Vậy x ∈ QR p . Khi chứng minh định lý này chúng ta có thể điều khẳng định sau: N ếu như x ∈ Z p mà tiêu chuẩn Euler không thỏa mãn thì: x ( p −1) / 2 ≡ −1 mod p . 2.10 Các ký hiệu Legendre và Jacobi Ký hiệu Legendre là một công cụ hữu ích để xem xét liệu một số nguyên a có là thặng dư bậc hai theo modulo của một số nguyên tố p hay không? 2.10.1 Ký hiệu Legendre Định nghĩa 2.7 Nếu p là số nguyên tố lẻ và a là một số tự nhiên, thì ký hiệu Legendre a   được xác định như sau:  p   o 0 nếu a|p; o 1 nếu a ∈ QR p ; o −1 nếu a ∈ QNR p Các tính chất của ký hiệu Legendre: Cho p là một số nguyên tố lẻ và a, b ∈ Z. Khi đó ký hiệu Legendre có các tính chất sau: 13
 ab   a  b  1.   ≡     p      b  p  a b 2. Nếu a ≡ b mod p , thì   =    p  p     1 3.   = 1  p   p −1  1  1 khi p ≡ 1 mod 4 4.  −  = ( −1)  p 2    − 1 khi p ≡ 3 mod 4 p 2 −1 2  1 khi p ≡ 1, 7 mod 8 5.   = ( −1)  p 8 =   − 1 khi p ≡ 3, 5 mod 8  p +1 3  6   1 khi p ≡ 1,11 mod 12 6.   = ( −1)   p  =   − 1 khi p ≡ 5, 7 mod 12  p−2 5  5   1 khi p ≡ 1, 4 mod 5 7.   = ( −1)   p  =   − 1 khi p ≡ 2, 3 mod 5  p +1  7  6   1 khi p ≡ 1, 3, 9,19, 25, 27 mod 28 8.   = ( −1)   p  =   − 1 khi p ≡ 5,11,13,15,17, 23 mod 28 q −1 p −1  q   p 9. Nếu p và q là các số nguyên tố lẻ thì   =  ( −1)  p  q 2 2     Tiêu chuẩn Euler có thể viết lại như sau: Euler chứng minh rằng với mọi số nguyên tố p và số a, 1≤ a ≤ p , a   ≡ a ( p −1) / 2 mod p .  p   Ví dụ 2.26: Tính ký hiệu Legendre 14
 12345     331   3  5  823  =     331  331  331   3  5  161  =     331  331  331   3  5  7  23  =      331  331  331  331   331   331   331   331  = ( − 1)  ( − 1)  ( − 1)  ( − 1)    3   5   7   23   1  1  2  9  = −      3  5  7  23  2  1  1  2  3  = −      3  5  7  23   1  1  2  9  = −      3  5  7  23  = −(1)(1)(1)(1) = −1 2.10.2 Ký hiệu Jacobi Định nghĩa 2.8 Ký hiệu Jacobi là tổng quát hóa của ký hiệu Legendre cho số tự nhiên lẻ n. Giả sử α α n = p1 1 p2 2 ... pαk k là dạng phân tích tiêu chuẩn của n. Với số nguyên a bất kỳ, ký hiệu Jacobi là: α α αk a  a   a   a  1 2   =     ...  p  p  p  , trong đó tất cả các ký hiệu bên vế phải là ký hiệu Legendre. n  1  2   k Các tính chất của ký hiệu Jacobi: Cho n là các số tự nhiên lẻ và a, b ∈ Z. Kí hiệu Jacobi có các tính chất sau: 1. Nếu n là số nguyên tố thì ký hiệu Jacobi là ký hiệu Legendre a 2.   ∈ { 0, 1, − 1} n a 3.   = 0 khi gcd( a, n ) ≠ 1 n  ab   a  b  4.   ≡     n   n  n   a   a  b   a  5.   ≡    , điều này dẫn tới  2  là 0 hoặc 1.  mn   m  n  n  a b 6. Nếu a ≡ b mod n , thì   =   n n 1 7.   = 1 n 15
n −1  1  1 khi n ≡ 1 mod 4 8.  −  = ( −1) 2   n − 1 khi n ≡ 3 mod 4 n 2 −1  1 khi n ≡ 1, 7 mod 8 9.   = ( −1) 2   8 = n − 1 khi n ≡ 3, 5 mod 8 m n ( m −1)( n −1) 10.   = ( )( − 1) 4 n m Ví dụ 2.27:  37035   294   2  147   =  = 1.    331   331   331  331   147   331   37   147  = − = = =   331   147   147   37  2  36   2   9   37   1  =   =     =   =   =1  37   37   37   9   9  Vì 331 là số nguyên tố nên từ đó 37035 là thặng dư bậc hai mod 331. Có hai tính chất đúng với ký hiệu Legendre nhưng không thể mở rộng cho ký hiệu Jacobi. + Khác với ký hiệu Legendre, ký hiệu Jacobi không cho bi ết li ệu a có phải là một thặng dư bậc a a hai theo modulo n hay không. Sự thực nếu a ∈ QRn thì   = 1 . Tuy nhiên   = 1 thì không có nghĩa n n là a ∈ QRn. a ( p −1) / 2 + Tính chất thứ hai không thể mở rộng gắn với đồng dư thức Euler   ≡ a  p mod p với số   nguyên tố p và số nguyên a bất kỳ. Một cách tự nhiên đồng dư thức Euler m ở rộng t ừ ký hi ệu a ( n −1) / 2 Legendre sang ký hiệu Jacobi là   ≡ a mod n với hợp số lẻ dương n. Tuy nhiên đồng dư thức n này là sai với ít nhất một nửa của các a mod n khi n là hợp số. 2.10.3 Bổ đề Gauss Nếu như p và q là 2 số nguyên tố khác 2 thì: q −1 p −1  q   p   =  ( −1)  p  q 2 2     Hay nói cách khác, nếu như một trong 2 số p hoặc q có dạng 4n + 1, thì p bình phương theo modulo q khi và chỉ khi q là bình phương theo modulo p. Nếu như cả 2 số p và q có dạng 4n + 3, thì p bình phương theo modulo q khi và chỉ khi q không là bình phương theo modulo p. Chứng minh: 16
Cho rằng p ≠ q , xét trường Fp q −1 và cũng biết đây là cũng là nhóm cyclic. Theo định lý nhỏ Fermat thì p q −1 − 1 chia hết cho q, nên trong nhóm tồn tại phần tử có bậc là q. Chúng ta xem nó là w, lúc này w ≠ 1 và w q = 1 trong Fp q −1 . Bây giờ xác định tổng Gauss: x y = ∑  w x   x∈Z q  q  Chúng ta sẽ chứng minh 2 điều khẳng định sau: Điều khẳng định 1: Ta có đẳng thức sau q −1 2 y = (−1) 2 q . Chứng minh:  xz   x (u − x )  y 2 = ∑   w x + z = ∑ wu ∑  q    x,z   u∈Z q x∈Z q  q  Từ tổng cuối cùng x lấy giá trị trong tập Z q \ { 0} . Khi x ≠ 0 , có: q −1  x(u − x)   − x 2  1 − ux −1   1 − ux −1    =    q  q   = (−1) 2   q .   q       Từ đây q −1 (−1) 2 y 2 = ∑ Cu wu , u∈Z q  1 − ux −1  với Cu = ∑   q .  x∈Z q \{ 0}   Rõ ràng 1 C0 = ∑   = q −1.   x∈Z q \{ 0}  q  Nếu như u ≠ 0 , thì s = 1 − ux −1 nhận các giá trị trong tập Z q \ {1} , cho nên:  s 1 1 Cu = ∑   −   = −  , q q q s∈Z q       0 Bởi vì   = 0 , còn trong Z q \ { 0} số lượng phần tử là bình phương và số lượng phần tử không là  p   chính phương là như nhau. Từ đây: u u ∑ Cu w = ( q − 1) − ∑w = q. u∈Z q u∈Z q \{ 0} Vậy đã chứng minh được điều khẳng định 1, chúng ta chứng minh tiếp điều khẳng định tiếp theo. Khẳng định thứ 2: Ta có đẳng thức sau: 17
 p y p −1 =   . q   Chứng minh: Theo định thức Newton rút ra: x  zp −1  z  p −1   p y p = ∑   w xp = ∑  q     q w =  q  y =  q  y ,   x∈Z q   x∈Z q       Vậy điều khẳng định thứ hai đã được chứng minh. Kết hợp hai khẳng định trên chúng ta có đẳng thức: p −1  p  q −1  2 q −1 p −1 q   = y p−1  ( −1) 2 q  ( −1) 2 2  q =   =  p        2.11 Tính căn bậc hai theo modulo số nguyên 2.11.1 Tính căn bậc hai theo modulo là số nguyên tố Cho p là số nguyên tố, thì có 2 phương án tính như sau: Phương án 1: p ≡ 3, 7 mod 4 Trong trường hợp này số p + 1 chia hết cho 4. Cho a ∈ QR p . Ta định nghĩa x như sau: p +1 x=a 4 mod p , Bởi vì theo tiêu chuẩn Euler, a ( p −1) / 2 ≡ 1 mod p , nên x 2 ≡ a ( p+1) / 2 ≡ a ( p−1) / 2 a ≡ a mod p . Như vậy số x là căn bậc hai của a theo modulo p. Phương án 2: p ≡ 5 mod 8 Trường hợp này thì p + 3 chia hết cho 8. Bởi vì số (p – 1)/2 là số chẵn, số –1 thỏa mãn tiêu chuẩn Euler và là thặng dư bậc hai. Cho a ∈ QR p , ta định nghĩa x như sau: p+ 3 x=a 8 (mod p) , Từ điều kiện a ( p −1) / 2 ≡ 1 mod p ⇒ a ( p −1) / 4 ≡ ±1 mod p ⇒ x 2 ≡ a ( p+3) / 4 ≡ a ( p −1) / 4 a ≡ a mod p . Điều này chứng tỏ x được định nghĩa ở trên là căn bậc hai của a hoặc –a. Nếu như số dương thì thuật toán giải quyết xong, nếu số âm thì: − x 2 ≡ ( − 1x) 2 ≡ a(mod p) p+3 Dẫn đến x = − 1a 8 (mod p) là số cần tìm. Như vậy bài toán dẫn đến tìm − 1 mod p . Cho số b là số không thặng dư bậc hai theo modulo p. Như vậy theo tiêu chuẩn Euler ta có (b ( p −1) / 4 ) 2 ≡ b ( p −1) / 2 ≡ −1 mod p , như vậy có thể thay − 1 mod p bằng b (b−1) / 4 mod p . Bởi vì p 2 − 1 = ( p + 1)( p − 1) = (8k + 6)(8k + 4) = 8( 4k '+3)(2k "+1) . Mà ta đã biết 2 ∈ QNR p , nên có thể dùng 2( p −1) / 4 thay cho −1 18
Ngoài ra trong trường hợp tổng quát ta có thể tính như sau: Từ tiêu chu ẩn Euler, đ ồng d ư th ức x 2 ≡ a mod p giải được khi và chỉ khi a là thặng dư bậc hai theo modulo p, có nghĩa: p −1 a 2 ≡ 1 mod p . Giả sử rằng đẳng thức cuối cùng là đúng. Cũng có thể giả sử rằng, chúng ta bi ết m ột s ố s ố N là không thặng dư bậc hai theo modulo p, có nghĩa là biết N với tính chất: p −1 N 2 ≡ −1 mod p . Giả sử p − 1 = 2 s l , l là số lẻ. Chúng ta sẽ tìm dãy số bs−1 , ..., b0 sao cho: i (a l bi2 ) 2 ≡ 1 mod p . i −1 Có thể đặt bs −1 = 1 . Nếu như bi đã tìm được và i > 0 , khai căn, chúng ta nhận được rằng (a l bi2 ) 2 so sánh với 1 và –1 theo modulo p. Trong trường hợp thứ nhất thiết lập quan hệ bi −1 = bi , trong trường p −1 l +1 hợp thứ hai thì b = b N 2i +1 . Cuối cùng ta có a l b0 ≡ 1 mod p , từ đây (a 2 2 b )2 ≡ a mod p , có nghĩa i −1 i 0 là: l +1 x ≡ ±a 2 b 0 mod p . 2.11.2 Tính căn bậc hai theo modulo là hợp số Chúng ta biết rằng n = pq , p và q là số nguyên tố, nhóm Z n đồng cấu với không gian Z p x Z q . * * * Bởi vì đồng cấu đảm bảo tính chất lệnh số học nên x 2 ≡ y mod n , thỏa mãn khi và chỉ khi biểu thức trên đúng với modulo p và modulo q. Điều này có nghĩa nếu triển khai được n ra thừa số thì căn bậc hai theo modulo n có thể tính được theo thuật toán sau: Thuật toán 2.1: Đầu vào 2 số nguyên tố p, q thỏa mãn điều kiện n = pq, số nguyên y ∈ QR p . Đầu ra: Căn bậc hai của y theo modulo n xp y (mod p) xq y (mod q) return vp * xp + vq * xq mod n v p = q −1 (mod p ). p vq = p −1 (mod q).q 2.12 Thuật toán bình phương và nhân Thuật toán “bình phương và nhân” để tính z = x b mod n . Đây là thuật toán tính hàm mũ hiệu quả, mà trong các chương tiếp theo chúng ta sẽ áp dụng nó để tính giá trị hàm mũ đ ối v ới tr ường h ợp s ố mũ lớn. 19
Trong thuật toán này, ta coi rằng số mũ b được biểu thị ở dạng nhị phân như sau: l b = ∑ bi 2i i =0 Trong đó bi = 0 hoặc 1, 0 ≤ i ≤ l − 1 . Thuật toán tính: z = x b mod n gồm các bước như sau: 1. z=1 2. for i = l – 1 down to 0 do 3. z = z 2 mod n 4. if bi = 1 then z = z * x mod n Ví dụ : Tính 97263533 mod 11413 3533 = 211+210+28+27+26+23+22+1 = 110111001101 z=1 i bi z mod 11413 11 1 12 * 9726 = 9726 10 1 97262 * 9726 = 2659 9 0 26592 = 5634 8 1 56342 * 9726 = 9167 7 1 91672 * 9726 = 4958 6 1 49582 * 9726 = 7783 5 0 77832 = 6298 4 0 62982 = 4629 3 1 46292 * 9726 = 10185 2 1 101852 * 9726 = 105 1 0 1052 = 11025 0 1 110252 * 9726 = 5761 Vậy 97263533 mod 11413 = 5761 20