OFFSET
1,1
LINKS
Robert Israel, Table of n, a(n) for n = 1..10000 (replaces erroneous b-file from Harvey P. Dale)
MAPLE
pent:= n -> n*(3*n-1)/2:
filter:= proc(n) local x, y, S;
S:= [isolve(pent(x) + y^2 = n)];
ormap(t -> subs(t, x) > 0 and subs(t, y) > 0, S)
end proc:
select(filter, map(pent, [$1 .. 30])); # Robert Israel, Jul 21 2025
CROSSREFS
KEYWORD
nonn
AUTHOR
Jonathan Vos Post, Jan 21 2008
EXTENSIONS
Added more terms from b-file. - David Radcliffe, Jul 21 2025
STATUS
approved
