%I #26 Jan 09 2023 12:33:02
%S 1,1,0,0,0,0,0,0,0,1,1,0,0,0,0,0,0,0,0,0,0,0,0,0,0,1,1,0,0,0,0,0,0,0,
%T 1,1,0,0,0,0,0,0,0,0,0,0,0,0,0,1,1,0,0,0,0,0,0,0,1,1,0,0,0,0,0,0,0,0,
%U 0,0,0,0,0,0,1,1,0,0,0,0,0,1,1,1,1,0,0,0,0,0,1,1,0,0,0,0,0,0,0,0,0,0,0,0,0
%N Number of partitions of n into distinct odd squares.
%C A167701 and A167702 give record values and where they occur: A167701(n)=a(A167702(n)) and a(m) < A167701(n) for m < A167702(n);
%C a(A167703(n)) = 0.
%H Reinhard Zumkeller, <a href="/A167700/b167700.txt">Table of n, a(n) for n = 0..10000</a>
%H Vaclav Kotesovec, <a href="/A167700/a167700.jpg">Graph - The asymptotic ratio</a>
%H <a href="/index/Su#ssq">Index entries for sequences related to sums of squares</a>.
%F a(n) = f(n,1,8) with f(x,y,z) = if x<y then 0^x else f(x-y,y+z,z+8) + f(x,y+z,z+8).
%F G.f.: Product_{k>=0} (1 + x^((2*k+1)^2)). - _Ilya Gutkovskiy_, Jan 11 2017
%F a(n) ~ exp(3 * 2^(-7/3) * Pi^(1/3) * (sqrt(2)-1)^(2/3) * Zeta(3/2)^(2/3) * n^(1/3)) * (sqrt(2)-1)^(1/3) * Zeta(3/2)^(1/3) / (2^(7/6) * sqrt(3) * Pi^(1/3) * n^(5/6)). - _Vaclav Kotesovec_, Sep 18 2017
%e a(50) = #{49+1} = 1;
%e a(130) = #{121+9, 81+49} = 2.
%t nmax = 100; CoefficientList[Series[Product[1 + x^((2*k-1)^2), {k, 1, Floor[Sqrt[nmax]/2] + 1}], {x, 0, nmax}], x] (* _Vaclav Kotesovec_, Sep 18 2017 *)
%o (Haskell)
%o a167700 = p a016754_list where
%o p _ 0 = 1
%o p (q:qs) m = if m < q then 0 else p qs (m - q) + p qs m
%o -- _Reinhard Zumkeller_, Mar 15 2014
%Y Cf. A033461, A016754, A167661, A000700, A111900.
%K nonn,look
%O 0,131
%A _Reinhard Zumkeller_, Nov 09 2009