login
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