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:
A026469 := proc(nmax) local a, candidat ; a := [1] ; 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: A026469(60) ; # R. J. Mathar, Nov 01 2006
CROSSREFS
KEYWORD
nonn
AUTHOR
EXTENSIONS
Corrected by Franklin T. Adams-Watters, Oct 25 2006
Definition corrected by N. J. A. Sloane, Nov 01 2006
STATUS
approved