login
A364252
a(n) = 1 if n has no other prime factors than 2 and/or Mersenne primes, otherwise 0.
3
1, 1, 1, 1, 0, 1, 1, 1, 1, 0, 0, 1, 0, 1, 0, 1, 0, 1, 0, 0, 1, 0, 0, 1, 0, 0, 1, 1, 0, 0, 1, 1, 0, 0, 0, 1, 0, 0, 0, 0, 0, 1, 0, 0, 0, 0, 0, 1, 1, 0, 0, 0, 0, 1, 0, 1, 0, 0, 0, 0, 0, 1, 1, 1, 0, 0, 0, 0, 0, 0, 0, 1, 0, 0, 0, 0, 0, 0, 0, 0, 1, 0, 0, 1, 0, 0, 0, 0, 0, 0, 0, 0, 1, 0, 0, 1, 0, 1, 0, 0, 0, 0, 0, 0, 0, 0, 0, 1, 0, 0, 0, 1, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 1, 0, 1, 1, 1
OFFSET
1
FORMULA
Fully multiplicative with a(p^e) = [p==2] + A036987(p), where [ ] is the Iverson bracket.
a(n) = [A336467(n) == 1].
a(n) >= A336923(n) >= A364251(n).
PROG
(PARI)
A209229(n) = (n && !bitand(n, n-1));
A364252(n) = { my(f=factor(n)); prod(k=1, #f~, (2==f[k, 1] || A209229(f[k, 1]+1))); };
CROSSREFS
Characteristic function of A219174.
Sequence in context: A257196 A176137 A290808 * A385888 A190239 A120529
KEYWORD
nonn,mult
AUTHOR
Antti Karttunen, Jul 16 2023
STATUS
approved