OFFSET
1,2
COMMENTS
{1,2,4,5,7,8,10,13,15,16,20,23,25,28,31,32,..}->{1,4,7,15,31,..}
LINKS
Chai Wah Wu, Table of n, a(n) for n = 1..10000
MATHEMATICA
q=66; q2=q^2+2; lst={}; Do[Do[Do[z=a^2+b^2+c^2; If[z<=q2, AppendTo[lst, z]], {c, b, 1, -1}], {b, a, 1, -1}], {a, q}]; lst; u=Union@lst; a=Complement[Range[q^2], u]; lst={}; Do[If[a[[n+1]]-a[[n]]==1, AppendTo[lst, a[[n]]]], {n, Length[a]-1}]; lst
CROSSREFS
KEYWORD
nonn
AUTHOR
Vladimir Joseph Stephan Orlovsky, May 30 2010
STATUS
approved