Cyclotomic polynomial

From HandWiki
Short description: Irreducible polynomial whose roots are nth roots of unity

In mathematics, the n-th cyclotomic polynomial, for any positive integer n, is the unique irreducible polynomial with integer coefficients that is a divisor of xn1 and is not a divisor of xk1 for any k<n. Its roots are all n-th primitive roots of unity e2iπkn, where k runs over the positive integers less than n and coprime to n (where i is the imaginary unit). In other words, the n-th cyclotomic polynomial is equal to

Φn(x)=gcd(k,n)=11kn1(xe2iπkn).

It may also be defined as the monic polynomial with integer coefficients that is the minimal polynomial over the field of the rational numbers of any primitive nth-root of unity (e2iπ/n is an example of such a root).

An important relation linking cyclotomic polynomials and primitive roots of unity is

dnΦd(x)=xn1,

showing that x is a root of xn1 if and only if it is a d-th primitive root of unity for some d that divides n.Cite error: Closing </ref> missing for <ref> tag

Φ1(x)=x1Φ2(x)=x+1Φ3(x)=x2+x+1Φ4(x)=x2+1Φ5(x)=x4+x3+x2+x+1Φ6(x)=x2x+1Φ7(x)=x6+x5+x4+x3+x2+x+1Φ8(x)=x4+1Φ9(x)=x6+x3+1Φ10(x)=x4x3+x2x+1Φ11(x)=x10+x9+x8+x7+x6+x5+x4+x3+x2+x+1Φ12(x)=x4x2+1Φ13(x)=x12+x11+x10+x9+x8+x7+x6+x5+x4+x3+x2+x+1Φ14(x)=x6x5+x4x3+x2x+1Φ15(x)=x8x7+x5x4+x3x+1Φ16(x)=x8+1Φ17(x)=x16+x15+x14+x13+x12+x11+x10+x9+x8+x7+x6+x5+x4+x3+x2+x+1Φ18(x)=x6x3+1Φ19(x)=x18+x17+x16+x15+x14+x13+x12+x11+x10+x9+x8+x7+x6+x5+x4+x3+x2+x+1Φ20(x)=x8x6+x4x2+1Φ21(x)=x12x11+x9x8+x6x4+x3x+1Φ22(x)=x10x9+x8x7+x6x5+x4x3+x2x+1Φ23(x)=x22+x21+x20+x19+x18+x17+x16+x15+x14+x13+x12+x11+x10+x9+x8+x7+x6+x5+x4+x3+x2+x+1Φ24(x)=x8x4+1Φ25(x)=x20+x15+x10+x5+1Φ26(x)=x12x11+x10x9+x8x7+x6x5+x4x3+x2x+1Φ27(x)=x18+x9+1Φ28(x)=x12x10+x8x6+x4x2+1Φ29(x)=x28+x27+x26+x25+x24+x23+x22+x21+x20+x19+x18+x17+x16+x15+x14+x13+x12+x11+x10+x9+x8+x7+x6+x5+x4+x3+x2+x+1Φ30(x)=x8+x7x5x4x3+x+1.

The case of the 105th cyclotomic polynomial is interesting because 105 is the least positive integer that is the product of three distinct odd prime numbers (3×5×7) and this polynomial is the first one that has a coefficient other than 1, 0, or −1:[1]

Φ105(x)=x48+x47+x46x43x422x41x40x39+x36+x35+x34+x33+x32+x31x28x26x24x22x20+x17+x16+x15+x14+x13+x12x9x82x7x6x5+x2+x+1.

Properties

Fundamental tools

The cyclotomic polynomials are monic polynomials with integer coefficients that are irreducible over the field of the rational numbers. Except for n equal to 1 or 2, they are palindromes of even degree.

The degree of Φn, or in other words the number of nth primitive roots of unity, is φ(n), where φ is Euler's totient function.

The fact that Φn is an irreducible polynomial of degree φ(n) in the ring [x] is a nontrivial result due to Gauss.[2] Depending on the chosen definition, it is either the value of the degree or the irreducibility which is a nontrivial result. The case of prime n is easier to prove than the general case, thanks to Eisenstein's criterion.

A fundamental relation involving cyclotomic polynomials is

xn1=1kn(xe2iπkn)=dn1kngcd(k,n)=d(xe2iπkn)=dnΦnd(x)=dnΦd(x).

