OFFSET
1,1
LINKS
Amiram Eldar, Table of n, a(n) for n = 1..10000
FORMULA
sigma_13(a(n)) mod phi(a(n)) = 0; sigma(a(n)) mod phi(a(n)) <> 0.
MATHEMATICA
q[n_] := Sign[Mod[DivisorSigma[{1, 13}, n], EulerPhi[n]]] == {1, 0}; Select[Range[24000], q] (* Amiram Eldar, Apr 11 2024 *)
PROG
(PARI) is(n) = {my(f = factor(n), phi = eulerphi(f)); (sigma(f) % phi) && !(sigma(f, 13) % phi); } \\ Amiram Eldar, Apr 11 2024
CROSSREFS
KEYWORD
nonn
AUTHOR
Labos Elemer, Dec 05 2002
STATUS
approved