A332720 - OEIS

%I #16 Jan 06 2021 02:12:31

%S 1,1,5,19,59,150,349,745,1515,2936,5514,10036,17851,31039,53006,88943,

%T 147057,239701,385885,613855,966137,1505137,2323124,3553914,5392315,

%U 8117758,12131618,18003740,26543030,38886999,56633453,82009410,118113488,169229009,241264461

%N Index position of {3}^n within the list of partitions of 3n in canonical ordering.

%C The canonical ordering of partitions is described in A080577.

%H Alois P. Heinz, <a href="/A332720/b332720.txt">Table of n, a(n) for n = 0..4000</a>

%H Wikipedia, <a href="https://www.wikipedia.org/wiki/integer_partition">Integer Partition</a>

%F a(n) ~ exp(Pi*sqrt(2*n)) / (4*3^(3/2)*n). - _Vaclav Kotesovec_, Feb 28 2020

%e a(2) = 5, because 33 has position 5 within the list of partitions of 6 in canonical ordering: 6, 51, 42, 411, 33, 321, 3111, 222, ... .

%p b:= proc(n) option remember;

%p      `if`(n=0, 1, b(n-1)+g(3*n, 2))

%p     end:

%p g:= proc(n, i) option remember; `if`(n=0 or i=1, 1,

%p      `if`(i<1, 0, g(n-i, min(n-i, i))+g(n, i-1)))

%p     end:

%p a:= n-> g(3*n$2)-b(n)+1:

%p seq(a(n), n=0..35);

%t b[n_] := b[n] = If[n == 0, 1, b[n - 1] + g[3n, 2]];

%t g[n_, i_] := g[n, i] = If[n == 0 || i == 1, 1, If[i < 1, 0, g[n - i, Min[n - i, i]] + g[n, i - 1]]];

%t a[n_] := g[3n, 3n] - b[n] + 1;

%t a /@ Range[0, 35] (* _Jean-François Alcover_, Jan 06 2021, after _Alois P. Heinz_ *)

%Y Cf. A000041, A080577, A322761, A330661, A332706, A332719.

%K nonn

%O 0,3

%A _Alois P. Heinz_, Feb 20 2020