On the Fourier transform of the greatest common divisor arxiv: v1 [math.nt] 16 Jan 2012

Transcription

1 On the Fourier transform of the greatest common ivisor arxiv: v1 [math.nt] 16 Jan 2012 Peter H. van er Kamp Department of Mathematics an Statistics La Trobe University Victoria 3086, Australia January 17, 2012 Abstract The iscrete Fourier transform of the greatest common ivisor m î[a](m) gc(, m)αm a, with α m a primitive m-th root of unity, is a multiplicative function that generalises both the gc-sum function an Euler s totient function. On the one han it is the Dirichlet convolution of the ientity with Ramanujan s sum, î[a] i c (a), an on the other han it can be written as a generalise convolution prouct, î[a] i a φ. We show that î[a](m) counts the number of elements in the set of orere pairs (i, j) such that i j a mo m. Furthermore we generalise a ozen nown ientities for the totient function, to ientities which involve the iscrete Fourier transform of the greatest common ivisor, incluing its partial sums, an its Lambert series. 1 Introuction In [15] iscrete Fourier transforms of functions of the greatest common ivisor were stuie, i.e. m ĥ[a](m) h(gc(, m))αm a, where α m is a primitive m-th root of unity. The main result in that paper is the ientity 1 ĥ[a] h c (a), where enote Dirichlet convolution, i.e. ĥ[a](m) m h( m )c (a), (1) 1 Similar results in the context of r-even function were obtaine earlier, see [10] for etails. 1

2 an c m (a) m gc(,m)1 α a m (2) enotes Ramanujan s sum. Ramanujan s sum generalises both Euler s totient function φ c (0) an the Möbius function µ c (1). Thus, ientity (1) generalizes the formula m h(gc(, m)) (h φ)(m), (3) alreay nown to Cesàro in The formula (1) shows that ĥ[a] is multiplicative if h is multiplicative (because c (a) is multiplicative an Dirichlet convolution preserves multiplicativity). Here we will tae h i to be the ientity function (of the greatest common ivisor) an stuy its Fourier transform. Obviously, as i : i(n) n is multiplicative, the function î[a] is multiplicative, for all a. Two special cases are well-nown. Taing a 0 we have î[0] P, where P(m) m gc(, m). (4) is nown as Pillai s arithmetical function or the gc-sum function. Seconly, by taing a 1 in (1), we fin that î[1] i µ equals φ, by Möbius inversion of Euler s ientity φ u i, where u µ 1 is the unit function efine by u(m) 1. Let Fa m enote the set of orere pairs (i, j) such that i j a mo m, the set of factorizations of a moulo m. We claim that î[a](m) counts its number of elements. Let us consier the mentione special cases. a 0 For given i {1, 2,..., m} the congruence i j 0 mo m yiels i gc(i, m) j 0 mo m gc(i, m), which has a unique solution moulo m/gc(i, m), an so there are gc(i, m) solutions moulo m. Hence, the total number of elements in F0 m is P(m). a 1 The totient function φ(m) counts the number of invertible congruence classes moulo m. As for every invertible congruence class i moulo m there is a unique j i 1 mo m such that i j 1 mo m, it counts the number of elements in the set F1 m. To prove the general case we employ a Kluyver-lie formula for î[a], that is, a formula similar to the formula for the Ramanujan sum function c (a) gc(a,) µ( ). (5) 2

3 attribute to Kluyver. Together the ientities (1) an (5) imply, cf. section 3, î[a](m) gc(a,m) φ( m ), (6) an we will show, in the next section, that the number of factorizations of a mo m is given by the same sum. The right han sies of (5) an (6) are particular instances of the so calle generalize Ramanujan sums [1], an both formulas follow as consequence of a general formula for the Fourier coefficients of these generalise Ramanujan sums [2, 3]. In section 3 we provie simple proofs for some of the nice properties of these sums. In particular we interpret the sums as a generalization of Dirichlet convolution. This interpretation lies at the heart of many of the generalise totient ientities we establish in section 4. 2 The number of factorizations of a mo m For given i, m N, enote g gc(i, m). If the congruence i j a mo m has a solution j, then g a an j i 1 a/g is unique mo m/g, so mo m there are g solutions. This yiels which can be written as #F m a m i1 gc(i,m) a #F m a a gc(i, m), m i1 gc(i,m) (7) If m then the sum m i1 gc(i,m) 1 is empty. Now let m. The only integers i which contribute to the sum are the multiples of,, where gc(, m/) 1. There are exactly φ(m/) of them. Therefore the right han sies of formulae (6) an (7) agree, an hence #F m a î[a](m). 3 A historical remar, an generalise Ramanujan sums It is well nown that Ramanujan was not the first who consiere the sum c m (a). Kluyver prove his formula (5) in 1906, twelve years before Ramanujan 3

