login
Denominator of the sum of inverse products of parts in all strict partitions of n.
8

%I #13 Feb 25 2020 08:14:40

%S 1,1,2,6,12,60,120,280,168,210,1680,2640,332640,4324320,8648640,

%T 43243200,900900,735134400,3150576,2618916300,83805321600,

%U 586637251200,586637251200,749592043200,2248776129600,642507465600,13492656777600,41644002400,53970627110400

%N Denominator of the sum of inverse products of parts in all strict partitions of n.

%H Alois P. Heinz, <a href="/A322381/b322381.txt">Table of n, a(n) for n = 0..1268</a>

%H Andreas B. G. Blobel, <a href="https://arxiv.org/abs/1904.07808">An Asymptotic Form of the Generating Function Prod_{k=1,oo} (1+x^k/k)</a>, arXiv:1904.07808 [math.CO], 2019.

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

%p b(n, i-1) +b(n-i, min(i-1, n-i))/i))

%p end:

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

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

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

%t a[n_] := Denominator[b[n, n]];

%t a /@ Range[0, 30] (* _Jean-François Alcover_, Feb 25 2020, after _Alois P. Heinz_ *)

%Y See A322380 for more information.

%K nonn,frac

%O 0,3

%A _Alois P. Heinz_, Dec 05 2018