OFFSET
1
COMMENTS
For all multiperfect numbers, a(A007691(n)) = 1.
LINKS
FORMULA
a(n) = 1 if A353762(n) = 1, otherwise 0.
MATHEMATICA
f[p_, e_] := (p - 1)*p^(e - 1)*If[p == 2, 1, NextPrime[p, -1]^e]; s[1] = 1; s[n_] := Times @@ f @@@ FactorInteger[n]; a[n_] := Boole[Divisible[s[DivisorSigma[1, n]], s[n]] ]; Array[a, 120] (* Amiram Eldar, May 10 2022 *)
PROG
CROSSREFS
KEYWORD
nonn
AUTHOR
Antti Karttunen, May 10 2022
STATUS
approved