OFFSET
1,2
LINKS
Robert Israel, Table of n, a(n) for n = 1..10000
MAPLE
N:= 1000: # for terms <= N
sort(convert({seq(i^2, i=1..floor(sqrt(N))), seq(i*(i+1)/2, i=1 .. (floor(sqrt(1+8*N))-1)/2),
seq(i*(3*i-1)/2, i=1..(floor(sqrt(1+24*N))+1)/6)}, list)); # Robert Israel, Nov 11 2025
PROG
(PARI) isok(m) = ispolygonal(m, 3) || ispolygonal(m, 4) || ispolygonal(m, 5);
CROSSREFS
KEYWORD
nonn
AUTHOR
Michel Marcus, Mar 13 2021
STATUS
approved