which means that each n-th root of unity is a primitive d-th root of unity for a unique d dividing n.

The Möbius inversion formula allows Φn(x) to be expressed as an explicit rational fraction:

Φn(x)=dn(xd1)μ(nd),

where μ is the Möbius function.

This provides a recursive formula for the cyclotomic polynomial Φn(x), which may be computed by dividing xn1 by the cyclotomic polynomials Φd(x) for the proper divisors d dividing n, starting from Φ1(x)=x1:

Φn(x)=xn1d<nd|nΦd(x).

This gives an algorithm for computing any Φn(x), provided integer factorization and division of polynomials are available. Many computer algebra systems, such as SageMath, Maple, Mathematica, and PARI/GP, have a built-in function to compute the cyclotomic polynomials.

Easy cases for computation

As noted above, if n = p is a prime number, then

Φp(x)=1+x+x2++xp1=k=0p1xk.

If n is an odd integer greater than one, then

Φ2n(x)=Φn(x).

In particular, if n = 2p is twice an odd prime, then (as noted above)

Φ2p(x)=1x+x2+xp1=k=0p1(x)k.

If n = pm is a prime power (where p is prime), then

Φpm(x)=Φp(xpm1)=k=0p1xkpm1.

More generally, if n = pmr with r relatively prime to p, then

Φpmr(x)=Φpr(xpm1).

These formulas may be applied repeatedly to get a simple expression for any cyclotomic polynomial Φn(x) in terms of a cyclotomic polynomial of square free index: If q is the product of the prime divisors of n (its radical), thenCite error: Closing </ref> missing for <ref> tag

Φnp(x)=Φn(xp)/Φn(x).

Integers appearing as coefficients

The problem of bounding the magnitude of the coefficients of the cyclotomic polynomials has been the object of a number of research papers.[3]

If n has at most two distinct odd prime factors, then Migotti showed that the coefficients of Φn are all in the set {1, −1, 0}.[4]

The first cyclotomic polynomial for a product of three different odd prime factors is Φ105(x); it has a coefficient −2 (see above). The converse is not true: Φ231(x)=Φ3×7×11(x) only has coefficients in {1, −1, 0}.

If n is a product of more different odd prime factors, the coefficients may increase to very high values. E.g., Φ15015(x)=Φ3×5×7×11×13(x) has coefficients running from −22 to 23; also Φ255255(x)=Φ3×5×7×11×13×17(x), the smallest n with 6 different odd primes, has coefficients of magnitude up to 532.

Let A(n) denote the maximum absolute value of the coefficients of Φn(x). It is known that for any positive k, the number of n up to x with A(n) > nk is at least c(k)⋅x for a positive c(k) depending on k and x sufficiently large. In the opposite direction, for any function ψ(n) tending to infinity with n we have A(n) bounded above by nψ(n) for almost all n.[5]

A combination of theorems of Bateman and Vaughan states that[3]: 10  on the one hand, for every ε>0, we have

A(n)<e(n(log2+ε)/(loglogn))

for all sufficiently large positive integers n, and on the other hand, we have

A(n)>e(n(log2)/(loglogn))

for infinitely many positive integers n. This implies in particular that univariate polynomials (concretely xn1 for infinitely many positive integers n) can have factors (like Φn) whose coefficients are superpolynomially larger than the original coefficients. This is not too far from the general Landau-Mignotte bound.

Gauss's formula

Let n be odd, square-free, and greater than 3. Then:[6][7]

4Φn(z)=An2(z)(1)n12nz2Bn2(z)

for certain polynomials An(z) and Bn(z) with integer coefficients, An(z) of degree φ(n)/2, and Bn(z) of degree φ(n)/2 − 2. Furthermore, An(z) is palindromic when its degree is even; if its degree is odd it is antipalindromic. Similarly, Bn(z) is palindromic unless n is composite and n ≡ 3 (mod 4), in which case it is antipalindromic.

The first few cases are

4Φ5(z)=4(z4+z3+z2+z+1)=(2z2+z+2)25z24Φ7(z)=4(z6+z5+z4+z3+z2+z+1)=(2z3+z2z2)2+7z2(z+1)2[6pt]4Φ11(z)=4(z10+z9+z8+z7+z6+z5+z4+z3+z2+z+1)=(2z5+z42z3+2z2z2)2+11z2(z3+1)2

