OFFSET
1
LINKS
FORMULA
For all n >= 1, a(n) <= A010055(n).
Sum_{k=1..n} a(k) ~ n/log n. - Bill McEachen and Charles R Greathouse IV, Sep 11 2022
MATHEMATICA
a[n_] := If[!CompositeQ[DivisorSigma[0, n]], 1, 0]; Array[a, 100] (* Amiram Eldar, Jul 21 2022 *)
PROG
(PARI) A355937(n) = ((1==n)||isprime(isprimepower(n)+1));
CROSSREFS
KEYWORD
nonn
AUTHOR
Antti Karttunen, Jul 21 2022
STATUS
approved
