%I #48 Nov 22 2023 12:16:26
%S 1,1,1,2,3,4,5,7,10,13,17,22,28,36,46,58,73,91,114,141,173,213,261,
%T 318,387,469,567,683,821,984,1176,1403,1671,1984,2351,2781,3284,3869,
%U 4550,5343,6264,7330,8565,9993,11642,13543,15733,18252,21148,24471,28282,32646,37640,43348,49867,57302,65776,75426,86405,98882
%N Number of partitions of n into at most 1 copy of 1, 2 copies of 2, 3 copies of 3, ... .
%C Also number of partitions into non-pronic numbers (cannot be written as i*(i+1)).
%C Convolution of A024940 and A225044. - _Vaclav Kotesovec_, Jan 02 2017
%H Alois P. Heinz, <a href="/A052335/b052335.txt">Table of n, a(n) for n = 0..10000</a> (first 129 terms from Reinhard Zumkeller)
%H Mircea Merca and Emil Simion, <a href="https://doi.org/10.3390/sym15112067">n-Color Partitions into Distinct Parts as Sums over Partitions</a>, Symmetry (2023) Vol. 15, Iss. 11.
%F G.f.: Product_{i>=1} (1-x^(i*(i+1)))/(1-x^i).
%F G.f.: (1+x) * (1+x^2+x^4) * (1+x^3+x^6+x^9) * (1+x^4+x^8+x^12+x^16) * ... (g.f. above, expanded). - _Joerg Arndt_, Apr 01 2014
%F G.f.: Product_{n>=1} (1 - q^(n*(n+1))) / Product_{n>=1} (1 - q^n). - _Joerg Arndt_, Apr 01 2014
%F a(n) = p(n,1,1) with p(n,t,k) = if t<0 then 0 else if k<=n then p(n-k,t-1,k)+p(n,k+1,k+1) else 0^n. - _Reinhard Zumkeller_, Jan 20 2010
%F a(n) ~ exp(Pi*sqrt(2*n/3) - 3^(1/4) * Zeta(3/2) * n^(1/4) / 2^(3/4) - 3*Zeta(3/2)^2/(32*Pi)) / sqrt(2*n). - _Vaclav Kotesovec_, Jan 01 2017
%e a(5)=4 because we have [5], [4,1], [3,2] and [2,2,1] ([3,1,1], [2,1,1,1] and [1,1,1,1,1] do not qualify).
%p g:=product((1-x^(j*(j+1)))/(1-x^j),j=1..53): gser:=series(g,x=0,55): seq(coeff(gser,x,n),n=0..49); # _Emeric Deutsch_, Mar 04 2006
%p # second Maple program:
%p with(numtheory):
%p a:= proc(n) option remember; `if`(n=0, 1, add(add(
%p `if`(issqr(4*d+1), 0, d), d=divisors(j))*a(n-j), j=1..n)/n)
%p end:
%p seq(a(n), n=0..80); # _Alois P. Heinz_, Apr 01 2014
%t CoefficientList[Series[Product[Sum[x^(i j ), {i, 0, j}], {j, 1, 49}], {x, 0, 49}], x]
%t (* Second program: *)
%t a[n_] := a[n] = If[n == 0, 1, Sum[Sum[If[IntegerQ @ Sqrt[4*d+1], 0, d], {d, Divisors[j]}]*a[n-j], {j, 1, n}]/n]; Table[a[n], {n, 0, 80}] (* _Jean-François Alcover_, Jan 30 2017, after _Alois P. Heinz_ *)
%o (PARI) N=66; q='q+O('q^N); Vec( prod(n=1,N, sum(k=0,n,q^(k*n)) ) ) \\ _Joerg Arndt_, Apr 01 2014
%Y Cf. A000009, A000041, A002378, A033461, A117144, A087153.
%K nonn
%O 0,4
%A _Christian G. Bower_, Dec 19 1999