4 publishe the novel iea of expressing arithmetical functions in the form of a series s a sc s (n) [11]. It is not well nown that Kluyver actually showe that c m (a) equals Von Sternec s function, introuce in [13], i.e. c m (a) µ( m gc(a,m) )φ(m) φ( m gc(a,m) ). (8) This relation is referre to in the literature as Höler s relation, cf. the remar on page 213 in [1]. However, Höler publishe this relation thirty years after Kluyver [7]. We refer to [1, Theorem 2], or [2, Theorem 8.8] for a generalisation of (8). The so calle generalize Ramanujan sums, f a g(m) f()g( m ), (9) gc(a,m) were introuce in [1]. The notation a is new, the sums are enote S(a; m) in [1], s m (a) in [2], an S f,h (a, m) in [3]. In the context of r-even functions [10] the sums are enote S f,g (a), an consiere as sequences of m-even functions, with argument a. We consier the sums as a sequence of functions with argument m, labele by a. We call f a g the a-convolution of f an g. The concept of a-convolution is a generalization of Dirichlet convolution as f 0 g f g. As we will see below, the function f a g is multiplicative (for all a) if f an g are, an the following inter-associative property hols, cf. [3, Theorem 4]. (f a g) h f a (g h). (10) We also aopt the notation f a i a f, an call this the Kluyver, or a- extension of f. Thus, we have f 0 i f, f 1 f, an formulas (5) an (6) become c m (a) µ a (m), an î[a] φ a, respectively. The ientity function I, efine by f I f, is given by I() [ 1], where the Iverson bracet is, with P a logical statement, { 1 if P, [P ] 0 if not P. Let us consier the function f a I. It is f a I() f()[ ] [ a]f(). gc(a,) Since the function [ a] is multiplicative, the function f a I is multiplicative if f is multiplicative. Also, we may write, cf. [3, eq. (9)], f a g(m) m[ a]f()g( m ) (f a I) g(m), which shows that f a g is multiplicative if f an g are. Also, the interassociativity property (10) now easily follows from the associativity of the Dirichlet convolution, (f a g) h ((f a I) g) h (f a I) (g h) f a (g h). 4

5 We note that the a-convolution prouct is neither associative, nor commutative. The inter-associativity an the commutativity of Dirichlet convolution imply that f a g (f g) a f g a. (11) Formula (6) states that the Fourier transform of the greatest common ivisor is the Kluyver extension of the totient function. We provie a simple proof. Proof [of (6)] Employing (1), (5) an (11) we have î[a] i c (a) i µ a (i µ) a φ a. The formula (6) also follows as a special case of the following formula for the Fourier coefficients of a-convolutions, f a g(m) m h (m)α a m, h g f i, (12) given in [1, 2]. The formula (12) combines with (1) an (5) to yiel a formula for functions of the greatest common ivisor, h[] : m h(gc(, m)), namely h[] (h µ) u. (13) Proof [of (13)] The Fourier coefficients of ĥ[a](m) are h[](m). But ĥ[a] h c (a) (i a µ) h i a (µ h), an so, using (12), the Fourier coefficients are also given by (h µ) u(m). For a Dirichlet convolution with a Fourier transform of a function of the greatest common ivisor we have f ĝ[a] f g[a]. (14) Proof [of (14)] f ĝ[a] f (g µ a ) (f g) µ a f g[a] Similarly, for an a-convolution with a Fourier transform of a function of the greatest common ivisor, f a ĝ[b] f a g[b]. (15) Proof [of (15)] f a ĝ[b] f a (g µ b ) (f a g) µ b f a g[b] 4 Generalise totient ientities The totient function is an important function in number theory, an relate fiels of mathematics. It is extensively stuie, connecte to many other notions an functions, an there exist numerous generalisation an extensions, cf. the chapter The many facets of Euler s totient in [12]. The Kluyver extension of the totient function is a very natural extension, an it is most surprising it has not been stuie before. In this section we generalise a ozen nown ientities for the totient function φ, to ientities which involve its Kluyver extension φ a, a..a. the iscrete Fourier transform of the greatest common ivisor. This inclues a generalisation of Euler s ientity, the partial sums of φ a, an its Lambert series. 5

6 4.1 The value of φ a at powers of primes We start by proviing a formula for the value of φ a at powers of primes. This epens only on the multiplicity of the prime in a. The formulae, with p prime, P(p ) ( + 1)p p 1, φ(p ) p p 1, of which the first one is Theorem 2.2 in [4], generalise to φ a (p ) { (p p 1 )(l + 1) if l <, ( + 1)p p 1 if l, (16) where l is the largest integer, or infinity, such that p l a. Proof [of (16)] We have which equals (16). φ a (p ) gc(p l,p ) min(l,) r0 { l 4.2 Partial sums of φ a /i To generalise the totient ientity φ( p ) p r φ(p r ) r0 p p 1 if l <, ( 1 r0 p p 1 ) + p if l, φ() to an ientity for φ a we first establish f 0 () µ() n. (17) f() n. (18) Proof [of (18)] Since there are n/ multiples of in the range [1, n] it follows that f i() f() f() n. As a corollary we obtain f a g 0 () 6 1 f a g() n. (19)

7 Proof [of (19)] Employing (10) we fin f a (g i)() (f a g) i() Now taing f i an g µ in (19) we fin φ a () f a g() n. c (a) n. (20) 4.3 Partial sums of P a /i expresse in terms of φ a Taing f i an g φ in (19) we fin P a () φ a () n. (21) Note that by taing either a 0 in (20), or a 1 in (21), we fin an ientity involving the totient function an the gc-sum function, 4.4 Partial sums of φ a P() φ() n. (22) To generalise the totient ientity, with n > 0, ( ) φ() µ() n 2 2, (23) we first establish ( f 0 () 1 n f() n ) f u(). (24) Proof [of (24)] We have, by changing variable l, (2f i f u)() f()( 2 1) 1 n/ l1 f() n 2. 1 f()(2l 1) 7

8 Note that this gives a nice proof of (23), taing f µ, as n I() [ > 0]. When f µ a, then (11) implies f i φ a, an f u I a, an therefore as a special case of (24) we obtain φ a () 1 [ n] + c (a) n 2 2. (25) a We remar that when n a we have a [ n] σ(a), where σ i u is the sum of ivisors function. 4.5 Generalisation of Euler s ientity Euler s ientity, φ u i, generalises to φ a () τ(gc(a, m))m, (26) m where τ u u is the number of ivisors function. Proof [of (26)] We have φ a u (φ u) a i a where i a (m) gc(a,m) m mτ(gc(a, m)). (27) 4.6 Partial sums of P a expresse in terms of φ a (an τ) If f φ a, then f i P a, an (24) becomes, using (26), ( P a () 1 n ) τ(gc(a, )) + φ a () n 2 2. (28) 4.7 Four ientities of Césaro Accoring to Dicson [5] the following three ientities were obtaine by Césaro: n n φ( n ) P(n), (29) n φ() φ()φ( n ) n n j1 n j1 1 gc(j, n), (30) φ(gc(j, n)). (31) 8

9 Ientity (29), which is Theorem 2.3 in [4], is obtaine by taing a 0 in (6), or h i in (3). It generalises to φ a ( n ) P a(n). (32) n Proof [of (32)] By taing f φ an g i in (11). Ientity (30) is obtaine by taing h 1/i in (3) an generalises to n n φ a() j1 gc(a,n) Proof [of (33)] By taing f φ an g 1/i in (11). 1 gc(j, n (33) ), Ientity (31) is also a special case of (3), with h φ. It generalises to φ a ()φ b ( n ) n Proof [of (34)] We have n j1 gc(a,n) φ b (gc(j, n )). (34) (i a φ) (i b φ) i a (i b (φ φ)) i a (i b φ[0]) i a φb [0], an evaluation at m yiels gc(a,m) m/ j1 φ b (gc(j, m )) The more general ientity (3) generalises to gc(a,m) j1 m φ b (gc(j, m )). m h a (gc(, m)) h φ a (m). (35) 4.8 Three ientities of Liouville Dicson [5, p ] states, amongst many others ientities that were presente by Liouville in the series [9], the following φ()τ( m ) σ(m), (36) m φ()σ[n + 1]( m ) mσ[n](m), (37) m m φ()τ( m2 2 ) θ( m ), (38) m 9

10 where σ[n] i[n] u, i[n](m) m n, an θ(m) is the number of ecompositions of m into two relatively prime factors. All three are of the form φ f g an therefore they gain significance ue to (3), thought Liouville might not have been aware of this. For example, (3) an (36) combine to yiel m τ(gc(, m)) σ(m). The three ientities are easily proven by substituting τ u u, σ[n] i[n] u, τ i[2] θ u, φ µ i, an using µ u I. They generalise to φ a ()τ( m ) σ a(m), (39) m φ a ()σ[n + 1]( m ) mu a σ[n](m), (40) m m φ a ()τ( m2 2 ) τ(gc(a, ))θ( m ). (41) m These generalisation are proven using the same substitutions, together with (11), or for the latter ientity, (10) an (27). 4.9 One ientity of Dirichlet Dicson [5] writes that Dirichlet [6], by taing partial sums on both sies of Euler s ientity, obtaine n ( ) n + 1 φ(). 2 By taing partial sums on both sies of equation (26) we obtain n φ a() ( n + 1 ). (42) 2 a Proof [of (42)] Summing the left han sie of (26) over m yiels φ a () m1 m n φ a() 1 an summing the right han sie of (26) over m yiels τ(gc(a, m))m m1 m1 gc(a,m) m a n/ a n ( n + 1 ) /2. 10

11 4.10 The Lambert series of φ a As shown by Liouville [9], cf. [5, p.120], the Lambert series of the totient function is given by x m φ(m) 1 x m x (1 x) 2. m1 The Lambert series for φ a is given by x m φ a (m) 1 x m p[a](x) x (1 x a ) 2, (43) m1 where the coefficients of p[a](x) 2a c[a]()x 1 are given by c[a] i a t[a], an t[a] is the piece-wise linear function t[a](n) a n a. The polynomials p[a] seem to be irreucible over Z an their zeros are in some sense close to the a-th roots of unity, see Figures 1 an 2. Figure 1: The roots of p[37] are epicte as crosses an the 37 th roots of unity as points. This figure shows that when a is prime the roots of p[a] that are close to 1 are closer to a th roots of unity. Proof [of (43)] Cesàro prove [5] x n f(n) 1 x n x n f(), n1 n n1 cf. exercise 31 to chapter 2 in [14]. By substituting (26) in this formula we fin x n φ a (n) 1 x n x n τ(gc(a, n))n. n1 n1 11

12 Multiplying the right han sie by (1 2x a + x 2a ) yiels ( x n τ(gc(a, n))n) 2( n1 + ( n1+2a na+1 x n τ(gc(a, n))(n 2a)) x n τ(gc(a, n))(n a)) c[a](n)x n, where τ(gc(a, n))n 0 < n a, c[a](n) τ(gc(a, n))(n 2(n a)) τ(gc(a, n))(2a n) a < n 2a, τ(gc(a, n))(n 2(n a) + n 2a) 0 n > 2a. Rewriting, using (27), the fact that gc(a, a + ) gc(a, a ), an iviing by x, yiels the result. n1 Figure 2: The roots of p[35] are epicte as crosses an the 35 th roots of unity as points. This figure shows that roots of p[a] are closest to primitive a th roots of unity A series relate to the Lambert series of φ a Liouville [9] also showe x m φ(m) 1 + x m (1 + x x2 ) (1 x 2 ) 2. m1 We show that x m φ a (m) 1 + x m q[a](x) x (1 x 2a ) 2, (44) m1 12

13 where q[a](x) 4a b[a]()x 1, with b[a] h[a] t[2a], h[a]() i a () 2[2 ]i a (/2). At the en of this section we show that 1 + x 2 ivies q[a](x) if a is o. Proof [of (44)] As the left han sie of (44) is obtaine from the left han sie of (43) by subtracting twice the same series with x replace by x 2, the same is true for the right han sie. Thus it follows that q[a](x) p[a](x)(1 + x a ) 2 2p[a](x 2 )x, an hence, that i a () 2[2 ]i a (/2) a, i b[a]() a (2a ) + 2i a ( a) 2[2 ]i a (/2) a < 2a, 2i a (3a ) + i a ( 2a) 2[2 ]i a (2a /2) 2a < 3a, i a (4a ) 2[2 ]i a (2a /2) 3a < 4a. The result follows ue to the ientities i a (2a ) + 2i a ( a) i a (), 2i a (3a ) + i a ( 2a) i a (4a ), which are easily verifie using (27). We can express the functions b[a] an h[a] in terms of an interesting fractal function. Let a function κ of two variables be efine recursively by 0 2 n, 2 a, or n 0, κ[a](n) κ[a/2](n/2) 2 n, 2 a, (45) τ(gc(a, n)) 2 n. The following properties are easily verifie using the efinition. We have an, with gc(a, b) 1, κ[a](2a + n) κ[a](2a n), (46) κ[an](bn) { 0 2 b, α(n) 2 b, (47) where α enotes the number of o ivisors function, i.e. for all an o m α(2 m) τ(m). (48) Property (47) is a quite remarable fractal property; from the origin in every irection we see either the zero sequence, or α, at ifferent scales. We claim that h[a] κ[a]i (49) follows from (45), (48), an (27). From (49) an (46) we obtain b[a] κ[a]t[2a]. (50) 13

14 We now prove that 1 +x 2 ivies q[a] when a is o. Noting that b[a](2a + ) b[a](2a ) an, when 2 a, b[a](2) 0, we therefore have q[a](x) 2a n1 b[a](2n 1)x 2n 2 a b[a](2a 2m + 1)x 2a 2m + b[a](2a + 2m 1)x 2a+2m 2 m1 a b[a](2a 2m + 1)x 2a 2m (1 + x 4m 2 ), m1 which vanishes at the points where x 2 1. Apart from the factor 1 + x 2 when a is o, the polynomial q[a] seems to be irreucible over Z an its zeros are in some sense close to the 2a th roots of 1 or, to the (a + 1) st roots of unity, see Figure 3. Figure 3: The roots of q[19] are epicte as boxes, the 38 th roots of 1 as points, an the 20 th roots of unity as crosses A perfect square Our last ientity generalises the faint fact that φ(1) 1. We have φ a (n) n 2. (51) a1 Proof [of (51)] For any lattice point (i, j) in the square [1, n] [1, n] the prouct i j mo n is congruent to some a in the range [1, n]. Acnowlegment This research has been fune by the Australian Research Council through the Centre of Excellence for Mathematics an Statistics of Complex Systems. 14

15 References [1] D.R. Anerson an T.M. Apostol, The evaluation of Ramanujan s sum an generalizations, Due Math. J. 20 (1953) [2] T.M. Apostol, Introuction to Analytic Number Theory, (1976) Springer- Verlag, New Yor. [3] T.M. Apostol, Arithmetical properties of generalize Ramanujan sums, Pacific J. Math. 41 (1972) [4] K.A. Broughan (2001), The gc-sum function, Journal of Integer Sequences 4, Art [5] L.E. Dicson, History of the theory of numbers, Carnegie Inst., Washington, D.C., 1919; reprinte by Chelsea, New Yor, [6] G.L. Dirichlet, Über ie Bestimmung er mittleren Werte in er Zahlentheorie, Abh. Aa. Wiss. Berlin, 1849; Also in Were, vol. 2, 1897, [7] O. Höler, Zur Theorie er Kreisteilungsgleichung K m (x) 0, Prace Matematyczno Fizyczne 43 (1936) [8] J.C. Kluyver, Some formulae concerning the integers less than n an prime to n, in: KNAW, Proceeings 9 I, Amsteram, (1906) [9] J. Liouville, Sur quelques séries et prouits infinis, J. Math. Pures Appl. 2 (1857) [10] L. Tóth an P. Hauanen, The iscrete Fourier transform of r-even functions, Acta Univ. Sapientiae Math. 3 (2011) [11] S. Ramanujan, On Certain Trigonometric Sums an their Applications in the Theory of Numbers, Transactions of the Cambrige Philosophical Society 22 (1918) Also in: Collecte papers of Srinivasa Ramanujan, E. G.H. Hary et al, Chelsea Publ. Comp., New Yor (1962) [12] J. Sanor an B. Crstici, Hanboo of Number Theory II, (2004) Kluwer Aca. Publ., Dorrecht. [13] R. Daublebsy Von Sternec, Ein Analogon zur aitiven Zahlentheorie, Sitzber, Aa. Wiss. Wien, Math. Naturw. Klasse, vol. I ll (Abt. IIa) (1902), [14] H. S. Wilf, Generatingfunctionology, Acaemic Press, 2n eition, [15] W. Schramm, The Fourier transform of functions of the greatest common ivisor, Integers 8 (2008) A50. 15