OFFSET
1,2
MAPLE
# return true if 'candid' is allowed (not a sum of squares)
A026469aux := proc(a, candid) local i, j; for j from 1 to nops(a) do for i from 1 to j do if (op(i, a))^2+(op(j, a))^2 = candid then RETURN(false); fi; od; od; RETURN(true); end:
A026468 := proc(nmax) local a, candidat; a := [1, 2]; while nops(a) < nmax do candidat := op(nops(a), a)+1; while A026469aux(a, candidat) = false do candidat := candidat+1; od; a := [op(a), candidat]; od: RETURN(a); end: A026468(60); # R. J. Mathar, Nov 01 2006
CROSSREFS
KEYWORD
nonn
AUTHOR
EXTENSIONS
Definition corrected by Ralf Stephan, Nov 01 2006
STATUS
approved