login
The OEIS is supported by the many generous donors to the OEIS Foundation.

 

Logo
Hints
(Greetings from The On-Line Encyclopedia of Integer Sequences!)
A294746 Number A(n,k) of compositions (ordered partitions) of 1 into exactly k*n+1 powers of 1/(k+1); square array A(n,k), n>=0, k>=0, read by antidiagonals. 21

%I #27 Sep 20 2019 05:13:09

%S 1,1,1,1,1,1,1,1,3,1,1,1,10,13,1,1,1,35,217,75,1,1,1,126,4245,8317,

%T 525,1,1,1,462,90376,1239823,487630,4347,1,1,1,1716,2019836,216456376,

%U 709097481,40647178,41245,1,1,1,6435,46570140,41175714454,1303699790001,701954099115,4561368175,441675,1

%N Number A(n,k) of compositions (ordered partitions) of 1 into exactly k*n+1 powers of 1/(k+1); square array A(n,k), n>=0, k>=0, read by antidiagonals.

%C Row r >= 2, is asymptotic to r^(r*n + 3/2) / (2*Pi*n)^((r-1)/2). - _Vaclav Kotesovec_, Sep 20 2019

%H Alois P. Heinz, <a href="/A294746/b294746.txt">Antidiagonals n = 0..45, flattened</a>

%F A(n,k) = [x^((k+1)^n)] (Sum_{j=0..k*n+1} x^((k+1)^j))^(k*n+1) for k>0, A(n,0) = 1.

%e A(3,1) = 13: [1/4,1/4,1/4,1/4], [1/2,1/4,1/8,1/8], [1/2,1/8,1/4,1/8], [1/2,1/8,1/8,1/4], [1/4,1/2,1/8,1/8], [1/4,1/8,1/2,1/8], [1/4,1/8,1/8,1/2], [1/8,1/2,1/4,1/8], [1/8,1/2,1/8,1/4], [1/8,1/4,1/2,1/8], [1/8,1/4,1/8,1/2], [1/8,1/8,1/2,1/4], [1/8,1/8,1/4,1/2].

%e Square array A(n,k) begins:

%e 1, 1, 1, 1, 1, 1, ...

%e 1, 1, 1, 1, 1, 1, ...

%e 1, 3, 10, 35, 126, 462, ...

%e 1, 13, 217, 4245, 90376, 2019836, ...

%e 1, 75, 8317, 1239823, 216456376, 41175714454, ...

%e 1, 525, 487630, 709097481, 1303699790001, 2713420774885145, ...

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

%p `if`(n<r, 0, `if`(r=0, `if`(n=0, p!, 0), add(

%p b(n-j, k*(r-j), p+j, k)/j!, j=0..min(n, r))))

%p end:

%p A:= (n, k)-> `if`(k=0, 1, b(k*n+1, 1, 0, k+1)):

%p seq(seq(A(n, d-n), n=0..d), d=0..10);

%t 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]}]]];

%t A[n_, k_] := If[k == 0, 1, b[k*n + 1, 1, 0, k + 1]];

%t Table[A[n, d - n], {d, 0, 10}, {n, 0, d}] // Flatten (* _Jean-François Alcover_, Apr 30 2018, after _Alois P. Heinz_ *)

%Y Columns k=0-10 give: A000012, A007178(n+1), A294850, A294851, A294852, A294853, A294854, A294855, A294856, A294857, A294858.

%Y Rows n=0+1, 2-10 give: A000012, A001700, A294982, A294983, A294984, A294985, A294986, A294987, A294988, A294989.

%Y Main diagonal gives: A294747.

%Y Cf. A294775.

%K nonn,tabl

%O 0,9

%A _Alois P. Heinz_, Nov 07 2017

Lookup | Welcome | Wiki | Register | Music | Plot 2 | Demos | Index | Browse | More | WebCam
Contribute new seq. or comment | Format | Style Sheet | Transforms | Superseeker | Recents
The OEIS Community | Maintained by The OEIS Foundation Inc.

License Agreements, Terms of Use, Privacy Policy. .

Last modified April 24 00:30 EDT 2024. Contains 371917 sequences. (Running on oeis4.)