1 / 21

Continued Fractions, Euclidean Algorithm and Lehmer’s Algorithm

Continued Fractions, Euclidean Algorithm and Lehmer’s Algorithm. Applied Symbolic Computation CS 567 Jeremy Johnson. TexPoint fonts used in EMF. Read the TexPoint manual before you delete this box.: A A A A A A A. Outline. Fast Fibonacci Continued Fractions Lehmer’s Algorithm

piper
Download Presentation

Continued Fractions, Euclidean Algorithm and Lehmer’s Algorithm

An Image/Link below is provided (as is) to download presentation Download Policy: Content on the Website is provided to you AS IS for your information and personal use and may not be sold / licensed / shared on other websites without getting consent from its author. Content is provided to you AS IS for your information and personal use only. Download presentation by click this link. While downloading, if for some reason you are not able to download a presentation, the publisher may have deleted the file from their server. During download, if you can't get a presentation, the file might be deleted by the publisher.

E N D

Presentation Transcript


  1. Continued Fractions, Euclidean Algorithm and Lehmer’s Algorithm Applied Symbolic Computation CS 567 Jeremy Johnson TexPoint fonts used in EMF. Read the TexPoint manual before you delete this box.: AAAAAAA

  2. Outline • Fast Fibonacci • Continued Fractions • Lehmer’s Algorithm • Analysis of the Euclidean Algorithm (bit complexity) • Assignment 1

  3. Euclidean Algorithm g = gcd(a,b) a1 = a; a2 = b; while (a2  0) a3 = a1 mod a2; a1 = a2; a2 = a3; } return a1;

  4. Remainder Sequence a1 = a, a2 = b a1 = q3 a2 + a3, 0  a3 < a2    ai = qi ai+1 + ai+2, 0  ai+2 < ai+1    al= ql al+1 gcd(a,b) = al+1

  5. Extended Euclidean Algorithm g = gcd(a,b,*x,*y) a1 = a; a2 = b; x1 = 1; x2 = 0; y1 = 0; y2 = 1; while (a2  0) a3 = a1 mod a2; q = floor(a1/a2); x3 = x1 – q*x2; y3 = y1 – q*y2; a1 = a2; a2 = a3; x1 = x2; x2 = x3; y1 = y2; y2 = y3; } return a1;

  6. Lehmer’s Algorithm u = 27182818, v = 10000000 u’ v’ q’ u’’ v’’ q’’ • 1001 2 2719 1000 2 1001 716 1 1000 719 1 716 285 2 719 281 2 285 146 1 281 157 1 146 139 1 157 124 1 139 7 19 124 33 3 u’/v’ < u/v < u’’/v’’

  7. Maximum Number of Divisions Theorem. Let a  b  0 and n = number of divisions required by the Euclidean algorithm to compute gcd(a,b). Then n < 2lg(a).

  8. Maximum Number of Divisions Theorem. The smallest pair of integers that require n divisions to compute their gcd is Fn+2 and Fn+1. Theorem. Let a  b  0 and n = number of divisions required by the Euclidean algorithm to compute gcd(a,b). Then n < 1.44lg(a).

  9. Average Number of Divisions Theorem. Let a  b  0 and n = average number of divisions required by the Euclidean algorithm to compute gcd(a,b). Then n  12ln(2)2/2lg(a)  0.584 lg(a). Theorem [Dixon] D(a,b)  ½ ln(a) for almost all pairs u  a  b  1 as u  

  10. Dominance and Codominance Definition. Let f, g be real valued functions on a common set S. • [Dominance] • [Codominance] • [Strict Dominance]

  11. Basic Properties Theorem. Let f, f1, f2, g, g1, and g2 be nonnegative real-valued functions on S and c>0.

  12. Integer Length Definition. A = i=0..maii, L(A) = m

  13. Basic Arithmetic Computing Times Theorem. Let A, M, D be the classical algorithms for addition, multiplication and division.

  14. Maximum Computing Time Theorem.

  15. Average Computing Time Theorem.

  16. Probability of Relative Primality p/d2 = 1  p = 1/(2) (z) = 1/nz (2) = 2/6

  17. Formal Proof Let qn be the number of 1 a,b n such that gcd(a,b) = 1. Then limnqn/n2 = 6/2

  18. Mobius Function • (1) = 1 • (p1    pt) = -1t • (n) = 0 if p2|n (ab) = (a)(b) if gcd(a,b) = 1.

  19. Mobius Inversion d|n (d) = 0 (n (n)ns)(n 1/ns) = 1

  20. Formal Proof qn = n(k)n/k2 limnqn/n2 = n (n)/n2= 1/(n1/n2) = 6/2

  21. Assignment 1 • Empirically investigate distribution of quotients in Euclidean algorithm • Implement and analyze classical division algorithm • Implement and analyze Lehmer’s algorithm • Study and summarize gcd algorithms in GMP

More Related