login
n appears once if n is the sum of 2 nonzero squares in 1 way, twice if n is the sum of 2 squares in 2 ways, 3 times if n is the sum of 2 squares 3 ways etc.
0

%I #28 Dec 01 2013 13:08:07

%S 2,5,5,8,10,10,13,13,17,17,18,20,20,25,25,26,26,29,29,32,34,34,37,37,

%T 40,40,41,41,45,45,50,50,50,52,52,53,53,58,58,61,61,65,65,65,65,68,68,

%U 72,73,73,74,74,80,80,82,82,85,85,85,85,89,89,90,90,97,97,98,100,100

%N n appears once if n is the sum of 2 nonzero squares in 1 way, twice if n is the sum of 2 squares in 2 ways, 3 times if n is the sum of 2 squares 3 ways etc.

%H <a href="/index/Su#ssq">Index entries for sequences related to sums of squares</a>

%F Does limit n -> infinity a(n)/n exist?

%F Answer: Yes, and it is equal to 4/Pi. - _Charles R Greathouse IV_, Dec 01 2013

%o (PARI) lista(nn) = for (n=1, nn, for (i=1, n-1, if (issquare(i) && issquare(n-i), print1(n, ", ")))) \\ _Michel Marcus_, Nov 30 2013

%o (PARI) list(lim)=my(v=List());for(n=2,lim,for(k=1,sqrtint(n-1),if(issquare(n-k^2),listput(v,n))));Vec(v) \\ _Charles R Greathouse IV_, Dec 01 2013

%Y Cf. A000404, A001105, A025285, A025302.

%K nonn

%O 1,1

%A _Benoit Cloitre_, Nov 09 2002

%E Typo in definition corrected by _Michel Marcus_, Dec 01 2013