Lucas's formula

Let n be odd, square-free and greater than 3. Then[7]

Φn(z)=Un2(z)(1)n12nzVn2(z)

for certain polynomials Un(z) and Vn(z) with integer coefficients, Un(z) of degree φ(n)/2, and Vn(z) of degree φ(n)/2 − 1. This can also be written

Φn((1)n12z)=Cn2(z)nzDn2(z).

If n is even, square-free and greater than 2 (this forces n/2 to be odd),

Φn2(z2)=Φ2n(z)=Cn2(z)nzDn2(z)

for Cn(z) and Dn(z) with integer coefficients, Cn(z) of degree φ(n), and Dn(z) of degree φ(n) − 1. Cn(z) and Dn(z) are both palindromic.

The first few cases are:

Φ3(z)=Φ6(z)=z2z+1=(z+1)23zΦ5(z)=z4+z3+z2+z+1=(z2+3z+1)25z(z+1)2Φ6/2(z2)=Φ12(z)=z4z2+1=(z2+3z+1)26z(z+1)2

Sister Beiter conjecture

The Sister Beiter conjecture is concerned with the maximal size (in absolute value) A(pqr) of coefficients of ternary cyclotomic polynomials Φpqr(x) where pqr are three odd primes.[8]

Cyclotomic polynomials over a finite field and over the p-adic integers

Over a finite field with a prime number p of elements, for any integer n that is not a multiple of p, the cyclotomic polynomial Φn factorizes into φ(n)d irreducible polynomials of degree d, where φ(n) is Euler's totient function and d is the multiplicative order of p modulo n. In particular, Φn is irreducible if and only if p is a primitive root modulo n, that is, p does not divide n, and its multiplicative order modulo n is φ(n), the degree of Φn.[9]

These results are also true over the p-adic integers, since Hensel's lemma allows lifting a factorization over the field with p elements to a factorization over the p-adic integers.

Polynomial values

If x takes any real value, then Φn(x)>0 for every n ≥ 3 (this follows from the fact that the roots of a cyclotomic polynomial are all non-real, for n ≥ 3).

For studying the values that a cyclotomic polynomial may take when x is given an integer value, it suffices to consider only the case n ≥ 3, as the cases n = 1 and n = 2 are trivial (one has Φ1(x)=x1 and Φ2(x)=x+1).

For n ≥ 2, one has

Φn(0)=1,
Φn(1)=1 if n is not a prime power,
Φn(1)=p if n=pk is a prime power with k ≥ 1.

The values that a cyclotomic polynomial Φn(x) may take for other integer values of x is strongly related with the multiplicative order modulo a prime number.

More precisely, given a prime number p and an integer b coprime with p, the multiplicative order of b modulo p, is the smallest positive integer n such that p is a divisor of bn1. For b > 1, the multiplicative order of b modulo p is also the shortest period of the representation of 1/p in the numeral base b (see Unique prime; this explains the notation choice).

The definition of the multiplicative order implies that, if n is the multiplicative order of b modulo p, then p is a divisor of Φn(b). The converse is not true, but one has the following.

If n > 0 is a positive integer and b > 1 is an integer, then (see below for a proof)

Φn(b)=2kgh,

where

  • k is a non-negative integer, always equal to 0 when b is even. (In fact, if n is neither 1 nor 2, then k is either 0 or 1. Besides, if n is not a power of 2, then k is always equal to 0)
  • g is 1 or the largest odd prime factor of n.
  • h is odd, coprime with n, and its prime factors are exactly the odd primes p such that n is the multiplicative order of b modulo p.

This implies that, if p is an odd prime divisor of Φn(b), then either n is a divisor of p − 1 or p is a divisor of n. In the latter case, p2 does not divide Φn(b).

Zsigmondy's theorem implies that the only cases where b > 1 and h = 1 are

Φ1(2)=1Φ2(2k1)=2kk>0Φ6(2)=3

It follows from above factorization that the odd prime factors of

Φn(b)gcd(n,Φn(b))

are exactly the odd primes p such that n is the multiplicative order of b modulo p. This fraction may be even only when b is odd. In this case, the multiplicative order of b modulo 2 is always 1.

