login

Year-end appeal: Please make a donation to the OEIS Foundation to support ongoing development and maintenance of the OEIS. We are now in our 61st year, we have over 378,000 sequences, and we’ve reached 11,000 citations (which often say “discovered thanks to the OEIS”).

A159458
Numbers m such that m^2 divides A159068(m), where A159068(m) = Sum_{k=1..m} binomial(m,k) * gcd(m,k).
1
1, 2, 3, 11, 33, 309, 665, 1461, 2323, 6789
OFFSET
1,2
COMMENTS
m divides A159068(m) for all positive integers m.
PROG
(PARI) is(m) = sum(k=1, m, binomial(m, k)*gcd(k, m))%m^2 == 0; \\ Jinyuan Wang, Aug 10 2021
CROSSREFS
Sequence in context: A300771 A296005 A374310 * A305846 A057838 A219497
KEYWORD
nonn,more
AUTHOR
Leroy Quet, Apr 12 2009
EXTENSIONS
4 more terms from R. J. Mathar, Apr 16 2009
a(10)=6789 from Ray Chandler, Jun 19 2009
STATUS
approved