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 tố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 tố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ố ngun 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,
P

ở đây tất cả các ap ≥ 0

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ố ngun 12 có thể trình bày như {a2 = 2, a3 = 1}.
Số ngun 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 ln ln 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.
Trong trường hợp chung:

2


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ố ngun 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 ln có thể tìm được q và r, phù hợp với
quan hệ tính tố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 hố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 .

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 .

4


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 


1 −
...1 −
φ (n) = n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,

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 .

5


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
Ta có một cơng thức quan trọng sau: aφ (n ) + 1 ≡ 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ố ngun 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
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 .

6


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) .

d
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 + iq 2 ) ⇒ (a + ib) .
Có nghĩa là:
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 − iq 2 ) ⇒ 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 tốn Euclide. Đây là thuật tố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
…
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 tốn Euclide): Tìm gcd(814, 187).
Giải: Khai triển dãy phép tính theo thuật tốn Euclide:
814 = 4 * 187 + 66
187 = 2 * 66 + 55
66 = 1 * 55 + 11

7

55 = 5 * 11 + 0
Vậy gcd(814, 187) = 11.
Thuật tốn có thể viết như sau:
Vào: Hai số ngun 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ì ln 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

…
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 tố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.

8


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ì ln 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.
Đị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ử hồ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 tốn Euclide đối với hai số a và n, có thể tính
aλ + µn = gcd(a, n) (2.3)

9


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),
và ở đây x =

λb
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
x+

ni
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 2x + 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

ở đây n = n1n2 ...nk và N i =

M
.
ni

Chứng minh:
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

10


...
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
Thặng dư bậc hai đóng vai trị rất quan trọng trong lý thuyết số. Ví dụ, thuật tốn phân tích số
ngun 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.

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 , 22 , 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 hồn tồn đúng.

Thường xuất hiện bài tố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 tốn này được giải quyết thơng qua định lý
sau:
*
Đị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:

12


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:
 ab   a  b 
1.   ≡   
 p
 
   b  p 
a b
2. Nếu a ≡ b mod p , thì   =  

 p  p
   
1
3.   = 1
 p
 
 1
4.  −  = ( −1)
 p


2
5.   = ( −1)
 p
 

p −1

2

p 2 −1
8

 p +1
6 



3
6.   = ( −1) 

 p
 

1 khi p ≡ 1 mod 4

− 1 khi p ≡ 3 mod 4

 1 khi p ≡ 1, 7 mod 8
=
− 1 khi p ≡ 3, 5 mod 8
 1 khi p ≡ 1,11 mod 12
=
− 1 khi p ≡ 5, 7 mod 12

13


 p−2
5 


 1 khi p ≡ 1, 4 mod 5
=
− 1 khi p ≡ 2, 3 mod 5

 p +1 
6 


 1 khi p ≡ 1, 3, 9,19, 25, 27 mod 28

=
− 1 khi p ≡ 5,11,13,15,17, 23 mod 28


5
7.   = ( −1) 
 p
 

7
8.   = ( −1) 
 p
 

 q   p
9. Nếu p và q là các số nguyên tố lẻ thì   =  ( −1)
 p  q
   

q −1 p −1
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
 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à:
α

α

1

2
 a
a  a   a 
=     ...
     
p
 n   p1   p2 
 k

αk






, trong đó tất cả các ký hiệu bên vế phải là ký hiệu Legendre.

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

14


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
 1
8.  −  = ( −1)
 n
2
9.   = ( −1)
 
n

n −1

2

1 khi n ≡ 1 mod 4


− 1 khi n ≡ 3 mod 4

n 2 −1
8

 1 khi n ≡ 1, 7 mod 8
=
− 1 khi n ≡ 3, 5 mod 8

( m −1)( n −1)
n
m
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ố ngun 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 hai
a
a
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 là a ∈
n
n
QRn.
a
( p −1) / 2
mod p với số
+ Tính chất thứ hai khơng thể mở rộng gắn với đồng dư thức Euler   ≡ a
 p
 
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 Legendre
a
( n −1) / 2
mod n với hợp số lẻ dương n. Tuy nhiên đồng dư thức này là sai
sang ký hiệu Jacobi là   ≡ a
n

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ì:

15


 q   p
  =  ( −1)
 p  q
   

q −1 p −1
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:
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 F p 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
y 2 = (−1)

q −1

2 q.

Chứng minh:
 xz 
 x (u − x ) 

y 2 = ∑   w x + z = ∑ wu ∑ 
q


q
x,z 
u∈Z q
x∈Z 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
 1 − ux −1 
 x(u − x)   − x 2  1 − ux −1 


 = (−1) 2 


=


 q .
q

q  q 

 





Từ đây
(−1)
với Cu =

q −1
2

y 2 = ∑ Cu wu ,
u∈Z q

 1 − ux −1 

∑ 
 q .
x∈Z q \{ 0} 


Rõ ràng
C0 =

1
∑   = 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    
 

16


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
∑ Cu w = ( q − 1) −

u∈Z q

u
∑w = 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:
 p
y p −1 =   .
q
 
Chứng minh: Theo định thức Newton rút ra:
 zp −1  z  p −1 
x
 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
  = y p−1
q
 

q −1


=  ( −1) 2 q 





p −1
2

= ( −1)

q −1 p −1

2 2 

q
 p

 

2.11Tí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ố ngun 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:
x=a

p +1
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:
x=a

p+ 3
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 tốn
giải quyết xong, nếu số âm thì:
− x 2 ≡ ( − 1x) 2 ≡ a(mod p)
Dẫn đến x = − 1a

p+3
8 (mod

p) là số cần tìm. Như vậy bài tốn dẫn đến tìm

− 1 mod p .

17


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

Ngồi ra trong trường hợp tổng qt 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:
a

p −1
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:
N

p −1
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 hợp
p −1

2
thứ hai thì b = b N 2i +1 . Cuối cùng ta có a l b0 ≡ 1 mod p , từ đây (a
i −1
i

x ≡ ±a

l +1
2 b
0

l +1
2 b )2
0

≡ a mod p , có nghĩa là:

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ố ngun 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 tốn sau:
Thuật tố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

18


2.12 Thuật tốn bình phương và nhân
Thuật tốn “bình phương và nhân” để tính z = x b mod n . Đây là thuật tố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.
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 tố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
11
10
9
8
7
6
5
4
3
2
1
0

bi
1
1

0
1
1
1
0
0
1
1
0
1

z mod 11413
12 * 9726 = 9726
97262 * 9726 = 2659
26592 = 5634
56342 * 9726 = 9167
91672 * 9726 = 4958
49582 * 9726 = 7783
77832 = 6298
62982 = 4629
46292 * 9726 = 10185
101852 * 9726 = 105
1052 = 11025
110252 * 9726 = 5761

Vậy 97263533 mod 11413 = 5761

19