More fun in 2015, Part 33

S = gcd ( k 1 , 2015 ) S=\sum\gcd(k-1,2015)

If the sum above is taken over all integers k k with 1 k 2015 1\leq k\leq 2015 and gcd ( k , 2015 ) = 1 \gcd(k,2015)=1 , find S S .


The answer is 11520.

This section requires Javascript.
You are seeing this because something didn't load right. We suggest you, (a) try refreshing the page, (b) enabling javascript if it is disabled on your browser and, finally, (c) loading the non-javascript version of this page . We're sorry about the hassle.

2 solutions

Set N N be a square-free number (and we will take N = 2015 = 5 13 31 N = 2015 = 5\cdot 13\cdot 31 ), and d d be a divisor of N N .

Lemma: gcd ( k 1 , N ) k is a multiple of d = gcd ( , N / d ) 0 < N / d . \langle \text{gcd}(k-1,N) | k\ \text{is a multiple of}\ d\rangle = \langle \text{gcd}(\ell,N/d) | 0 \leq \ell < N/d\rangle. (This is an equality of "bags", i.e. the lists of numbers are identical except for the order.)

Proof : Apply the Chinese Remainder Theorem to the system of equations { x 1 mod d , x mod N / d . \left\{\begin{array}{ll} x \equiv -1\ \text{mod}\ d, \\ x \equiv \ell\ \text{mod}\ N/d. \end{array}\right. to be guaranteed a unique solution 0 x = k 1 < N 0 \leq x = k-1 < N for every value of \ell . (It is essential that N N is square-free, so that d d and N / d N/d are coprime.) This defines a bijection between the values of k k and \ell . Moreover, gcd ( k 1 , N ) = gcd ( x , N ) = gcd ( x , N / d ) = gcd ( , N / d ) . \gcd(k-1,N) = \gcd(x,N) = \gcd(x,N/d) = \gcd(\ell,N/d).\ \ \ \square

Now define Σ ( n ) : = = 0 n 1 gcd ( , n ) \Sigma(n) := \sum_{\ell = 0}^{n-1} \text{gcd}(\ell, n) for square-free integers n n .

Lemma: Σ \Sigma is a multiplicative function and for prime numbers p p , Σ ( p ) = 2 p 1. \Sigma(p) = 2p-1. (I leave the proof to you.)

Now to solve the problem we observe that the set of k k values over which we sum is { multiples of 1 } { multiples of 5 } { multiples of 13 } { multiples of 31 } + { multiples of 65 } + { multiples of 155 } + { multiples of 403 } { multiples of 2015 } \{\text{multiples of 1}\} - \{\text{multiples of 5}\} - \{\text{multiples of 13}\} - \{\text{multiples of 31}\} \\ + \{\text{multiples of 65}\} + \{\text{multiples of 155}\} + \{\text{multiples of 403}\} - \{\text{multiples of 2015}\} (E.g. after subtracting multiples of 5 and multiples of 13 we have subtracted all multiples of 65 twice, so we need to add them back in once, and so on. This is a common counting strategy.)

We obtain the desired sum by adding and subtraction various sums along the same lines; and note that k is multiple of d gcd ( k 1 , 2015 ) = = 0 2015 / d 1 gcd ( , 2015 / d ) = Σ ( 2015 / d ) . \sum_{k\ \text{is multiple of}\ d} \text{gcd}(k-1,2015) = \sum_{\ell = 0}^{2015/d-1} \text{gcd}(\ell, 2015/d) = \Sigma(2015/d). Thus we get as result Σ ( 2015 ) Σ ( 403 ) Σ ( 155 ) Σ ( 65 ) + Σ ( 31 ) + Σ ( 13 ) + Σ ( 5 ) Σ ( 1 ) = 9 25 61 25 61 9 61 9 25 + 61 + 25 + 9 1 = 11520 . \Sigma(2015) - \Sigma(403) - \Sigma(155) - \Sigma(65) + \Sigma(31) + \Sigma(13) + \Sigma(5) - \Sigma(1) \\ = 9\cdot 25\cdot 61 - 25\cdot 61 - 9\cdot 61 - 9\cdot 25 + 61 + 25 + 9 - 1 = \boxed{11520}.

Very nice! You are proving Menon's identity in the case when N N is square-free, which makes the proof more transparent (compare with the general proof ) . Your sum in the last two lines is simply ( 9 1 ) ( 25 1 ) ( 61 1 ) = p N ( Σ ( p ) 1 ) = p N ( 2 p 2 ) = 2 m p N ( p 1 ) = τ ( N ) ϕ ( N ) (9-1)(25-1)(61-1)=\prod_{p|N}(\Sigma(p)-1)=\prod_{p|N}(2p-2)=2^{m}\prod_{p|N}(p-1)=\tau(N)\phi(N) as in Menon's identity, where m m is the number of prime factors of the square free integer N N .

Stated succinctly, the sum we seek is the Dirichlet convolution of the Mobius function with the gcd-sum function, d N μ ( N / d ) Σ ( d ) \sum_{d|N}\mu(N/d)\Sigma(d) .

Otto Bretscher - 5 years, 6 months ago
Otto Bretscher
Dec 5, 2015

By Menon's elegant identity, the answer is τ ( 2015 ) ϕ ( 2015 ) \tau(2015)\phi(2015) , where τ \tau denotes the number of positive integer factors. Now 2015 = 5 13 31 2015=5*13*31 so that τ ( 2015 ) ϕ ( 2015 ) = 2 3 1440 = 11520 \tau(2015)\phi(2015)=2^3*1440=11520 .

Menon Identity would be a nice way.

But is there any way to solve it using elementary way? I tried it but failed.

Dev Sharma - 5 years, 6 months ago

Log in to reply

Now that school is out, I had a bit of time to think about this. For a number like 2015, with just 3 (distinct) prime factor, it's not hard to find this sum by inspection. All we need is some careful counting and the Euclidean algorithm (EA) that allows us to write 1 as a linear combination of any two co-prime integers.

Let me try to work this out in an "elementary way".

It might be slightly easier to find C = gcd ( k 1 , 2015 ) C=\sum\gcd(k-1,2015) where we sum over all k k with gcd ( k , 2015 ) 1 \gcd(k,2015)\neq 1 . Then we have S = 13725 C S=13725-C where 13725 is the gcd-sum of 2015.

Let p 1 = 5 , p 2 = 13 , p 3 = 31 p_1=5,p_2=13,p_3=31 be the prime factors of 2015 2015 . All integers we consider will be between 1 and 2015, inclusive.

How often does the summand p 1 p 2 p_1p_2 appear in C C ? By the EA, there is exactly one multiple j p 3 jp_3 such that j p 3 1 jp_3-1 is divisible by p 1 p 2 p_1p_2 . Thus the summand p 1 p 2 p_1p_2 appears exactly once in C C ; likewise for p 1 p 3 p_1p_3 and p 2 p 3 p_2p_3 .

How often does the summand p 1 p_1 appear in C C ? There are p 2 1 p_2-1 multiples j p 3 jp_3 such that gcd ( j p 3 1 , N ) = p 1 \gcd(jp_3-1,N)=p_1 , and, likewise, there are p 3 1 p_3-1 multiples j p 2 jp_2 such that gcd ( j p 2 1 , N ) = p 1 \gcd(jp_2-1,N)=p_1 . With one double count for the multiples of p 2 p 3 p_2p_3 , the summand p 1 p_1 appears p 2 + p 3 3 p_2+p_3-3 times in C C . Likewise for p 2 p_2 and p 3 p_3 , mutatis mutandis.

The remaining summands of S S are all 1; we count 483 of them.

Now we are ready for the final tally: C = 3 ( 5 13 + 5 31 + 13 31 ) 3 ( 5 + 13 + 31 ) + 483 1 = 2205 C=3(5*13+5*31+13*31)-3(5+13+31)+483*1=2205 and S = 13725 2205 = 11520 S=13725-2205=\boxed{11520}

Otto Bretscher - 5 years, 6 months ago

0 pending reports

×

Problem Loading...

Note Loading...

Set Loading...