OFFSET
1,1
COMMENTS
For p prime, contains p^2 if and only if p == 1 (mod 6). - Robert Israel, May 14 2019
LINKS
Robert Israel, Table of n, a(n) for n = 1..10000
MAPLE
N:= 1000:
V:= Vector(N):
for d from 1 to N do
r:= numtheory:-sigma(d)/numtheory:-tau(d);
C:= [seq(i, i=d..N, d)];
V[C]:= V[C] +~ r
od:
select(t -> not(isprime(t)) and V[t]::integer, [$2..N]); # Robert Israel, May 14 2019
MATHEMATICA
q[n_] := CompositeQ[n] && IntegerQ @ DivisorSum[n, Divide @@ DivisorSigma[{1, 0}, #] &]; Select[Range[200], q] (* Amiram Eldar, Jun 08 2022 *)
CROSSREFS
KEYWORD
easy,nonn
AUTHOR
Benoit Cloitre, Apr 06 2002
STATUS
approved
