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

%I #41 Feb 24 2022 18:22:36

%S 1,1,1,5,7,37,79,173,101,127,1033,1571,200069,2564519,5126711,

%T 25661369,532393,431100529,1855391,1533985991,48977868113,

%U 342880481117,342289639579,435979161889,1308720597671,373092965489,7824703695283,24141028973,31250466692609

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

%C a(n)/A322381(n) = A007838(n)/A000142(n) is the probability that a random permutation of [n] has distinct cycle sizes. - _Geoffrey Critzer_, Feb 23 2022

%H Alois P. Heinz, <a href="/A322380/b322380.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.

%H Philippe Flajolet and Robert Sedgewick, <a href="http://algo.inria.fr/flajolet/Publications/AnaCombi/anacombi.html">Analytic Combinatorics</a>, Cambridge Univ. Press, 2009, page 137.

%H A. Knopfmacher and J. N. Ridley, <a href="http://dx.doi.org/10.1137/0406031">Reciprocal sums over partitions and compositions</a>, SIAM J. Discrete Math. 6 (1993), no. 3, 388-399.

%H D. H. Lehmer, <a href="http://matwbn.icm.edu.pl/ksiazki/aa/aa21/aa21123.pdf">On reciprocally weighted partitions</a>, Acta Arithmetica XXI (1972), 379-388.

%H D. Zeilberger, N. Zeilberger, <a href="http://sites.math.rutgers.edu/~zeilberg/mamarim/mamarimhtml/fcp.html">Fractional Counting of Integer Partitions</a>, 2018.

%F Limit_{n->infinity} a(n)/A322381(n) = exp(-gamma) = A080130.

%F Sum_{n>=0} a(n)/A322381(n)*x^n = Product_{i>=1} (1 + x^i/i). - _Geoffrey Critzer_, Feb 23 2022

%e 1/1, 1/1, 1/2, 5/6, 7/12, 37/60, 79/120, 173/280, 101/168, 127/210, 1033/1680, 1571/2640, 200069/332640, 2564519/4324320, 5126711/8648640, ... = A322380/A322381

%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-> numer(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_] := Numerator[b[n, n]];

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

%Y Denominators: A322381.

%Y Cf. A000009, A000142, A007838, A022629, A080130, A177208, A177209, A322364, A322365, A323290, A323291, A323339, A323340.

%K nonn,frac

%O 0,4

%A _Alois P. Heinz_, Dec 05 2018