OFFSET
1,1
LINKS
Robert Israel, Table of n, a(n) for n = 1..10000
EXAMPLE
a(1) = 2 = 1+1.
25 is a term because 25 = 16 + 9 = 10 + 15 = 4 + 21, and 16, 9, 4 are squares, 10, 15, 21 are triangular numbers.
MAPLE
N:= 1000: # for terms <= N
SS:= {seq(seq(i^2+j^2, i=1..min(j, floor(sqrt(N-j^2)))), j=1..floor(sqrt(N)))}:
TS:= {seq(seq(i*(i+1)/2 + j^2, j = 1..floor(sqrt(N-i*(i+1)/2))), i=1..floor((sqrt(1+8*N)-1)/2))}:
TT:= {seq(seq(i*(i+1)/2 + j*(j+1)/2, i = 1 .. min(j, floor((sqrt(1+8*(N-j*(j+1)/2))-1)/2))), j=1..floor((sqrt(1+8*N)-1)/2))}:
sort(convert(SS intersect TS intersect TT, list)); # Robert Israel, Mar 19 2026
CROSSREFS
KEYWORD
nonn
AUTHOR
Alex Ratushnyak, Nov 03 2015
STATUS
approved