There are many pairs (n, b) with b > 1 such that Φn(b) is prime. In fact, Bunyakovsky conjecture implies that, for every n, there are infinitely many b > 1 such that Φn(b) is prime. See OEISA085398 for the list of the smallest b > 1 such that Φn(b) is prime (the smallest b > 1 such that Φn(b) is prime is about γφ(n), where γ is Euler–Mascheroni constant, and φ is Euler's totient function). See also OEISA206864 for the list of the smallest primes of the form Φn(b) with n > 2 and b > 1, and, more generally, OEISA206942, for the smallest positive integers of this form.


Applications

Using Φn, one can give an elementary proof for the infinitude of primes congruent to 1 modulo n,[10] which is a special case of Dirichlet's theorem on arithmetic progressions.

Periodic recursive sequences

The constant-coefficient linear recurrences which are periodic are precisely the power series coefficients of rational functions whose denominators are products of cyclotomic polynomials.

In the theory of combinatorial generating functions, the denominator of a rational function determines a linear recurrence for its power series coefficients. For example, the Fibonacci sequence has generating function

F(x)=F1x+F2x2+F3x3+=x1xx2,

and equating coefficients on both sides of

F(x)(1xx2)=x

gives

FnFn1Fn2=0

for

n2

.

Any rational function whose denominator is a divisor of xn1 has a recursive sequence of coefficients which is periodic with period at most n. For example,

P(x)=1+2xΦ6(x)=1+2x1x+x2=n0Pnxn=1+3x+2x2x33x42x5+x6+3x7+2x8+

has coefficients defined by the recurrence

PnPn1+Pn2=0

for

n2

, starting from

P0=1,P1=3

. But

1x6=Φ6(x)Φ3(x)Φ2(x)Φ1(x)

, so we may write

P(x)=(1+2x)Φ3(x)Φ2(x)Φ1(x)1x6=1+3x+2x2x33x42x51x6,

which means

PnPn6=0

for

n6

, and the sequence has period 6 with initial values given by the coefficients of the numerator.

See also

References

  1. Brookfield, Gary (2016), "The coefficients of cyclotomic polynomials", Mathematics Magazine 89 (3): 179–188, doi:10.4169/math.mag.89.3.179 
  2. Lang, Serge (2002), Algebra, Graduate Texts in Mathematics, 211 (Revised third ed.), New York: Springer-Verlag, ISBN 978-0-387-95385-4 
  3. 3.0 3.1 Sanna, Carlo (2021), "A Survey on Coefficients of Cyclotomic Polynomials", arXiv:2111.04034 [math.NT]
  4. Isaacs, Martin (2009), Algebra: A Graduate Course, AMS Bookstore, p. 310, ISBN 978-0-8218-4799-2 
  5. Maier, Helmut (2008), "Anatomy of integers and cyclotomic polynomials", in De Koninck, Jean-Marie; Granville, Andrew; Luca, Florian, Anatomy of integers. Based on the CRM workshop, Montreal, Canada, March 13-17, 2006, CRM Proceedings and Lecture Notes, 46, Providence, RI: American Mathematical Society, pp. 89–95, ISBN 978-0-8218-4406-9 
  6. Gauss, DA, Articles 356-357
  7. 7.0 7.1 Riesel, Hans (1994), Prime Numbers and Computer Methods for Factorization (2nd ed.), Boston: Birkhäuser, pp. 309-316, 436, 443, ISBN 0-8176-3743-5 
  8. Beiter, Marion (April 1968), "Magnitude of the Coefficients of the Cyclotomic Polynomial Fpqr(x)", The American Mathematical Monthly 75 (4): 370–372, doi:10.2307/2313416 
  9. Lidl, Rudolf; Niederreiter, Harald (2008), Finite Fields (2nd ed.), Cambridge University Press, p. 65 .
  10. S. Shirali. Number Theory. Orient Blackswan, 2004. p. 67. ISBN 81-7371-454-1

Further reading

Gauss's book Disquisitiones Arithmeticae [Arithmetical Investigations] has been translated from Latin into French, German, and English. The German edition includes all of his papers on number theory: all the proofs of quadratic reciprocity, the determination of the sign of the Gauss sum, the investigations into biquadratic reciprocity, and unpublished notes.