COMPUTING π(x): THE MEISSEL, LEHMER, LAGARIAS, MILLER, ODLYZKO METHOD
|
|
- Sophie Lane
- 6 years ago
- Views:
Transcription
1 MATHEMATICS OF COMPUTATION Volume 65, Number 213 Januar 1996, Pages COMPUTING π(): THE MEISSEL, LEHMER, LAGARIAS, MILLER, ODLYZKO METHOD M. DELEGLISE AND J. RIVAT Abstract. Let π() denote the number of primes. Our aim in this paper is to present some refinements of a combinatorial method for computing single values of π(), initiated b the German astronomer Meissel in 1870, etended and simplified b Lehmer in 1959, and improved in 1985 b Lagarias, Miller and Odlzko. We show that it is possible to compute π() ino( 2/3 log 2 )time and O( 1/3 log 3 log log ) space. The algorithm has been implemented and used to compute π(10 18 ). 1. Introduction One of the oldest problems in mathematics is to compute π(), the eact number of primes. The most obvious method for computing π() is to find and count all primes p, for instance b the sieve of Eratosthenes. According to the Prime Number Theorem (1) π() log,. Therefore, such a method cannot compute π() with less than about log operations. Despite its time compleit, the sieve of Eratosthenes has been for a ver long time the practical wa to compute π(). In the second half of the 19th centur, the astronomer Meissel discovered a practicable combinatorial method that is faster than finding all primes. He used his algorithm to compute b hand π(10 8 )and π(10 9 ) (which turned out to be too small b 56) [4, 5, 6, 7]. In 1959, Lehmer etended and simplified Meissel s method. He used an IBM 701 computer to obtain the value of π(10 10 ) (his value was shown [1] to be too large b 1). In 1985, Lagarias, Miller and Odlzko [2] adapted the Meissel-Lehmer method and proved that it is possible to compute π() with O( 2/3 log ) operations using O( 1/3 log 2 log log ) space. The used their algorithm to compute several values of π() upto= The also corrected the value of π(10 13 ) given in [1], which was too small b 941. In 1987, Lagarias and Odlzko [3] described a completel different method, based on numerical integration of certain integral transforms of the Riemann ζ-function, for computing π(), using O( 1/2+ε )timeando( 1/4+ε ) space for each ε>0. Received b the editor Januar 12, 1994 and, in revised form, December 1, Mathematics Subject Classification. Primar 11N05, 11Y c 1996 American Mathematical Societ
2 236 M. DELEGLISE AND J. RIVAT Despite its asmptotic superiorit, this algorithm has never been implemented. Its authors noted [2] that the implied constants are probabl large, and therefore that it would not be competitive with their version of the Meissel-Lehmer method for In this paper we describe a modified form of the algorithm presented in [2] which computes π() usingo( 2/3 log 2 )timeando(1/3 log 3 log log ) space. 2. Outline of the method For clarit we will describe the whole method given in [2], in order to introduce the quantities needed for the analsis. For the convenience of the reader we adopt the notations used in [2], and follow as long as possible the same approach. In particular, 3, 4, 5 are close to [2]. The idea that man special leaves (see below, 6) could be computed at the same time, saving much computation (and a log factor in the compleit) was also present in [2]. We develop this idea further, and show that it is possible to compute more special leaves at the same time, saving another log factor in the compleit (see 6.2 and below). 3. The Meissel-Lehmer method Let p 1,p 2,p 3,... denote the primes 2, 3, 5,... numbered in increasing order. Let (, a) denote the partial sieve function, which counts numbers with all prime factors greater than p a : (2) (, a) =#{n ; p n p>p a } and let (3) P k (, a) =#{n ; n=q 1 q 2 q k, q 1,...,q k >p a }, which counts numbers with eactl k prime factors, all larger than p a.weset P 0 (, a) =1. If we sort the numbers b the number of their prime factors greater than p a, we obtain the following identit: (, a) =P 0 (, a)+p 1 (, a)+ +P k (, a)+, where the sum on the right has onl finitel man nonzero terms, because P k (, a) = 0forp k a >. Let denote an integer such that 1/3 1/2,andleta=π(). From P 1 (, a) =π() aand P k (, a) =0fork 3 we deduce (4) π() =(, a)+a 1 P 2 (, a). Hence, for computing π() it remains to compute (, a) andp 2 (, a). 4. Computing P 2 (, a) B (3) we have to count all pairs (p, q) of prime numbers such that <p q and. We first remark that p [ +1, ]. Furthermore, for each p, wehaveq [p, /p]. Thus, P 2 (, a) = ( ) (5) π π(p)+1. p <p
3 COMPUTING π() 237 When p [ +1, ], we have p [1, ]. Hence, P 2(, a) can be computed b completel sieving the interval [1, ] and then adding up π( p ) π(p) + 1 for all primes p [ +1, ]. In order to reduce the space compleit of the above method, we can work with blocks of length L. ForL=, we can compute P 2 (, a) in O( log log ) timeando() space. 5. The sieve machiner for computing (, a) For b a the set of all integers whose prime factors are >p b 1 is composed of two classes: 1. those that are multiples of p b, 2. those not divisible b p b. The first class has ( p b,b 1) elements, while the second has (, b) elements. Hence we conclude: Lemma 5.1. The function satisfies the following identities: (6) (7) (u, 0) = [u], (, b) =(, b 1),b 1. p b A straightforward method for computing (, a) can be deduced from this lemma: it suffices to appl repeatedl the recurrence (7) until we get terms of the form (u, 0), which are eas to compute using (6). One ma think of this process as creating a rooted binar tree starting with the root node (, a);seefig.1.using this method, we obtain the following formula: (, a) = 1 n P + (n) [ µ(n), n] where µ(n) denotes the Möbius function and P + (n) denotes the greatest prime factor of n. (, a) (, a 1) ( p, a 1) a (, a 2) ( p, a 2) a 1 ( p, a 2) (, a 2) a p a p a 1 Figure 1. A binar tree for computing (, a): the sum of the terminal nodes is (, a)
4 238 M. DELEGLISE AND J. RIVAT Unfortunatel, this sum has too man terms for our purpose: as 1/3,ifwe onl count the n s which are the product of three primes, we get at least about such terms. log 3 In order to limit the growth of the tree, we must replace the trivial truncation rule, Truncation Rule 1. Do not split a node µ(n)( n,b) if b =0 with the more powerful: Truncation Rule 2. Do not split a node µ(n)( n,b) if either of the following holds: 1. b =0and n, 2. n>. We are now able to define two clases of leaves: 1. ordinar leaves are those of the form µ(n)( n, 0) satisfing n, 2. special leaves are those of the form µ(n)( n,b 1) satisfing n>with n = mp b and m. We conclude that: Lemma 5.2. We have (8) (, a) =S 0 +S, where S 0 is the contribution of the ordinar leaves, S 0 = [ (9) µ(n), n] n and S is the contribution of the special leaves, S = ( ) (10) µ(n) n n,π(δ(n)) 1, δ(n) <n where δ(n) denotes the smallest prime factor of n. The computation of S 0 can be achieved in O( log log ) time, which is negligible. It remains to compute S. (11) We have S = p 6. Computing S δ(m)>p m <mp µ(m) mp,π(p) 1.
5 COMPUTING π() 239 We now write with S 1 = S 2 = S 3 = 1 3 <p δ(m)>p m <mp 1 4 <p 1 3 p 1 4 S = S 1 + S 2 + S 3 δ(m)>p m <mp δ(m)>p m <mp µ(m) mp,π(p) 1, µ(m) µ(m) mp,π(p) 1, mp,π(p) 1. First observe that the m s involved in S 1 and S 2 are all prime: otherwise, since δ(m) >p> 1/4,wewouldhavem>p 2 >, a contradiction with m. Moreover, the condition pm is true when mp > 1/2. Hence we have S 1 =,π(p) 1, S 2 = 1 3 <p 1 4<p 1 3,π(p) Computing S 1. Since <1/3 <p, we have,π(p) 1 =1. Hence all terms involved in S 1 are equal to 1. So we have to count all pairs (p, q) such that Thus, 1/3 <. S 1 = (π() π(1/3 ))(π() π( 1/3 ) 1). 2 This takes constant time to compute S Computing S 2. We have S 2 = 1/4 <p 1/3,π(p) 1.
6 240 M. DELEGLISE AND J. RIVAT We split S 2 in two parts, depending on q>/p 2 or q /p 2 : S 2 = U + V with and U = V = 1/4 <p 1/3 1/4 <p 1/3 q> p 2 q p 2,π(p) 1,π(p) Computing U. The condition q>/p 2 implies p 2 >/q / and p> /. Thus, U= <p 1/3 q> p 2,π(p) 1. From /p 2 <qwe deduce / < p and (/, π(p) 1) = 1. Each term in the sum U equals 1. Hence, U = # {q; p } 2 <q. <p 1/3 Thus, U = <p 1/3 ( ) π() π p 2. Since /p 2 <,thesumucan be calculated in O() operationsoncewehave tabulated π(t) for t Computing V. For each term involved in V we have p <1/2 <p 2. Hence,,π(p) 1 =1+π (π(p) 1) Thus, with =2 π(p)+π V = V 1 + V 2 ( ).
7 V 1 = V 2 = 1/4 <p 1/3 COMPUTING π() 241 p<q min( p 2,) (2 π(p)), π 1/4 <p 1/3 p<q min( p 2,). Computing V 1 can be achieved in O( 1/3 ) time once we have tabulated π(t) for t. In order to speed up the computation of V 2, we observe that for each p we can split the summation over q into sums over q on intervals where the function q π( ) is constant. Thus, we onl need the length of these intervals, and the set of values of q where q π( ) is changing. More precisel, we first split V 2 in two parts in order to simplif the condition q min(/p 2,): V 2 = We now write 1/4 <p π + <p 1/3 V 2 = W 1 + W 2 + W 3 + W 4 + W 5 p<q p 2 π. with W 1 = W 2 = W 3 = W 4 = W 5 = 1/4 <p 2 <p 2 <p 2 <p 1/3 <p 1/3 p<q p p <q π p<q p π, p <q p 2 π π,, π,. Computing W 1 and W 2. These two quantities need values of π(/) with< / < 1/2. The are computed simultaneousl with a sieve of the interval [1, ]. The sieving is done b blocks, and for each block we sum π(/) for the pairs (p, q) subjected to the conditions of the sums W 1 or W 2 andsuchthat/ lies in the block. Computing W 3. For each p we speed up the computation of the sum over q b computing in O(1) operations the sums of the π(/) for the values of q for which π(/) is constant. When we obtain a new value of q, we compute π(/) with
8 242 M. DELEGLISE AND J. RIVAT the table of values of π(t) fort. Then a table of all primes gives t such that π(t) <π(t+1)=π( ). We then deduce the net value of q for which π(/) is changed. Computing W 4. We simpl sum over (p, q). There would be no advantage to proceed as for W 3 since most of the values π(/) are distinct. Computing W 5. We proceed as for W Computing S 3 We sieve the interval [1, ] successivel b all primes less than 1/4. Assoonaswe have sieved b p k 1,wesumall µ(m)( mp k,k 1) for all squarefree m [/p, ] such that δ(m) > p k. This computation can be done b blocks, see [2]. The main idea is that we maintain a binar tree (as eplained in [2, pp. 545, 546]) in connection with the interval we are sieving, to keep track of the intermediate results after sieving b all primes up to a given prime. It is then possible to know the number of unsieved elements in the interval less than a given value, using onl O(log ) operations. 8. Time and space compleit The time and space significant computations are: 1. The computation of P 2 (, a), 2. The computation of W 1,W 2,W 3,W 4,W 5, 3. The computation of S Cost of computing P 2 (, ). We have alread seen that it costs O( log log ) time and O() space Cost of computing W 1,W 2,W 3,W 4,W 5. For W 1 and W 2 the sieve costs O( log log ) timeando() space, working b blocks of length. ThetimenecessartocomputethesumW 1 is about π 2 π() =O log 2. ThetimenecessartocomputethesumW 2 is about O <p 2 π ( ) 3/4 = O p 1/4 log 2. In W 3,foreachpwe have /p. Hence, π(/) takes at most π( /p) values. For each such value it costs constant time to determine the number of q s concerned. Hence, the time necessar to compute the sum W 3 is about O <p 2 π ( ) 3/4 = O p 1/4 log 2.
9 COMPUTING π() 243 ThetimenecessartocomputethesumW 4 is about ( ) 2/3 O π = O p log 2. <p 1/3 We proceed for W 5 as for W 3, and the time necessar to compute the sum W 5 is about ( ) 2/3 O π = O log 2. <p 1/ Cost of computing S 3. The sieve. Owing to the necessit of quickl retrieving the values (u, b), we have to maintain a data structure such that each access costs O(log ) instead of O(1) in a normal sieve. Hence the cost is O( log log log ). The sum. For each term in the sum we have to access the data structure mentioned above, doing O(log ) operations. It remains to count the terms in the sum. All leaves are of the form ±(/mp b,b 1) with m and b<π( 1/4 ). Hence, the number of these leaves is O(π( 1/4 )). The total cost of S 3 is O log log log + 1/ Total cost. We have described an algorithm taking O() space and ( O log log + ) log log log + 1/4 + 2/3 log 2 time. If we choose = 1/3 log 3 log log, the time compleit is O( 2/3 log 2 ). 9. Practical considerations We describe here some modifications which improve the time of computation without changing the asmptotic compleit. In the truncation rule 2, we ma replace b some z>. It is possible to prove that the time compleit for computing S 3 then becomes 1/4 O log log log + z log + z3/2. This also gives a good wa for checking the computations b changing the value of z. For clarit we chose to split the sum S at 1/4, but in fact we onl need to have p <p2. One can take advantage of this, but the asmptotic compleit remains the same. Precomputing the sieving b the first primes 2, 3, 5 saves some more time.
10 244 M. DELEGLISE AND J. RIVAT 10. Results The algorithm has been implemented in C + +. All the computations were done using a HP 730 workstation (SPEC92INT=47.8). The 64-bit integers were emulated b the long long tpe of GNU C/C++ Compiler. For comparison we tried our program for some specific values on a DEC Alpha 3000 Model 600 at 175 Mhz (which has 64-bit integers, SPEC92INT=114). The latter turned out to be more than three times faster. The difference could be greater because our program was optimized for a 32-bit computer, which is a drawback on a 64-bit computer. We confirmed all the values alread computed in [2]. Table 1 gives the new values compared with the corresponding values of Li() = 0 dt log t Table 1. Results and times of computation on HP-730 π() Li() π() R() π() Time (s)
11 COMPUTING π() 245 and R() = n=1 µ(n) n Li(1/n ). References 1. J. Bohman, On the number of primes less than a given limit, BIT12 (1972), MR 48 # J. C. Lagarias, V. S. Miller, and A. M. Odlzko, Computing π(): The Meissel-Lehmer method, Math. Comp. 44 (1985), MR 86h: J. C. Lagarias and A. M. Odlzko, Computing π(): An analtic method, J. Algorithms 8 (1987), MR 88k: E. D. F. Meissel, Über die Bestimmung der Primzahlenmenge innerhalb gegebener Grenzen, Math. Ann. 2 (1870), , Berechnung der Menge von Primzahlen, welche innerhalb der ersten hundert Millionen natürlicher Zahlen vorkommen, Math. Ann. 3 (1871), , Über Primzahlenmengen, Math. Ann. 21 (1883), , Berechnung der Menge von Primzahlen, welche innerhalb der ersten Milliarde natürlicher Zahlen vorkommen, Math. Ann. 25 (1885), Département de Mathématiques, Université Lon 1, 43 Blvd. du 11 Novembre 1918, Villeurbanne Cede, France address: deleglis@lmdi.univ-lon1.fr address: rivat@caissa.univ-lon1.fr
COMPUTING π(x): THE MEISSEL-LEHMER METHOD. J. C. Lagarias AT&T Bell Laboratories Murray Hill, New Jersey
COMPUTING π(x): THE MEISSEL-LEHMER METHOD J. C. Lagarias AT&T Bell Laboratories Murray Hill, New Jersey V. S. Miller IBM Watson Research Center Yorktown Heights, New York A. M. Odlyzko AT&T Bell Laboratories
More informationCOUNTING PRIMES IN RESIDUE CLASSES
MATHEMATICS OF COMPUTATION Volume 73, Number 247, Pages 1565 1575 S 0025-5718(04)01649-7 Article electronically published on February 25, 2004 COUNTING PRIMES IN RESIDUE CLASSES MARC DELÉGLISE, PIERRE
More informationExplicit estimates on the summatory functions of the Moebius function with coprimality restrictions
Eplicit estimates on the summator functions of the Moebius function with coprimalit restrictions Olivier Ramaré, CNRS, Laboratoire Paul Painlevé, Université Lille, 59 55 Villeneuve d Ascq, France Email:
More informationMA123, Chapter 8: Idea of the Integral (pp , Gootman)
MA13, Chapter 8: Idea of the Integral (pp. 155-187, Gootman) Chapter Goals: Understand the relationship between the area under a curve and the definite integral. Understand the relationship between velocit
More informationComputing π(x): the combinatorial method
REVISTA DO DETUA, VOL. 4, N 6, MARCH 2006 759 Computing π: the combinatorial method Tomás Oliveira e Silva Abstract This article presents a self-contained description of the method proposed by Deléglise
More informationFermat numbers and integers of the form a k + a l + p α
ACTA ARITHMETICA * (200*) Fermat numbers and integers of the form a k + a l + p α by Yong-Gao Chen (Nanjing), Rui Feng (Nanjing) and Nicolas Templier (Montpellier) 1. Introduction. In 1849, A. de Polignac
More informationMAT389 Fall 2016, Problem Set 2
MAT389 Fall 2016, Problem Set 2 Circles in the Riemann sphere Recall that the Riemann sphere is defined as the set Let P be the plane defined b Σ = { (a, b, c) R 3 a 2 + b 2 + c 2 = 1 } P = { (a, b, c)
More informationAppendix 2 Complex Analysis
Appendix Complex Analsis This section is intended to review several important theorems and definitions from complex analsis. Analtic function Let f be a complex variable function of a complex value, i.e.
More informationRANK AND PERIOD OF PRIMES IN THE FIBONACCI SEQUENCE. A TRICHOTOMY
RANK AND PERIOD OF PRIMES IN THE FIBONACCI SEQUENCE. A TRICHOTOMY Christian Ballot Université de Caen, Caen 14032, France e-mail: ballot@math.unicaen.edu Michele Elia Politecnico di Torino, Torino 10129,
More informationSection 3.1. ; X = (0, 1]. (i) f : R R R, f (x, y) = x y
Paul J. Bruillard MATH 0.970 Problem Set 6 An Introduction to Abstract Mathematics R. Bond and W. Keane Section 3.1: 3b,c,e,i, 4bd, 6, 9, 15, 16, 18c,e, 19a, 0, 1b Section 3.: 1f,i, e, 6, 1e,f,h, 13e,
More informationMethods for Advanced Mathematics (C3) Coursework Numerical Methods
Woodhouse College 0 Page Introduction... 3 Terminolog... 3 Activit... 4 Wh use numerical methods?... Change of sign... Activit... 6 Interval Bisection... 7 Decimal Search... 8 Coursework Requirements on
More informationFACULTY OF MATHEMATICAL STUDIES MATHEMATICS FOR PART I ENGINEERING. Lectures AB = BA = I,
FACULTY OF MATHEMATICAL STUDIES MATHEMATICS FOR PART I ENGINEERING Lectures MODULE 7 MATRICES II Inverse of a matri Sstems of linear equations Solution of sets of linear equations elimination methods 4
More information10.2 The Unit Circle: Cosine and Sine
0. The Unit Circle: Cosine and Sine 77 0. The Unit Circle: Cosine and Sine In Section 0.., we introduced circular motion and derived a formula which describes the linear velocit of an object moving on
More informationMAT 1275: Introduction to Mathematical Analysis. Graphs and Simplest Equations for Basic Trigonometric Functions. y=sin( x) Function
MAT 275: Introduction to Mathematical Analsis Dr. A. Rozenblum Graphs and Simplest Equations for Basic Trigonometric Functions We consider here three basic functions: sine, cosine and tangent. For them,
More informationLEGENDRE TYPE FORMULA FOR PRIMES GENERATED BY QUADRATIC POLYNOMIALS
Ann. Sci. Math. Québec 33 (2009), no 2, 115 123 LEGENDRE TYPE FORMULA FOR PRIMES GENERATED BY QUADRATIC POLYNOMIALS TAKASHI AGOH Deicate to Paulo Ribenboim on the occasion of his 80th birthay. RÉSUMÉ.
More informationON SUMS OF SEVEN CUBES
MATHEMATICS OF COMPUTATION Volume 68, Number 7, Pages 10 110 S 005-5718(99)01071-6 Article electronically published on February 11, 1999 ON SUMS OF SEVEN CUBES F. BERTAULT, O. RAMARÉ, AND P. ZIMMERMANN
More informationTHE NUMBER OF GOLDBACH REPRESENTATIONS OF AN INTEGER
PROCEEDINGS OF THE AMERICAN MATHEMATICAL SOCIETY Volume 140, Number 3, March 01, Pages 795 804 S 000-9939011)10957- Article electronicall published on Jul 0, 011 THE NUMBER OF GOLDBACH REPRESENTATIONS
More informationSums of Squares. Bianca Homberg and Minna Liu
Sums of Squares Bianca Homberg and Minna Liu June 24, 2010 Abstract For our exploration topic, we researched the sums of squares. Certain properties of numbers that can be written as the sum of two squares
More information4Cubic. polynomials UNCORRECTED PAGE PROOFS
4Cubic polnomials 4.1 Kick off with CAS 4. Polnomials 4.3 The remainder and factor theorems 4.4 Graphs of cubic polnomials 4.5 Equations of cubic polnomials 4.6 Cubic models and applications 4.7 Review
More informationCourse Requirements. Course Mechanics. Projects & Exams. Homework. Week 1. Introduction. Fast Multipole Methods: Fundamentals & Applications
Week 1. Introduction. Fast Multipole Methods: Fundamentals & Applications Ramani Duraiswami Nail A. Gumerov What are multipole methods and what is this course about. Problems from phsics, mathematics,
More informationJournal of Inequalities in Pure and Applied Mathematics
Journal of Inequalities in Pure and Applied Mathematics http://jipam.vu.edu.au/ Volume 2, Issue 2, Article 15, 21 ON SOME FUNDAMENTAL INTEGRAL INEQUALITIES AND THEIR DISCRETE ANALOGUES B.G. PACHPATTE DEPARTMENT
More informationProperties of Limits
33460_003qd //04 :3 PM Page 59 SECTION 3 Evaluating Limits Analticall 59 Section 3 Evaluating Limits Analticall Evaluate a it using properties of its Develop and use a strateg for finding its Evaluate
More informationChapter 4 Analytic Trigonometry
Analtic Trigonometr Chapter Analtic Trigonometr Inverse Trigonometric Functions The trigonometric functions act as an operator on the variable (angle, resulting in an output value Suppose this process
More informationEvaluation of some integrals of sums involving the Möbius function
International Journal of Computer Mathematics Vol. 84, o. 4, April 2007, 469 475 Evaluation of some integrals of sums involving the Möbius function TADEJ KOTIK* Faculty of Electrical Engineering, University
More informationThe Maze Generation Problem is NP-complete
The Mae Generation Problem is NP-complete Mario Alviano Department of Mathematics, Universit of Calabria, 87030 Rende (CS), Ital alviano@mat.unical.it Abstract. The Mae Generation problem has been presented
More information14.1 Systems of Linear Equations in Two Variables
86 Chapter 1 Sstems of Equations and Matrices 1.1 Sstems of Linear Equations in Two Variables Use the method of substitution to solve sstems of equations in two variables. Use the method of elimination
More informationChapter 6. Self-Adjusting Data Structures
Chapter 6 Self-Adjusting Data Structures Chapter 5 describes a data structure that is able to achieve an epected quer time that is proportional to the entrop of the quer distribution. The most interesting
More informationSection 1.2: A Catalog of Functions
Section 1.: A Catalog of Functions As we discussed in the last section, in the sciences, we often tr to find an equation which models some given phenomenon in the real world - for eample, temperature as
More informationarxiv: v2 [math.nt] 29 Jul 2017
Fibonacci and Lucas Numbers Associated with Brocard-Ramanujan Equation arxiv:1509.07898v2 [math.nt] 29 Jul 2017 Prapanpong Pongsriiam Department of Mathematics, Faculty of Science Silpakorn University
More informationLESSON 35: EIGENVALUES AND EIGENVECTORS APRIL 21, (1) We might also write v as v. Both notations refer to a vector.
LESSON 5: EIGENVALUES AND EIGENVECTORS APRIL 2, 27 In this contet, a vector is a column matri E Note 2 v 2, v 4 5 6 () We might also write v as v Both notations refer to a vector (2) A vector can be man
More informationPolynomial and Rational Functions
Name Date Chapter Polnomial and Rational Functions Section.1 Quadratic Functions Objective: In this lesson ou learned how to sketch and analze graphs of quadratic functions. Important Vocabular Define
More informationArithmetic, Algebra, Number Theory
Arithmetic, Algebra, Number Theory Peter Simon 21 April 2004 Types of Numbers Natural Numbers The counting numbers: 1, 2, 3,... Prime Number A natural number with exactly two factors: itself and 1. Examples:
More informationC) x m A) 260 sq. m B) 26 sq. m C) 40 sq. m D) 364 sq. m. 7) x x - (6x + 24) = -4 A) 0 B) all real numbers C) 4 D) no solution
Sample Departmental Final - Math 46 Perform the indicated operation. Simplif if possible. 1) 7 - - 2-2 + 3 2 - A) + - 2 B) - + 4-2 C) + 4-2 D) - + - 2 Solve the problem. 2) The sum of a number and its
More informationThe Fibonacci Quarterly 44(2006), no.2, PRIMALITY TESTS FOR NUMBERS OF THE FORM k 2 m ± 1. Zhi-Hong Sun
The Fibonacci Quarterly 44006, no., 11-130. PRIMALITY TESTS FOR NUMBERS OF THE FORM k m ± 1 Zhi-Hong Sun eartment of Mathematics, Huaiyin Teachers College, Huaian, Jiangsu 3001, P.R. China E-mail: zhsun@hytc.edu.cn
More informationFinding Divisors, Number of Divisors and Summation of Divisors. Jane Alam Jan
Finding Divisors, Number of Divisors and Summation of Divisors by Jane Alam Jan Finding Divisors If we observe the property of the numbers, we will find some interesting facts. Suppose we have a number
More informationFoundations of Databases
Foundations of Databases (Slides adapted from Thomas Eiter, Leonid Libkin and Werner Nutt) Foundations of Databases 1 Quer optimization: finding a good wa to evaluate a quer Queries are declarative, and
More informationc 2016 Society for Industrial and Applied Mathematics
SIAM J. SCI. COMPUT. Vol. 38, No. 2, pp. A765 A788 c 206 Societ for Industrial and Applied Mathematics ROOTS OF BIVARIATE POLYNOMIAL SYSTEMS VIA DETERMINANTAL REPRESENTATIONS BOR PLESTENJAK AND MICHIEL
More informationCarmichael numbers with a totient of the form a 2 + nb 2
Carmichael numbers with a totient of the form a 2 + nb 2 William D. Banks Department of Mathematics University of Missouri Columbia, MO 65211 USA bankswd@missouri.edu Abstract Let ϕ be the Euler function.
More informationRiemann s explicit formula
Garrett 09-09-20 Continuing to review the simple case (haha!) of number theor over Z: Another example of the possibl-suprising application of othe things to number theor. Riemann s explicit formula More
More informationComputation of Csiszár s Mutual Information of Order α
Computation of Csiszár s Mutual Information of Order Damianos Karakos, Sanjeev Khudanpur and Care E. Priebe Department of Electrical and Computer Engineering and Center for Language and Speech Processing
More informationPartial differentiation. Background mathematics review
Partial differentiation Background mathematics review David Miller Partial differentiation Partial derivatives Background mathematics review David Miller Partial derivative Suppose we have a function f(,)
More informationDiscrete Structures Lecture Primes and Greatest Common Divisor
DEFINITION 1 EXAMPLE 1.1 EXAMPLE 1.2 An integer p greater than 1 is called prime if the only positive factors of p are 1 and p. A positive integer that is greater than 1 and is not prime is called composite.
More informationAN ANALYSIS OF A FRACTAL KINETICS CURVE OF SAVAGEAU
ANZIAM J. 452003, 26 269 AN ANALYSIS OF A FRACTAL KINETICS CURVE OF SAVAGEAU JOHN MALONEY and JACK HEIDEL Received 26 June, 2000; revised 2 Januar, 200 Abstract The fractal kinetics curve derived b Savageau
More informationPrimes in arithmetic progressions to large moduli
CHAPTER 3 Primes in arithmetic progressions to large moduli In this section we prove the celebrated theorem of Bombieri and Vinogradov Theorem 3. (Bombieri-Vinogradov). For any A, thereeistsb = B(A) such
More informationChapter 1. Introduction to prime number theory. 1.1 The Prime Number Theorem
Chapter 1 Introduction to prime number theory 1.1 The Prime Number Theorem In the first part of this course, we focus on the theory of prime numbers. We use the following notation: we write f( g( as if
More informationUnit 3 NOTES Honors Common Core Math 2 1. Day 1: Properties of Exponents
Unit NOTES Honors Common Core Math Da : Properties of Eponents Warm-Up: Before we begin toda s lesson, how much do ou remember about eponents? Use epanded form to write the rules for the eponents. OBJECTIVE
More informationComplex Zeros of the Modified Bessel Function Kn(ZT
MATHEMATICS OF COMPUTATION, VOLUME 26, NUMBER 120, OCTOBER 1972 Complex Zeros of the Modified Bessel Function Kn(ZT B R. Parnes Abstract. The complex zeros of Kn(Z) are computed for integer orders n =
More informationSecond-Order Linear Differential Equations C 2
C8 APPENDIX C Additional Topics in Differential Equations APPENDIX C. Second-Order Homogeneous Linear Equations Second-Order Linear Differential Equations Higher-Order Linear Differential Equations Application
More informationOn the counting function for the generalized Niven numbers
On the counting function for the generalized Niven numbers Ryan Daileda, Jessica Jou, Robert Lemke-Oliver, Elizabeth Rossolimo, Enrique Treviño Abstract Given an integer base q 2 and a completely q-additive
More information2 Ordinary Differential Equations: Initial Value Problems
Ordinar Differential Equations: Initial Value Problems Read sections 9., (9. for information), 9.3, 9.3., 9.3. (up to p. 396), 9.3.6. Review questions 9.3, 9.4, 9.8, 9.9, 9.4 9.6.. Two Examples.. Foxes
More information1. Introduction. Let P and Q be non-zero relatively prime integers, α and β (α > β) be the zeros of x 2 P x + Q, and, for n 0, let
C O L L O Q U I U M M A T H E M A T I C U M VOL. 78 1998 NO. 1 SQUARES IN LUCAS SEQUENCES HAVING AN EVEN FIRST PARAMETER BY PAULO R I B E N B O I M (KINGSTON, ONTARIO) AND WAYNE L. M c D A N I E L (ST.
More informationTrigonometric. equations. Topic: Periodic functions and applications. Simple trigonometric. equations. Equations using radians Further trigonometric
Trigonometric equations 6 sllabusref eferenceence Topic: Periodic functions and applications In this cha 6A 6B 6C 6D 6E chapter Simple trigonometric equations Equations using radians Further trigonometric
More informationOn the counting function for the generalized. Niven numbers
On the counting function for the generalized Niven numbers R. C. Daileda, Jessica Jou, Robert Lemke-Oliver, Elizabeth Rossolimo, Enrique Treviño Introduction Let q 2 be a fied integer and let f be an arbitrary
More informationStudy of some equivalence classes of primes
Notes on Number Theory and Discrete Mathematics Print ISSN 3-532, Online ISSN 2367-8275 Vol 23, 27, No 2, 2 29 Study of some equivalence classes of primes Sadani Idir Department of Mathematics University
More informationMAT 127: Calculus C, Fall 2010 Solutions to Midterm I
MAT 7: Calculus C, Fall 00 Solutions to Midterm I Problem (0pts) Consider the four differential equations for = (): (a) = ( + ) (b) = ( + ) (c) = e + (d) = e. Each of the four diagrams below shows a solution
More informationSpotlight on the Extended Method of Frobenius
113 Spotlight on the Extended Method of Frobenius See Sections 11.1 and 11.2 for the model of an aging spring. Reference: Section 11.4 and SPOTLIGHT ON BESSEL FUNCTIONS. Bessel functions of the first kind
More informationOn Prime Divisors of Binomial Coefficients
MATHEMATICS OF COMPUTATION VOLUME 51, NUMBER 183 JULY 1988, PAGES 325-329 On Prime Divisors of Binomial Coefficients By Pierre Goetgheluck Abstract. This paper, using computational and theoretical methods,
More information18.785: Analytic Number Theory, MIT, spring 2007 (K.S. Kedlaya) Brun s combinatorial sieve
18.785: Analytic Number Theory, MIT, spring 2007 (K.S. Kedlaya) Brun s combinatorial sieve In this unit, we describe a more intricate version of the sieve of Eratosthenes, introduced by Viggo Brun in order
More informationALGEBRA+NUMBER THEORY +COMBINATORICS
ALGEBRA+NUMBER THEORY +COMBINATORICS COMP 321 McGill University These slides are mainly compiled from the following resources. - Professor Jaehyun Park slides CS 97SI - Top-coder tutorials. - Programming
More informationPrimes at a Glance. By R. K. Guy, C. B. Lacampagne, and J. L. Selfridge
mathematics of computation volume 4x. number 177 january 1987, pages 18-2(12 Primes at a Glance By R. K. Guy, C. B. Lacampagne, and J. L. Selfridge To Dan Shanks: May he stay on form und again become the
More informationInitial Value Problems for. Ordinary Differential Equations
Initial Value Problems for Ordinar Differential Equations INTRODUCTION Equations which are composed of an unnown function and its derivatives are called differential equations. It becomes an initial value
More informationToda s Theorem: PH P #P
CS254: Computational Compleit Theor Prof. Luca Trevisan Final Project Ananth Raghunathan 1 Introduction Toda s Theorem: PH P #P The class NP captures the difficult of finding certificates. However, in
More informationThe incomplete gamma functions. Notes by G.J.O. Jameson. These notes incorporate the Math. Gazette article [Jam1], with some extra material.
The incomplete gamma functions Notes by G.J.O. Jameson These notes incorporate the Math. Gazette article [Jam], with some etra material. Definitions and elementary properties functions: Recall the integral
More informationRiemann s ζ-function
Int. J. Contemp. Math. Sciences, Vol. 4, 9, no. 9, 45-44 Riemann s ζ-function R. A. Mollin Department of Mathematics and Statistics University of Calgary, Calgary, Alberta, Canada, TN N4 URL: http://www.math.ucalgary.ca/
More information8 Differential Calculus 1 Introduction
8 Differential Calculus Introduction The ideas that are the basis for calculus have been with us for a ver long time. Between 5 BC and 5 BC, Greek mathematicians were working on problems that would find
More informationProvable intractability. NP-completeness. Decision problems. P problems. Algorithms. Two causes of intractability:
Provable intractabilit -completeness Selected from Computer and Intractabilit b MR Gare and DS Johnson (979) Two causes of intractabilit: The problem is so difficult that an exponential amount of time
More informationCSE 546 Midterm Exam, Fall 2014
CSE 546 Midterm Eam, Fall 2014 1. Personal info: Name: UW NetID: Student ID: 2. There should be 14 numbered pages in this eam (including this cover sheet). 3. You can use an material ou brought: an book,
More informationPure Core 1. Revision Notes
Pure Core Revision Notes Ma 06 Pure Core Algebra... Indices... Rules of indices... Surds... 4 Simplifing surds... 4 Rationalising the denominator... 4 Quadratic functions... 5 Completing the square....
More informationFINITE ABELIAN GROUPS Amin Witno
WON Series in Discrete Mathematics and Modern Algebra Volume 7 FINITE ABELIAN GROUPS Amin Witno Abstract We detail the proof of the fundamental theorem of finite abelian groups, which states that every
More informationKEY IDEAS. Chapter 1 Function Transformations. 1.1 Horizontal and Vertical Translations Pre-Calculus 12 Student Workbook MHR 1
Chapter Function Transformations. Horizontal and Vertical Translations A translation can move the graph of a function up or down (vertical translation) and right or left (horizontal translation). A translation
More informationSEVERAL PROOFS OF THE IRREDUCIBILITY OF THE CYCLOTOMIC POLYNOMIALS
SEVERAL PROOFS OF THE IRREDUCIBILITY OF THE CYCLOTOMIC POLYNOMIALS STEVEN H. WEINTRAUB ABSTRACT. We present a number of classical proofs of the irreducibility of the n-th cyclotomic polynomial Φ n (x).
More informationEULER S THEOREM KEITH CONRAD
EULER S THEOREM KEITH CONRAD. Introduction Fermat s little theorem is an important property of integers to a prime modulus. Theorem. (Fermat). For prime p and any a Z such that a 0 mod p, a p mod p. If
More informationCourse 15 Numbers and Their Properties
Course Numbers and Their Properties KEY Module: Objective: Rules for Eponents and Radicals To practice appling rules for eponents when the eponents are rational numbers Name: Date: Fill in the blanks.
More informationAssociative property
Addition Associative property Closure property Commutative property Composite number Natural numbers (counting numbers) Distributive property for multiplication over addition Divisibility Divisor Factor
More informationmathematics of computation, volume 27, number 124, October 1973 The First Occurrence of Large Gaps Between Successive Primes By Richard P.
mathematics of computation, volume 27, number 124, October 1973 The First Occurrence of Large Gaps Between Successive Primes By Richard P. Brent Abstract. A table of the first occurrence of a string of
More informationUNCORRECTED SAMPLE PAGES. 3Quadratics. Chapter 3. Objectives
Chapter 3 3Quadratics Objectives To recognise and sketch the graphs of quadratic polnomials. To find the ke features of the graph of a quadratic polnomial: ais intercepts, turning point and ais of smmetr.
More informationDefinitions, Theorems and Exercises. Abstract Algebra Math 332. Ethan D. Bloch
Definitions, Theorems and Exercises Abstract Algebra Math 332 Ethan D. Bloch December 26, 2013 ii Contents 1 Binary Operations 3 1.1 Binary Operations............................... 4 1.2 Isomorphic Binary
More information5.6. Differential equations
5.6. Differential equations The relationship between cause and effect in phsical phenomena can often be formulated using differential equations which describe how a phsical measure () and its derivative
More informationRepresentation of Functions by Power Series. Geometric Power Series
60_0909.qd //0 :09 PM Page 669 SECTION 9.9 Representation of Functions b Power Series 669 The Granger Collection Section 9.9 JOSEPH FOURIER (768 80) Some of the earl work in representing functions b power
More informationComputing Bernoulli Numbers Quickly
Computing Bernoulli Numbers Quickly Kevin J. McGown December 8, 2005 The Bernoulli numbers are defined via the coefficients of the power series expansion of t/(e t ). Namely, for integers m 0 we define
More informationTHE JACOBI SYMBOL AND A METHOD OF EISENSTEIN FOR CALCULATING IT
THE JACOBI SYMBOL AND A METHOD OF EISENSTEIN FOR CALCULATING IT STEVEN H. WEINTRAUB ABSTRACT. We present an exposition of the asic properties of the Jacoi symol, with a method of calculating it due to
More informationCullen Numbers in Binary Recurrent Sequences
Cullen Numbers in Binary Recurrent Sequences Florian Luca 1 and Pantelimon Stănică 2 1 IMATE-UNAM, Ap. Postal 61-3 (Xangari), CP 58 089 Morelia, Michoacán, Mexico; e-mail: fluca@matmor.unam.mx 2 Auburn
More informationFACTORS OF GENERALIZED FERMAT NUMBERS
mathematics of computation volume 64, number 20 january, pages -40 FACTORS OF GENERALIZED FERMAT NUMBERS HARVEY DUBNER AND WILFRID KELLER Abstract. Generalized Fermât numbers have the form Fb
More informationPrime Pythagorean triangles
1 2 3 47 6 23 11 Journal of Integer Sequences, Vol. 4 (2001), Article 01.2.3 Prime Pythagorean triangles Harvey Dubner 449 Beverly Road, Ridgewood, New Jersey 07450 Tony Forbes Department of Pure Mathematics,
More informationAlgebra II Notes Unit Six: Polynomials Syllabus Objectives: 6.2 The student will simplify polynomial expressions.
Algebra II Notes Unit Si: Polnomials Sllabus Objectives: 6. The student will simplif polnomial epressions. Review: Properties of Eponents (Allow students to come up with these on their own.) Let a and
More informationLecture 4: Exact ODE s
Lecture 4: Exact ODE s Dr. Michael Doughert Januar 23, 203 Exact equations are first-order ODE s of a particular form, and whose methods of solution rel upon basic facts concerning partial derivatives,
More informationSolving Linear Systems
1.4 Solving Linear Sstems Essential Question How can ou determine the number of solutions of a linear sstem? A linear sstem is consistent when it has at least one solution. A linear sstem is inconsistent
More informationSystems of Linear Equations: Solving by Graphing
8.1 Sstems of Linear Equations: Solving b Graphing 8.1 OBJECTIVE 1. Find the solution(s) for a set of linear equations b graphing NOTE There is no other ordered pair that satisfies both equations. From
More informationChapter 8 Notes SN AA U2C8
Chapter 8 Notes SN AA U2C8 Name Period Section 8-: Eploring Eponential Models Section 8-2: Properties of Eponential Functions In Chapter 7, we used properties of eponents to determine roots and some of
More informationMATH Line integrals III Fall The fundamental theorem of line integrals. In general C
MATH 255 Line integrals III Fall 216 In general 1. The fundamental theorem of line integrals v T ds depends on the curve between the starting point and the ending point. onsider two was to get from (1,
More information2.5 CONTINUITY. a x. Notice that Definition l implicitly requires three things if f is continuous at a:
SECTION.5 CONTINUITY 9.5 CONTINUITY We noticed in Section.3 that the it of a function as approaches a can often be found simpl b calculating the value of the function at a. Functions with this propert
More informationExponentials and Logarithms Review Part 2: Exponentials
Eponentials and Logaritms Review Part : Eponentials Notice te difference etween te functions: g( ) and f ( ) In te function g( ), te variale is te ase and te eponent is a constant. Tis is called a power
More informationChapter 1. Introduction to prime number theory. 1.1 The Prime Number Theorem
Chapter 1 Introduction to prime number theory 1.1 The Prime Number Theorem In the first part of this course, we focus on the theory of prime numbers. We use the following notation: we write f g as if lim
More informationSPEEDING FERMAT S FACTORING METHOD
MATHEMATICS OF COMPUTATION Volume 68, Number 228, Pages 1729 1737 S 0025-5718(99)01133-3 Article electronically published on March 1, 1999 SPEEDING FERMAT S FACTORING METHOD JAMES MCKEE Abstract. A factoring
More informationHarmonic sets and the harmonic prime number theorem
Harmonic sets and the harmonic prime number theorem Version: 9th September 2004 Kevin A. Broughan and Rory J. Casey University of Waikato, Hamilton, New Zealand E-mail: kab@waikato.ac.nz We restrict primes
More informationSOLVING FERMAT-TYPE EQUATIONS x 5 + y 5 = dz p
MATHEMATICS OF COMPUTATION Volume 79, Number 269, January 2010, Pages 535 544 S 0025-5718(09)02294-7 Article electronically published on July 22, 2009 SOLVING FERMAT-TYPE EQUATIONS x 5 + y 5 = dz p NICOLAS
More informationINSTRUCTOR S SOLUTIONS MANUAL SECTION 17.1 (PAGE 902)
INSTRUCTOR S SOLUTIONS MANUAL SECTION 7 PAGE 90 CHAPTER 7 ORDINARY DIFFEREN- TIAL EQUATIONS Section 7 Classifing Differential Equations page 90 = 5: st order, linear, homogeneous d d + = : nd order, linear,
More informationDefinition: div Let n, d 0. We define ndiv d as the least integer quotient obtained when n is divided by d. That is if
Section 5. Congruence Arithmetic A number of computer languages have built-in functions that compute the quotient and remainder of division. Definition: div Let n, d 0. We define ndiv d as the least integer
More informationDETERMINISTIC METHODS TO FIND PRIMES
DETERMINISTIC METHODS TO FIND PRIMES D.H.J. POLYMATH Abstract. Given a large positive integer N, how quickly can one construct a prime number larger than N (or between N and 2N)? Using probabilistic methods,
More informationSome elementary explicit bounds for two mollifications of the Moebius function
Some elementary explicit bounds for two mollifications of the Moebius function O. Ramaré August 22, 213 Abstract We prove that the sum { /d 1+ε is bounded by, uniformly in x 1, r and ε >. We prove a similar
More information