A306024 Number A(n,k) of length-n restricted growth strings (RGS) with growth <= k and first element in [k]; square array A(n,k), n>=0, k>=0, read by antidiagonals.
1, 1, 0, 1, 1, 0, 1, 2, 2, 0, 1, 3, 7, 5, 0, 1, 4, 15, 31, 15, 0, 1, 5, 26, 95, 164, 52, 0, 1, 6, 40, 214, 717, 999, 203, 0, 1, 7, 57, 405, 2096, 6221, 6841, 877, 0, 1, 8, 77, 685, 4875, 23578, 60619, 51790, 4140, 0, 1, 9, 100, 1071, 9780, 67354, 297692, 652595, 428131, 21147, 0
Offset: 0
A305963 Number of length-n restricted growth strings (RGS) with growth <= n and fixed first element.
1, 1, 3, 22, 305, 6756, 216552, 9416240, 530764089, 37498693555, 3235722405487, 334075729235172, 40587204883652869, 5722676826879812177, 925590727478445526747, 170032646641380554970304, 35173161711207720944899921, 8132124409499796317194563900
Offset: 0
Keywords
Examples
a(2) = 3: 11, 12, 13. a(3) = 22: 111, 112, 113, 114, 121, 122, 123, 124, 125, 131, 132, 133, 134, 135, 136, 141, 142, 143, 144, 145, 146, 147.
Links
- Alois P. Heinz, Table of n, a(n) for n = 0..266
Programs
-
Maple
b:= proc(n, k, m) option remember; `if`(n=0, 1, add(b(n-1, k, max(m, j)), j=1..m+k)) end: a:= n-> b(n$2, 1-n): seq(a(n), n=0..20); # second Maple program: a:= n-> `if`(n=0, 1, (n-1)!*coeff(series(exp(x+add( (exp(j*x)-1)/j, j=1..n)), x, n), x, n-1)): seq(a(n), n=0..20);
-
Mathematica
b[n_, k_, m_] := b[n, k, m] = If[n == 0, 1, Sum[b[n - 1, k, Max[m, j]], {j, 1, m + k}]]; a[n_] := b[n, n, 1 - n]; Table[a[n], {n, 0, 20}] (* Jean-François Alcover, Mar 21 2022, after Alois P. Heinz *)
Formula
a(n) = (n-1)! * [x^(n-1)] exp(x + Sum_{j=1..n} (exp(j*x)-1)/j) for n > 0, a(0) = 1.
a(n) = A305962(n,n).
Comments
Examples
Links
Crossrefs
Programs
Maple
Mathematica
Formula