OFFSET
1,2
COMMENTS
Numbers k such that k | A000182(k).
All the powers of 2 are terms.
EXAMPLE
MATHEMATICA
Select[Range[1000], Divisible[((-4)^# - (-16)^#) * BernoulliB[2*#]/(2*#), #] &]
PROG
(PARI) is(n) = (((-4)^n - (-16)^n) * bernfrac(2*n) / (2*n)) % n == 0;
CROSSREFS
KEYWORD
nonn,more
AUTHOR
Amiram Eldar, May 17 2024
STATUS
approved