This is a front-end for the Online Encyclopedia of Integer Sequences, made by Christian Perfect. The idea is to provide OEIS entries in non-ancient HTML, and then to think about how they're presented visually. The source code is on GitHub.
%I A147599 #19 Nov 22 2020 09:00:50 %S A147599 1,0,0,1,0,0,0,1,0,0,1,1,0,0,1,1,0,0,2,1,0,1,2,1,0,1,3,1,0,2,3,1,0,3, %T A147599 4,1,1,4,4,1,1,5,5,1,2,7,5,1,3,8,6,1,5,10,6,2,6,12,7,2,9,14,7,3,11,16, %U A147599 8,4,15,19,8,6,18,21,9,8,23,24,10,11,27,27,11,14,34,30,12,19,39,33,14,24,47 %N A147599 Expansion of Product_{i>=1} (1+x^(4*i-1)). %C A147599 Number of partitions into distinct parts 4*k+3. %C A147599 Convolution of A147599 and A169975 is A000700. - _Vaclav Kotesovec_, Jan 18 2017 %H A147599 Vaclav Kotesovec, <a href="/A147599/b147599.txt">Table of n, a(n) for n = 0..10000</a> %F A147599 G.f. sum(n>=0, x^(2*n^2+n) / prod(k=1,n, 1-x^(4*k))) - _Joerg Arndt_, Mar 10 2011. %F A147599 a(n) ~ exp(sqrt(n/3)*Pi/2) / (4*6^(1/4)*n^(3/4)) * (1 - (3*sqrt(3)/(4*Pi) + Pi/(192*sqrt(3))) / sqrt(n)). - _Vaclav Kotesovec_, Jan 18 2017 %t A147599 nmax = 200; CoefficientList[Series[Product[(1 + x^(4*k - 1)), {k, 1, nmax}], {x, 0, nmax}], x] (* _Vaclav Kotesovec_, Jan 18 2017 *) %t A147599 nmax = 200; poly = ConstantArray[0, nmax + 1]; poly[[1]] = 1; poly[[2]] = 0; Do[If[Mod[k, 4] == 3, Do[poly[[j + 1]] += poly[[j - k + 1]], {j, nmax, k, -1}]; ], {k, 2, nmax}]; poly (* _Vaclav Kotesovec_, Jan 18 2017 *) %Y A147599 Cf. A000700, A169975, A170956-A170965. %Y A147599 Cf. A262928, A281243, A281244, A281245. %K A147599 nonn %O A147599 0,19 %A A147599 _N. J. A. Sloane_, Aug 29 2010