OFFSET
1,2
COMMENTS
Any terms for n > 25 exceed 10^10. - Lucas A. Brown, Dec 24 2022
a(25) <= 8147739600, a(26) <= 26771144400, a(27) <= 36082846800, a(28) <= 80313433200. - Jon E. Schoenfield, Dec 16 2022
LINKS
Lucas A. Brown, Python program.
Eric Weisstein's World of Mathematics, Square Pyramidal Number.
EXAMPLE
a(3) = 30 because 30 has 3 square pyramidal divisors {1, 5, 30} and this is the smallest such number.
PROG
(PARI) issqpyr(n) = my(m = sqrtnint(3*n, 3)); n==m*(m+1)*(2*m+1)/6; \\ A253903
a(n) = my(k=1); while (sumdiv(k, d, issqpyr(d)) != n, k++); k; \\ Michel Marcus, Nov 21 2022
CROSSREFS
KEYWORD
nonn,more
AUTHOR
Ilya Gutkovskiy, Nov 21 2022
EXTENSIONS
a(15) from Michel Marcus, Nov 21 2022
a(16)-a(20) from Jinyuan Wang, Nov 28 2022
a(21)-a(22) from Lucas A. Brown, Dec 14 2022
a(23)-a(24) from Lucas A. Brown, Dec 18 2022
a(25) from Lucas A. Brown, Dec 22 2022
STATUS
approved