OFFSET
1,1
EXAMPLE
7 is a term because 7 is not in A000408 and 7^2 = 49 = 2^2 + 3^2 + 6^2.
PROG
(PARI) isA000408(n) = my(a, b) ; a=1 ; while(a^2+1<n, b=1 ; while(b<=a && a^2+b^2<n, if(issquare(n-a^2-b^2), return(1) ) ; b++ ; ) ; a++ ; ) ; return(0) ;
lista(nn) = for(n=1, nn, if(!isA000408(n) && isA000408(n^2), print1(n, ", ")));
CROSSREFS
KEYWORD
nonn,easy
AUTHOR
Altug Alkan, Jun 16 2016
STATUS
approved