cp's OEIS Frontend

A333174 a(n) = Sum_{k=0..n} r_4(k^2 + 1), where r_4(k) is the number of ways of writing k as a sum of 4 squares (A000118).

%I A333174 #8 Mar 10 2020 02:48:02
%S A333174 8,32,80,224,368,704,1008,1752,2424,3432,4248,5736,7176,9768,11352,
%T A333174 14088,16152,20472,23944,28312,31528,37576,42280,50056,54680,62216,
%U A333174 67640,78296,85880,96008,103784,116552,126968,142808,152888,167624,178008,197880,212616,230904
%N A333174 a(n) = Sum_{k=0..n} r_4(k^2 + 1), where r_4(k) is the number of ways of writing k as a sum of 4 squares (A000118).
%H A333174 Amiram Eldar, <a href="/A333174/b333174.txt">Table of n, a(n) for n = 0..10000</a>
%H A333174 R. Sitaramachandrarao and P. V. Krishnaiah, <a href="https://doi.org/10.1016/0022-314X(86)90086-7">On the sums Sigma_{n<=x} A(f(n)) and Sigma_{p<=x} A(f(p))</a>, Journal of Number Theory, Vol. 23, No. 2 (1986), pp. 149-168.
%F A333174 a(n) ~ (40*G/Pi^2) * n^3, where G is Catalan's constant (A006752).
%e A333174 a(0) = r_4(0^2 + 1) = r_4(1) = A000118(1) = 8.
%e A333174 a(1) = r_4(0^2 + 1) + r_4(1^1 + 1) = r_4(1) + r_4(2) = A000118(1) + A000118(2) = 8 + 24 = 32.
%t A333174 Accumulate @ Table[SquaresR[4, k^2 + 1], {k, 0, 100}]
%Y A333174 Partial sums of A333173.
%Y A333174 Cf. A000118, A002522, A006752.
%K A333174 nonn
%O A333174 0,1
%A A333174 _Amiram Eldar_, Mar 09 2020