OFFSET
1,3
LINKS
Amiram Eldar, Table of n, a(n) for n = 1..10000
EXAMPLE
5967 is a decagonal number whose sum of digits 5+9+6+7 = 27 is also decagonal.
MATHEMATICA
sod[n_] := Plus @@ IntegerDigits[n]; decQ[n_] := n == 0 || IntegerQ[(Sqrt[64*n + 36] + 6)/16]; Select[Table[n*(4*n - 3), {n, 0, 300}], decQ[sod[#]] &] (* Amiram Eldar, Feb 06 2021 *)
PROG
(PARI) isok(n) = ispolygonal(n, 10) && ispolygonal(sumdigits(n), 10); \\ Michel Marcus, Feb 06 2021
CROSSREFS
KEYWORD
base,nonn
AUTHOR
Luc Stevens (lms022(AT)yahoo.com), May 28 2006
EXTENSIONS
a(1)=0 and a missing term inserted by Amiram Eldar, Feb 06 2021
STATUS
approved