A015762 Numbers n such that phi(n) | sigma_4(n).

1, 2, 3, 6, 249, 498

Offset: 1

Author

Robert G. Wilson v

Keywords

• nonn

Comments

sigma_4(n) is the sum of the 4th powers of the divisors of n (A001159).

sigma_{8j+4}(x)/phi(x) is an integer for j=0..500, x=1,2,3,6,249,498, and this is conjectured to hold for possible larger terms of A015762 and all j. Compare with comments to A015759, A091285, A015770. - Labos Elemer, May 27 2004

For any odd n in this sequence, 2n is also in the sequence, since phi(2n) = phi(n) and sigma_4(2n) = 17 sigma_4(n). More generally, if gcd(m,n) = 1 and m and n both are in this sequence, then mn is also in the sequence. No odd prime > 3 can be in the sequence, since if p = 2r + 1, then sigma_4(p) = 8r(2r^3 + 4r^2 + 3r + 1) + 2 is divisible by phi(p) = 2r only for r = 1. The term a(5) = 3*83 is the only odd semiprime term with a factor < 10^5. - M. F. Hasler, Aug 21 2017

a(7) > 3*10^11, if it exists. - Giovanni Resta, Aug 23 2017

Crossrefs

Cf. A020492, A015759, A015761, A015763, A015764, A015765, A015766, A015767, A015768, A015769, A015770, A015771, A015773, A015774, A094470.

Programs

Magma

[n: n in [1..1000]| DivisorSigma(4, n) mod EulerPhi(n) eq 0]; // Vincenzo Librandi, Aug 22 2017

Mathematica

Select[Range[500],Divisible[DivisorSigma[4,#],EulerPhi[#]]&] (* Harvey P. Dale, Dec 16 2012 *)

PARI

select( is(n)=sigma(n,4)%eulerphi(n)==0, [1..10^4])  \\ M. F. Hasler, Aug 21 2017