|
%I #29 May 06 2022 13:13:51
%S 1,1,4,4,12,18,44,72,158,288,604,1146,2332,4528,9126,17944,35940,
%T 71130,142132,282344,563630,1121936,2239060,4462530,8906236,17764160,
%U 35458774,70761520,141272876,282025466,563159588,1124543256,2245918406,4485670168,8960061076
%N Number of compositions of n with the multiplicity of the first part odd.
%H Alois P. Heinz, <a href="/A331606/b331606.txt">Table of n, a(n) for n = 1..1000</a>
%H M. Archibald, A. Blecher, A. Knopfmacher, M. E. Mays, <a href="https://cs.uwaterloo.ca/journals/JIS/VOL23/Archibald/arch3.html">Inversions and Parity in Compositions of Integers</a>, J. Int. Seq., Vol. 23 (2020), Article 20.4.1.
%F G.f.: Sum_{i>=1} (1-x)*x^i/(2*(-2*x^(i+1)+2*x^i-2*x+1)) + x/(2*(1-2*x)).
%F a(n) = A011782(n) - A331609(n). - _Alois P. Heinz_, Jan 23 2020
%e For n=3, a(4)=4 as we count 4, 3+1, 1+3 and 2+1+1.
%p b:= proc(n, p, t) option remember; `if`(n=0, t,
%p add(b(n-j, p, `if`(p=j, 1-t, t)), j=1..n))
%p end:
%p a:= n-> add(b(n-j, j, 1), j=1..n):
%p seq(a(n), n=1..38); # _Alois P. Heinz_, Jan 23 2020
%t gf[x_] := x/(2 (1 - 2 x)) + Sum[(1 - x) x^i/(2 (-2 x^(i + 1) + 2 x^i - 2 x + 1)) , {i, 1, 40}]; CL := CoefficientList[Series[gf[x], {x, 0, 35}], x];
%t Drop[CL, 1] (* _Peter Luschny_, Jan 23 2020 *)
%Y Cf. A011782, A331609 (similar, with even).
%K nonn
%O 1,3
%A _Arnold Knopfmacher_, Jan 22 2020
|
