Treat yourself with a hard number theory problem

n n is a positive integer, and d ( n ) d(n) is defined as the number of positive factors of n . n.
For example, 12 12 has six positive factors 1 , 2 , 3 , 4 , 6 , 1,2,3,4,6, and 12 12 in all, so d ( 12 ) = 6. d(12)=6.

Find the sum of all the positive integers n n satisfying ( d ( n ) ) 3 = 4 n \big(d(n)\big)^3 = 4n .


The answer is 2130.

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.

3 solutions

Haosen Chen
Feb 15, 2018

In this solution, I will use a part of Bernoulli's inequality: ( 1 + x ) r 1 + r x \displaystyle (1+x)^{r}\geq 1+rx for every integer r 0 \displaystyle r ≥ 0 and every real number x 1 \displaystyle x ≥ -1 .

\bullet https://en.wikipedia.org/wiki/Bernoulli%27s_inequality …………………………………………………………………………………………………………………………………………………………………………………………………

1. Prime factorization: n = i = 1 k p i α i \displaystyle n=\prod_{i=1}^{k} p_i^{\alpha_i} ,here prime factors p 1 < p 2 < < p k \ p_{1}<p_{2}<…<p_{k} . 4 n 4n is a perfect cube \Rightarrow p 1 = 2 p_{1}=2 , α 1 = 3 β 1 + 1 \alpha_{1} = 3\beta_{1} +1 , α 2 = 3 β 2 \alpha_{2} = 3\beta_{2} ,…, α k = 3 β k \alpha_{k} = 3\beta_{k} , here β i N , 1 i k \beta_{i}\in N, \forall 1\le i\le k . From d ( n ) 3 = 4 n d(n)^{3} =4n ,we get d ( n ) = ( 3 β 1 + 2 ) ( 3 β 2 + 1 ) ( 3 β k + 1 ) \displaystyle d(n)={(3\beta_{1}+2)(3\beta_{2}+1)…(3\beta_{k}+1)} = 2 β 1 + 1 i = 2 k p i β i {\displaystyle =2^{\beta_{1}+1}\displaystyle \prod_{i=2}^{k} p_{i}^{\beta_{i}}} , (#) , since we have the formula d ( n ) = i = 1 k ( α i + 1 ) \displaystyle d(n)=\prod_{i=1}^{k} (\alpha_i +1) .

2. Rewrite (#) : 3 β 1 + 2 2 β 1 + 1 \displaystyle \frac{3\beta_{1}+2}{2^{\beta_{1}+1}} = p 2 p 3 p k ( 3 β 2 + 1 ) ( 3 β 3 + 1 ) ( 3 β k + 1 ) \displaystyle =\frac{p_{2}p_{3}…p_{k}}{(3\beta_{2}+1)(3\beta_{3}+1)…(3\beta_{k}+1)} .

Because d ( n d(n ) is not a multiple of 3 3 ,we see p 2 5 p_{2}\ge 5 . So when i 2 i\ge 2 , holds p i ( 1 + 4 ) β i 1 + 4 β i \displaystyle p_{i}\ge (1+4)^{\beta_{i}}\ge 1+4\beta_{i} ( ) (\divideontimes) by Bernoulli's inequality.

Apply ( ) (\divideontimes) into (#) and get 3 β 1 + 2 2 β 1 + 1 1 \displaystyle \frac{3\beta_{1}+2}{2^{\beta_{1}+1}} \ge 1 β 1 2 \Rightarrow \beta_{1}\le 2

3. Two cases to be consider.

① When β 1 = 0 , 2 \beta_{1}=0,2 , substitute into (#) and so p 2 p 3 p k ( 3 β 2 + 1 ) ( 3 β 3 + 1 ) ( 3 β k + 1 ) = 1 \displaystyle \frac{p_{2}p_{3}…p_{k}}{(3\beta_{2}+1)(3\beta_{3}+1)…(3\beta_{k}+1)}=1 which implies that ( ) (\divideontimes) takes equal sign for every i 2 i\ge 2 i.e. β i = 0 , 2 i k \beta_{i}=0,\forall 2\le i\le k . Thus,in this case n = 2 n=2 or 2 7 2^7 .

② When β 1 = 1 \beta_{1}=1 , substitute into (#) and so 5 4 = p 2 p 3 p k ( 3 β 2 + 1 ) ( 3 β 3 + 1 ) ( 3 β k + 1 ) \displaystyle\frac{5}{4}= \frac{p_{2}p_{3}…p_{k}}{(3\beta_{2}+1)(3\beta_{3}+1)…(3\beta_{k}+1)} .

WLOG β i 1 , 2 i k \beta_{i}\ge 1, \forall 2\le i\le k ,since there must exist β j 1 ( 2 j k ) \beta_{j}\ge 1 (2\le j\le k) to reach the 5 4 \frac{5}{4} .

Assume β j 2 ( 2 j k ) \exists \beta_{j}\ge 2 (2\le j\le k) ,then 5 4 > \displaystyle\frac{5}{4}> p j β j ( 3 β j + 1 ) > \displaystyle\frac{p_{j}^{\beta_{j}}}{(3\beta_{j}+1)}> 1 + 4 β j 3 β j + 1 > 5 4 \displaystyle\frac{1+4\beta_{j}}{3\beta_{j}+1}>\frac{5}{4} ,which is a contradiction. Therefore, β i = 1 , 2 i k \beta_{i}=1,\forall 2\le i\le k .

Assume p i 7 ( 2 i k ) \exists p_{i} \ge 7 (2\le i\le k) ,then 5 4 = p 2 p 3 p k ( 3 β 2 + 1 ) ( 3 β 3 + 1 ) ( 3 β k + 1 ) \displaystyle\frac{5}{4}=\frac{p_{2}p_{3}…p_{k}}{(3\beta_{2}+1)(3\beta_{3}+1)…(3\beta_{k}+1)} \ge p i 3 + 1 \displaystyle\frac{p_{i}}{3+1} 7 4 \displaystyle\ge \frac{7}{4} > > 5 4 \displaystyle \frac{5}{4} ,which is a contradiction.

Thus in this case, n = 2 4 5 3 \displaystyle n=2^{4}5^{3} . …………………………………………………………………………………………………………………………………………………………………………………………………

So, A = { 2 , 128 , 2000 } \displaystyle A=\{2,128,2000\} .Answer is 2130 \boxed{2130}

Wonderful problem buddy I liked it very much. Although failed to solve it.

Shreyansh Mukhopadhyay - 3 years, 3 months ago
Leonel Castillo
Feb 14, 2018

My solution will be based on the product formula for the divisor function. It is known that if n = p i a i n = \prod p_i ^{a_i} then d ( n ) = ( a i + 1 ) d(n) = \prod (a_i + 1) . This product representation turns the equation d ( n ) 3 = 4 n d(n)^3 = 4n into ( a i + 1 ) 3 = 4 p i a i \prod (a_i + 1)^3 = 4 \prod p_i ^{a_i} . Notice that if the left hand side is a perfect cube, that implies that the right hand side is a perfect cube too. But 4 = 2 2 4 = 2^2 is not a perfect cube. This means that one of the p i p_i must be 2, and it must be raised to a 3 k + 1 3k + 1 power so that 2 2 × 2 3 k + 1 = 2 3 k + 3 = ( 2 k + 1 ) 3 2^2 \times 2^{3k + 1} = 2^{3k + 3} = (2^{k+1})^3 . This means that p 1 = 2 p_1 = 2 and a 1 = 3 k + 1 a_1 = 3k + 1 and if we plug this into our equation we get ( 3 k + 2 ) 3 p i 2 ( a i + 1 ) 3 = 2 3 k + 3 p i 2 p i a i (3k + 2)^3 \prod_{p_i \neq 2} (a_i + 1)^3 = 2^{3k + 3} \prod_{p_i \neq 2} p_i ^{a_i} . A final deduction that will be useful later is that all of the powers of the other primes, the a i a_i , must also be multiples of 3 3 in order for the right hand side to be a perfect cube. We should now branch out into two cases.

Case 1: n n has no prime factors other than 2 Then we have to solve the equation ( 3 k + 2 ) 3 = 2 3 k + 3 3 k + 2 = 2 k + 1 (3k+2)^3 = 2^{3k + 3} \implies 3k + 2 = 2^{k+1} . This is a simple diophantine equation that can be solved simply by realizing that the left-hand side grows linearly, while the right-hand side grows exponentially and only the cases k = 0 , 1 , 2 k=0,1,2 must be checked. When checking we find that equality only holds when k = 0 , 2 k=0,2 which give us n = 2 2 , 2 3 × 2 + 1 n = 2^2, 2^{3\times 2 + 1} or, in simpler terms, 2 , 128 A 2, 128 \in A .

Case 2: n n has other prime factors In this case, a clever idea would be to find which prime factors it can be. The first technique that should come to mind is an argument by size, given that few other strategies can rule out infinitely many primes. Rewrite the equation as ( 3 k + 2 2 k + 1 ) 3 = p i 2 p i a i ( a i + 1 ) 3 \left( \frac{3k + 2}{2^{k+1}} \right)^3 = \prod_{p_i \neq 2} \frac{p_i^{a_i}}{(a_i + 1)^3} . Now suppose that n n has a prime factor bigger than 5 5 . Then p i 2 p i a i ( a i + 1 ) 3 7 3 4 3 3 3 4 3 = 2.260986328 \prod_{p_i \neq 2} \frac{p_i^{a_i}}{(a_i + 1)^3} \geq \frac{7^3}{4^3} \frac{3^3}{4^3} = 2.260986328 . But we can easily prove that ( 3 k + 2 2 k + 1 ) 3 < 2 \left( \frac{3k + 2}{2^{k+1}} \right)^3 < 2 for integer values of k k so equality will be impossible. This means that if n n has other prime factors, they can only be 3 3 or 5 5 .

I will now prove that n n can't have both 3 3 and 5 5 as prime factors. Suppose it does, then p i 2 p i a i ( a i + 1 ) 3 3 3 4 3 × 5 3 4 3 = 0.823974609 \prod_{p_i \neq 2} \frac{p_i^{a_i}}{(a_i + 1)^3} \geq \frac{3^3}{4^3} \times \frac{5^3}{4^3} = 0.823974609 . We know that ( 3 k + 2 2 k + 1 ) 3 \left( \frac{3k + 2}{2^{k+1}} \right)^3 is decreasing and for k > 2 k > 2 we have ( 3 k + 2 2 k + 1 ) 3 \left( \frac{3k + 2}{2^{k+1}} \right)^3 < 0.823974609). This means that if solutions exist, it must be that k = 1 , 2 k=1,2 . We can plug in those values and in both cases we will find a contradiction.

So we know that n n has only one other prime factor. Let's suppose it is 3 3 , then we have to solve ( 3 k + 2 2 k + 1 ) 3 = 3 a ( a + 1 ) 3 \left( \frac{3k + 2}{2^{k+1}} \right)^3 = \frac{3^a}{(a + 1)^3} . If a 6 a \geq 6 , a size argument proves equality is impossible. This means that a = 3 a=3 and then checking we see that no valid k k exists.

We now do the same for 5 5 and find that the equation ( 3 k + 2 2 k + 1 ) 3 = 5 a ( a + 1 ) 3 \left( \frac{3k + 2}{2^{k+1}} \right)^3 = \frac{5^a}{(a + 1)^3} has unique solution a = 3 , k = 1 a=3, k=1 which corresponds to n = 2 4 5 3 = 2000 n = 2^4 5^3 = 2000 .

Nice! Thank you,Leonel Castillo. And I will post my solution later(based on similar idea but less calculation).

Haosen Chen - 3 years, 3 months ago
Giorgos K.
Feb 21, 2018

just searched the first 10.000 integers and got the solution in under 1 sec with Mathematica

Select[Range@10000, DivisorSigma[0,#]^3==4#&]

{2, 128, 2000}

0 pending reports

×

Problem Loading...

Note Loading...

Set Loading...