login
a(n) = number of partitions of n into an odd number of parts, the least being 2; also a(n+2) = number of partitions of n into an even number of parts, each >=2.
5

%I #18 May 15 2023 11:13:35

%S 0,0,1,0,0,0,1,1,2,2,4,4,6,7,11,12,17,20,28,33,44,52,69,82,105,126,

%T 161,191,239,286,355,423,520,618,755,896,1084,1285,1549,1829,2190,

%U 2583,3079,3621,4297,5041,5960,6977,8214,9595,11264,13123,15353,17854,20828

%N a(n) = number of partitions of n into an odd number of parts, the least being 2; also a(n+2) = number of partitions of n into an even number of parts, each >=2.

%H Alois P. Heinz, <a href="/A027188/b027188.txt">Table of n, a(n) for n = 0..1000</a>

%F G.f.: sum(n>=0, q^(2*n+1)/prod(k=1..n, 1-q^(2*k)) * q^(2*n+1)/prod(k=1..n, 1-q^(2*k-1)) ). [_Joerg Arndt_, Feb 27 2014]

%F a(n+2) + A027194(n+2) = A002865(n). - _R. J. Mathar_, Jun 16 2018

%F G.f.: x^2 * Sum_{k>=0} x^(4*k)/Product_{j=1..2*k} (1-x^j). - _Seiichi Manyama_, May 15 2023

%p b:= proc(n, i, t) option remember; `if`(n=0, t, `if`(i<2, 0,

%p b(n, i-1, t) +`if`(i>n, 0, b(n-i, i, 1-t))))

%p end:

%p a:= n-> b(n-2$2, 1):

%p seq(a(n), n=0..80); # _Alois P. Heinz_, Feb 27 2014

%t b[n_, i_, t_] := b[n, i, t] = If[n==0, t, If[i<2, 0, b[n, i-1, t] + If[i>n, 0, b[n-i, i, 1-t]]]]; a[n_] := b[n-2, n-2, 1]; Table[a[n], {n, 0, 80}] (* _Jean-François Alcover_, Apr 08 2015, after _Alois P. Heinz_ *)

%o (PARI) gf=sum(n=0, N, q^(2*n+1)/prod(k=1,n, 1-q^(2*k)) * q^(2*n+1)/prod(k=1,n, 1-q^(2*k-1)) );

%o concat([0], Vec(gf) ) \\ _Joerg Arndt_, Feb 27 2014

%K nonn

%O 0,9

%A _Clark Kimberling_