OFFSET
1
LINKS
FORMULA
a(n) = [A355934(n) == 1], where [ ] is the Iverson bracket.
MATHEMATICA
f[p_, e_] := ((q = NextPrime[p])^(e + 1) - 1)/(q - 1); a[1] = 1; a[n_] := If[Divisible[Times @@ f @@@ FactorInteger[n], DivisorSigma[1, n]], 1, 0]; Array[a, 100] (* Amiram Eldar, Jul 22 2022 *)
PROG
CROSSREFS
KEYWORD
nonn
AUTHOR
Antti Karttunen, Jul 22 2022
STATUS
approved