|
%I #19 May 29 2022 18:06:59
%S 1,1,1,1,2,1,2,1,2,2,2,1,3,3,2,2,3,2,3,2,3,3,4,2,4,3,3,3,4,3,4,4,4,4,
%T 4,2,5,5,4,3,6,4,5,4,5,5,5,3,6,6,5,5,6,4,5,5,5,6,8,4,6,6,7,5,7,5,7,7,
%U 6,6,6,5,8,8,6,5,10,6,8,6,7,7,8,5,8,10,7,8,8,6,8,7,9,9,11,5,8,10,7,7
%N Number of partitions of 2*n-1 into four squares.
%C This is a bisection of A002635.
%C While A002635 contains each positive integer infinitely often, here a number can appear only finitely many times.
%C By the Jacobi theorem, a(n) >= A000203(n)/48 >= (1+n)/48, which implies the previous comment. - _Robert Israel_, Nov 21 2017
%H Robert Israel, <a href="/A295218/b295218.txt">Table of n, a(n) for n = 1..10000</a>
%H Wikipedia, <a href="https://en.wikipedia.org/wiki/Jacobi%27s_four-square_theorem"> Jacobi four square theorem</a>
%p N:= 100: # to get a(1)...a(N)
%p V:= Array(0..2*N-1):
%p for a from 0 while 4*a^2 <= 2*N-1 do
%p for b from a while a^2 + 3*b^2 <= 2*N-1 do
%p for c from b while a^2 + b^2 + 2*c^2 <= 2*N-1 do
%p for d from c while a^2 + b^2 + c^2 + d^2 <= 2*N-1 do
%p t:= a^2 + b^2 + c^2 + d^2;
%p V[t]:= V[t]+1
%p od od od od:
%p seq(V[2*i-1],i=1..N); # _Robert Israel_, Nov 21 2017
%Y Cf. A002635, A295159.
%K nonn
%O 1,5
%A _Franklin T. Adams-Watters_, Nov 17 2017
|
