A294989 Number of compositions (ordered partitions) of 1 into exactly 10n+1 powers of 1/(n+1).

1, 68958747, 27212315953140892, 39880061006390454401626995, 110656003660578876500875377620844376, 423205992807070499591372608204571223421862945, 1944053748730350081768848916806347783770184147756976500

Offset

0,2

Links

Alois P. Heinz, Table of n, a(n) for n = 0..101

Formula

a(n) ~ 10^(10*n + 3/2) / (2*Pi*n)^(9/2). - Vaclav Kotesovec, Sep 20 2019

Maple

b:= proc(n, r, p, k) option remember;
      `if`(n<r, 0, `if`(r=0, `if`(n=0, p!, 0), add(
       b(n-j, k*(r-j), p+j, k)/j!, j=0..min(n, r))))
    end:
a:= n-> (k-> `if`(n=0, 1, b(k*n+1, 1, 0, n+1)))(10):
seq(a(n), n=0..10);

Mathematica

b[n_, r_, p_, k_] := b[n, r, p, k] = If[n < r, 0, If[r == 0, If[n == 0, p!, 0], Sum[b[n - j, k*(r - j), p + j, k]/j!, {j, 0, Min[n, r]}]]];
a[n_] := If[n == 0, 1, b[#*n + 1, 1, 0, n + 1]]&[10];
Table[a[n], {n, 0, 10}] (* Jean-François Alcover, May 21 2018, translated from Maple *)

Crossrefs

Row n=10 of A294746.

Sequence in context: A251615 A321706 A172608 * A204405 A186618 A186597

Adjacent sequences: A294986 A294987 A294988 * A294990 A294991 A294992

Keyword

nonn

Author

Alois P. Heinz, Nov 12 2017

Status

approved