|
%I #18 Jan 29 2021 19:02:31
%S 1,3,9,19,39,72,128,216,354,563,876,1335,1998,2946,4284,6154,8742,
%T 12294,17129,23667,32442,44151,59682,80169,107054,142167,187812,
%U 246895,323058,420852,545958,705438,908043,1164609,1488504,1896174,2407836,3048255,3847716
%N Number of colored integer partitions of n such that three colors are used and parts differ by size or by color.
%C With offset 0 convolution cube of A000009(k+1). - _George Beck_, Jan 29 2021
%H Vaclav Kotesovec, <a href="/A327381/b327381.txt">Table of n, a(n) for n = 3..10000</a> (terms 3..5000 from Alois P. Heinz)
%H Wikipedia, <a href="https://en.wikipedia.org/wiki/Partition_(number_theory)">Partition (number theory)</a>
%F a(n) ~ exp(Pi*sqrt(n)) / (8 * n^(3/4)). - _Vaclav Kotesovec_, Sep 14 2019
%F G.f.: (-1 + Product_{m >= 1} (1 + x^m))^3. - _George Beck_, Jan 29 2021
%p b:= proc(n, i, k) option remember; `if`(n=0, 1, `if`(i<1, 0, add((t->
%p b(t, min(t, i-1), k)*binomial(k, j))(n-i*j), j=0..min(k, n/i))))
%p end:
%p a:= n-> (k-> add(b(n$2, k-i)*(-1)^i*binomial(k, i), i=0..k))(3):
%p seq(a(n), n=3..45);
%t b[n_, i_, k_] := b[n, i, k] = If[n==0, 1, If[i<1, 0, Sum[With[{t = n - i j}, b[t, Min[t, i - 1], k] Binomial[k, j]], {j, 0, Min[k, n/i]}]]];
%t a[n_] := With[{k = 3}, Sum[b[n, n, k-i] (-1)^i Binomial[k, i], {i, 0, k}]];
%t a /@ Range[3, 45] (* _Jean-François Alcover_, Dec 15 2020, after _Alois P. Heinz_ *)
%Y Column k=3 of A308680.
%Y Cf. A000009.
%K nonn
%O 3,2
%A _Alois P. Heinz_, Sep 03 2019
|
