Greatest Common Divisor (Discrete Math)

Download as ppt, pdf, or txt
Download as ppt, pdf, or txt
You are on page 1of 29

Greatest Common Divisor

3 Gallon Jug

5 Gallon Jug

Lecture 8: Sep 29

This Lecture
In this lecture we will learn the Euclidean algorithm for computing greatest common divisor (GCD), which is one of the earliest important algorithms. Then we use the Euclidean algorithm to derive an important result in number theory, which is the basic in elementary number theory.

Quotient remainder theorem Greatest common divisor & Euclidean algorithm Linear combination and GCD, extended Euclidean algorithm Prime factorization and other applications

The Quotient-Remainder Theorem


For b > 0 and any a, there are unique numbers q ::= quotient(a,b), r ::= remainder(a,b), such that a = qb + r and 0 e r < b.

We also say

q = a div b

and

r = a mod b.

When b=2, this says that for any a, there is a unique q such that a=2q or a=2q+1.

When b=3, this says that for any a, there is a unique q such that a=3q or a=3q+1 or a=3q+2.

The Quotient-Remainder Theorem


For b > 0 and any a, there are unique numbers q ::= quotient(a,b), r ::= remainder(a,b), such that a = qb + r and 0 e r < b.

Given any b, we can divide the integers into many blocks of b numbers. For any a, there is a unique position for a in this line. q = the block where a is in r = the offset in this block

-b

2b

kb

a (k+1)b

Clearly, given a and b, q and r are uniquely defined.

This Lecture

Quotient remainder theorem Greatest common divisor & Euclidean algorithm Linear combination and GCD, extended Euclidean algorithm Prime factorization and other applications

Common Divisors
c is a common divisor of a and b means c|a and c|b. gcd(a,b) ::= the greatest common divisor of a and b.

Say a=8, b=10, then 1,2 are common divisors, and gcd(8,10)=2. Say a=10, b=30, then 1,2,5,10 are common divisors, and gcd(10,30)=10. Say a=3, b=11, then the only common divisor is 1, and gcd(3,11)=1.

Claim. If p is prime, and p does not divide a, then gcd(p,a) = 1.

Greatest Common Divisors


Given a and b, how to compute gcd(a,b)? Can try every number, but can we do it more efficiently?

Lets say a>b. 1. If a=kb, then gcd(a,b)=b, and we are done. 2. Otherwise, by the Quotient-Remainder Theorem, a = qb + r for r>0.

Greatest Common Divisors


Lets say a>b. 1. If a=kb, then gcd(a,b)=b, and we are done. 2. Otherwise, by the Quotient-Remainder Theorem, a = qb + r for r>0.

a=12, b=8 => 12 = 8 + 4 a=21, b=9 => 21 = 2x9 + 3 a=99, b=27 => 99 = 3x27 + 18

gcd(12,8) = 4 gcd(21,9) = 3 gcd(99,27) = 9

gcd(8,4) = 4 gcd(9,3) = 3 gcd(27,18) = 9

Euclid: gcd(a,b) = gcd(b,r)!

Euclids GCD Algorithm


a = qb + r Euclid: gcd(a,b) = gcd(b,r)

gcd(a,b) if b = 0, then answer = a. else write a = qb + r answer = gcd(b,r)

Example 1
gcd(a,b) if b = 0, then answer = a. else write a = qb + r answer = gcd(b,r)

GCD(102, 70) = GCD(70, 32) = GCD(32, 6) = GCD(6, 2) = GCD(2, 0)

102 = 70 + 32 70 = 2x32 + 6 32 = 5x6 + 2 6 = 3x2 + 0 Return value: 2.

Example 2
gcd(a,b) if b = 0, then answer = a. else write a = qb + r answer = gcd(b,r)

GCD(252, 189) = GCD(189, 63) = GCD(63, 0)

252 = 1x189 + 63 189 = 3x63 + 0

Return value: 63.

Example 3
gcd(a,b) if b = 0, then answer = a. else write a = qb + r answer = gcd(b,r)

GCD(662, 414) = GCD(414, 248) = GCD(248, 166) = GCD(166, 82) = GCD(82, 2) = GCD(2, 0)

662 = 1x414 + 248 414 = 1x248 + 166 248 = 1x166 + 82 166 = 2x82 + 2 82 = 41x2 + 0 Return value: 2.

Correctness of Euclids GCD Algorithm


a = qb + r Euclid: gcd(a,b) = gcd(b,r)

When r = 0: Then gcd(b, r) = gcd(b, 0) = b since every number divides 0. But a = qb so gcd(a, b) = b = gcd(b, r), and we are done.

Correctness of Euclids GCD Algorithm


a = qb + r When r > 0: Let d be a common divisor of b, r b = k1d and r = k2d for some k1, k2. a = qb + r = qk1d + k2d = (qk1 + k2)d Let d be a common divisor of a, b a = k3d and b = k1d for some k1, k3. r = a qb = k3d qk1d = (k3 qk1)d => d is a divisor of r => d is a divisor of a Euclid: gcd(a,b) = gcd(b,r)

So d is a common factor of a, b iff d is a common factor of b, r d = gcd(a, b) iff d = gcd(b, r)

How fast is Euclids GCD Algorithm?


Naive algorithm: try every number, Then the running time is about 2b iterations.

Euclids algorithm: In two iterations, then b is decreased by half. (why?) Then the running time is about 2log2(b) iterations.

