A070037 Nonprime numbers k such that sigma(k) == k+1 (mod phi(k)).
1, 4, 15, 900, 903, 28611063
Offset: 1
Examples
Below 30000000 only 5 composite numbers were found: C = {4,15,900,902,28611063}, Chowla(C) = {2,8,1920,504,17600976}, phi(C) = {2,8,240,504,17600976}, quotient = {1,1,8,1,1}.
Programs
-
Mathematica
s2[x_] := DivisorSigma[1, x]-x-1 e0[x_] := EulerPhi[x] Do[s=s2[n]/e0[n]; If[IntegerQ[s]&&!PrimeQ[n], Print[{n, s2[n], e0[n], s}]], {n, 1, 1000000}]
-
PARI
is(n)=!isprime(n) && Mod(sigma(n),eulerphi(n))==n+1 \\ Charles R Greathouse IV, Jun 06 2013
Comments