%I #16 May 19 2019 11:34:11
%S 0,1,2,4,5,8,9,10,13,18,25,26,29,34,50,64,65,68,73,89,128,169,170,173,
%T 178,194,233,338,441,442,445,450,466,505,610,882,1156,1157,1160,1165,
%U 1181,1220,1325,1597,2312,3025,3026,3029,3034,3050,3089,3194,3466
%N Sums of two squares of Fibonacci numbers.
%H T. D. Noe, <a href="/A045702/b045702.txt">Table of n, a(n) for n = 1..1000</a>
%H <a href="/index/Su#ssq">Index entries for sequences related to sums of squares</a>
%F a(4)=5 because 5 = 1^2 + 2^2;
%F a(8)=13 because 13 = 2^2 + 3^2.
%t Take[Union[Total/@Tuples[Fibonacci[Range[0,15]]^2,{2}]],60] (* _Harvey P. Dale_, Jan 17 2011 *)
%o (PARI) list(lim)=my(sq=sqrtint(lim\=1),v=List(),f=List([0,1]),t); while((t=f[#f]+f[#f-1])<=sq, listput(f,t)); f=apply(sqr,f); for(i=1,#f, for(j=1,i, t=f[i]+f[j]; if(t>lim, break); listput(v,t))); Set(v) \\ _Charles R Greathouse IV_, Jun 27 2017
%K nonn,easy,nice
%O 1,3
%A _Felice Russo_
|