Exponentially faster!!

This Lecture

Quotient remainder theorem Greatest common divisor & Euclidean algorithm Linear combination and GCD, extended Euclidean algorithm Prime factorization and other applications

Linear Combination vs Common Divisor


Greatest common divisor d is a common divisor of a and b if d|a and d|b gcd(a,b) = greatest common divisor of a and b

Smallest positive integer linear combination d is an integer linear combination of a and b if d=sa+tb for integers s,t. spc(a,b) = smallest positive integer linear combination of a and b

Theorem: gcd(a,b) = spc(a,b)

Linear Combination vs Common Divisor Theorem: gcd(a,b) = spc(a,b)


For example, the greatest common divisor of 52 and 44 is 4. And 4 is a linear combination of 52 and 44: 6 52 + (7) 44 = 4 Furthermore, no linear combination of 52 and 44 is equal to a smaller positive integer. To prove the theorem, we will prove: gcd(a,b) <= spc(a,b) spc(a,b) <= gcd(a,b) gcd(a,b) | spc(a,b) Write gcd as a positive integer linear combination

GCD <= SPC


3. If d | a and d | b, then d | sa + tb for all s and t. Proof of (3) d | a => a = dk1 d | b => b = dk2 sa + tb = sdk1 + tdk2 = d(sk1 + tk2) => d|(sa+tb)

Let d = gcd(a,b). By definition, d | a and d | b. GCD | SPC Let f = spc(a,b) = sa+tb

By (3), d | f. This implies d <= f. That is gcd(a,b) <= spc(a,b).

Extended GCD Algorithm


How can we write gcd(a,b) as an integer linear combination? This can be done by extending the Euclideans algorithm.

Example: a = 259, b=70 259 = 370 + 49 70 = 149 + 21 49 = a 3b 21 = 70 - 49 21 = b (a-3b) = -a+4b 49 = 221 + 7 7 = 49 - 221 7 = (a-3b) 2(-a+4b) = 3a 11b 21 = 73 + 0 done, gcd = 7

Extended GCD Algorithm


Example: a = 899, b=493 899 = 1 493 + 406 493 = 1 406 + 87 so 406 = a - b so 87 = 493 406 = b (a-b) = -a + 2b 406 = 4 87 + 58 so 58 = 406 - 4 87 = (a-b) 4(-a+2b) = 5a - 9b 87 = 1 58 + 29 so 29 = 87 1 58 = (-a+2b) - (5a-9b) = -6a + 11b 58 = 2 29 + 0 done, gcd = 29

This Lecture

Quotient remainder theorem Greatest common divisor & Euclidean algorithm Linear combination and GCD, extended Euclidean algorithm Prime factorization and other applications

Application of the Theorem Theorem: gcd(a,b) = spc(a,b)


Why is this theorem useful?

(1) we can now write down gcd(a,b) as some concrete equation, (i.e. gcd(a,b) = sa+tb for some integers s and t), and this allows us to reason about gcd(a,b) much easier. (2) If we can find integers s and t so that sa+tb=c, then we can conclude that gcd(a,b) <= c. In particular, if c=1, then we can conclude that gcd(a,b)=1.

Prime Divisibility Theorem: gcd(a,b) = spc(a,b)


Lemma: p prime and p|a b implies p|a or p|b. pf: say p does not divide a. so gcd(p,a)=1. So by the Theorem, there exist s and t such that sa + tp =1

(sa)b + (tp)b = b Hence p|b p|ab p|p p|ai for some i.

Cor : If p is prime, and p| a1 a2 am then Proof: by induction and the Lemma.

Fundamental Theorem of Arithmetic

Every integer, n>1, has a unique factorization into primes: p0 p1 p0 p1 pk pk = n

Example: 61394323221 = 3337111137373753

Unique Factorization
Theorem: There is a unique factorization.

proof: suppose, by contradiction, that there are numbers with two different factorization. By the well-ordering principle, we choose the smallest such n >1: n = p1 p2 pk = q1 q2 qm Since n is smallest, we must have that pi { qj all i,j (Otherwise, we can obtain a smaller counterexample.) Since p1|n = q1 q2 qm, so by Cor., p1|qi for some i. contradiction!

Since both p1 = qi are prime numbers, we must have p1 = qi.

Application of the Theorem Theorem: gcd(a,b) = spc(a,b)


Claim. If gcd(a,b)=1 and gcd(a,c)=1, then gcd(a,bc)=1.

By the Theorem, there exist s,t,u,v such that sa + tb = 1 ua + vc = 1 Multiplying, we have (sa + tb)(ua + vc) = 1 saua + savc + tbua + tbvc = 1 (sau + svc + tbu)a + (tv)bc = 1 By the Theorem, since spc(a,bc)=1, we have gcd(a,bc)=1

Die Hard (Optional)


Use a 3 gallon jug and a 5 gallon jug to fill in exactly 4 gallon of water.

3 Gallon Jug

5 Gallon Jug

Theorem: gcd(a,b) = spc(a,b)


Using this theorem, we can completely settle this problem: for any two jug sizes and for any target, we can either find a solution or show that no solutions exist (details omitted).

Quick Summary

Make sure you understand the Euclidean algorithm and the extended Euclidean algorithm. Also make sure to understand the relation between GCD and SPC. It is the basic of all the elementary number theory well see.

You might also like