A370419
A(n, k) = 2^n*Pochhammer(k/2, n). Square array read by ascending antidiagonals.
Original entry on oeis.org
1, 0, 1, 0, 1, 1, 0, 3, 2, 1, 0, 15, 8, 3, 1, 0, 105, 48, 15, 4, 1, 0, 945, 384, 105, 24, 5, 1, 0, 10395, 3840, 945, 192, 35, 6, 1, 0, 135135, 46080, 10395, 1920, 315, 48, 7, 1, 0, 2027025, 645120, 135135, 23040, 3465, 480, 63, 8, 1
Offset: 0
The array starts:
[0] 1, 1, 1, 1, 1, 1, 1, 1, 1, ...
[1] 0, 1, 2, 3, 4, 5, 6, 7, 8, ...
[2] 0, 3, 8, 15, 24, 35, 48, 63, 80, ...
[3] 0, 15, 48, 105, 192, 315, 480, 693, 960, ...
[4] 0, 105, 384, 945, 1920, 3465, 5760, 9009, 13440, ...
[5] 0, 945, 3840, 10395, 23040, 45045, 80640, 135135, 215040, ...
.
Seen as the triangle T(n, k) = A(n - k, k):
[0] 1;
[1] 0, 1;
[2] 0, 1, 1;
[3] 0, 3, 2, 1;
[4] 0, 15, 8, 3, 1;
[5] 0, 105, 48, 15, 4, 1;
[6] 0, 945, 384, 105, 24, 5, 1;
.
From _Werner Schulte_, Mar 07 2024: (Start)
Illustrating the LU decomposition of A:
/ 1 \ / 1 1 1 1 1 ... \ / 1 1 1 1 1 ... \
| 0 1 | | 1 2 3 4 ... | | 0 1 2 3 4 ... |
| 0 3 2 | * | 1 3 6 ... | = | 0 3 8 15 24 ... |
| 0 15 18 6 | | 1 4 ... | | 0 15 48 105 192 ... |
| 0 105 174 108 24 | | 1 ... | | 0 105 384 945 1920 ... |
| . . . | | . . . | | . . . |. (End)
-
A := (n, k) -> 2^n*pochhammer(k/2, n):
for n from 0 to 5 do seq(A(n, k), k = 0..9) od;
T := (n, k) -> A(n - k, k): seq(seq(T(n, k), k = 0..n), n = 0..9);
# Using the exponential generating functions of the columns:
EGFcol := proc(k, len) local egf, ser, n; egf := (1 - 2*x)^(-k/2);
ser := series(egf, x, len+2): seq(n!*coeff(ser, x, n), n = 0..len) end:
seq(lprint(EGFcol(n, 9)), n = 0..8);
# Using the generating polynomials for the rows:
P := (n, x) -> local k; add(Stirling1(n, k)*(-2)^(n - k)*x^k, k=0..n):
seq(lprint([n], seq(P(n, k), k = 0..8)), n = 0..5);
# Implementing the comment of Werner Schulte about the LU decomposition of A:
with(LinearAlgebra):
L := Matrix(7, 7, (n, k) -> A371025(n - 1, k - 1)):
U := Matrix(7, 7, (n, k) -> binomial(n - 1, k - 1)):
MatrixMatrixMultiply(L, Transpose(U)); # Peter Luschny, Mar 08 2024
-
A370419[n_, k_] := 2^n*Pochhammer[k/2, n];
Table[A370419[n-k, k], {n, 0, 10}, {k, 0, n}] (* Paolo Xausa, Mar 06 2024 *)
-
def A(n, k): return 2**n * rising_factorial(k/2, n)
for n in range(6): print([A(n, k) for k in range(9)])
A288952
Number of relaxed compacted binary trees of right height at most one with empty sequences between branch nodes on level 0.
Original entry on oeis.org
1, 0, 1, 2, 15, 92, 835, 8322, 99169, 1325960, 19966329, 332259290, 6070777999, 120694673748, 2594992240555, 59986047422378, 1483663965460545, 39095051587497488, 1093394763005554801, 32347902448449172530, 1009325655965539561231, 33125674098690460236620
Offset: 0
- Muniru A Asiru, Table of n, a(n) for n = 0..100
- Antoine Genitrini, Bernhard Gittenberger, Manuel Kauers and Michael Wallner, Asymptotic Enumeration of Compacted Binary Trees, arXiv:1703.10031 [math.CO], 2017.
- Michael Wallner, A bijection of plane increasing trees with relaxed binary trees of right height at most one, arXiv:1706.07163 [math.CO], 2017.
Cf.
A001147 (relaxed compacted binary trees of right height at most one).
Cf.
A082161 (relaxed compacted binary trees of unbounded right height).
-
a := [1,0];; for n in [3..10^2] do a[n] := (n-2)*a[n-1] + (n-2)^2*a[n-2]; od; a; # Muniru A Asiru, Jan 26 2018
-
a:=proc(n) option remember: if n=0 then 1 elif n=1 then 0 elif n>=2 then (n-1)*procname(n-1)-(n-1)^2*procname(n-2) fi; end:
seq(a(n),n=0..100); # Muniru A Asiru, Jan 26 2018
-
Fold[Append[#1, (#2 - 1) Last[#1] + #1[[#2 - 1]] (#2 - 1)^2] &, {1, 0}, Range[2, 21]] (* Michael De Vlieger, Jan 28 2018 *)
A288953
Number of relaxed compacted binary trees of right height at most one with minimal sequences between branch nodes except after the last branch node on level 0.
Original entry on oeis.org
1, 1, 3, 10, 51, 280, 1995, 15120, 138075, 1330560, 14812875, 172972800, 2271359475, 31135104000, 471038042475, 7410154752000, 126906349444875, 2252687044608000, 43078308695296875, 851515702861824000, 17984171447178811875, 391697223316439040000
Offset: 0
Denote by L the leaf and by o nodes. Every node has exactly two out-going edges or pointers. Internal edges are denoted by - or |. Pointers are omitted and may point to any node further right. The root is at level 0 at the very left.
The general structure is
L-o-o-o-o-o-o-o-o
| | | | |
o o o o o.
For n=0 the a(0)=1 solution is L.
For n=1 the a(1)=1 solution is L-o.
For n=2 the a(2)=3 solutions are
L-o-o L-o
|
o
2 + 1 solutions of this shape with pointers.
Cf.
A288954 (variation with additional initial sequence).
Cf.
A177145 (variation without final sequence).
Cf.
A001147 (relaxed compacted binary trees of right height at most one).
Cf.
A082161 (relaxed compacted binary trees of unbounded right height).
A153271
Triangle T(n, k) = Product_{j=0..k} (j*n + prime(m)), with T(n, 0) = prime(m) and m = 3, read by rows.
Original entry on oeis.org
5, 5, 30, 5, 35, 315, 5, 40, 440, 6160, 5, 45, 585, 9945, 208845, 5, 50, 750, 15000, 375000, 11250000, 5, 55, 935, 21505, 623645, 21827575, 894930575, 5, 60, 1140, 29640, 978120, 39124800, 1838865600, 99298742400, 5, 65, 1365, 39585, 1464645, 65909025, 3493178325, 213083877825, 14702787569925
Offset: 0
Triangle begins as:
5;
5, 30;
5, 35, 315;
5, 40, 440, 6160;
5, 45, 585, 9945, 208845;
5, 50, 750, 15000, 375000, 11250000;
5, 55, 935, 21505, 623645, 21827575, 894930575;
Sequences related to m values:
-
m:=3;
function T(n,k)
if k eq 0 then return NthPrime(m);
else return (&*[j*n + NthPrime(m): j in [0..k]]);
end if; return T; end function;
[T(n,k): k in [0..n], n in [0..10]]; // G. C. Greubel, Dec 03 2019
-
m:=3; seq(seq(`if`(k=0, ithprime(m), mul(j*n + ithprime(m), j=0..k)), k=0..n), n=0..10); # G. C. Greubel, Dec 03 2019
-
T[n_, k_, m_]:= If[k==0, Prime[m], Product[j*n + Prime[m], {j,0,k}]];
Table[T[n,k,3], {n,0,10}, {k,0,n}]//Flatten
-
T(n,k) = my(m=3); if(k==0, prime(m), prod(j=0,k, j*n + prime(m)) ); \\ G. C. Greubel, Dec 03 2019
-
def T(n, k):
m=3
if (k==0): return nth_prime(m)
else: return product(j*n + nth_prime(m) for j in (0..k))
[[T(n, k) for k in (0..n)] for n in (0..10)] # G. C. Greubel, Dec 03 2019
A288954
Number of relaxed compacted binary trees of right height at most one with minimal sequences between branch nodes except before the first and after the last branch node on level 0.
Original entry on oeis.org
1, 1, 3, 13, 79, 555, 4605, 42315, 436275, 4894155, 60125625, 794437875, 11325612375, 172141044075, 2793834368325, 48009995908875, 874143494098875, 16757439016192875, 338309837281040625, 7157757510792763875, 158706419654857449375, 3673441093896736036875
Offset: 2
Cf.
A288953 (variation without initial sequence).
Cf.
A177145 (variation without initial and final sequence).
Cf.
A001147 (relaxed compacted binary trees of right height at most one).
Cf.
A082161 (relaxed compacted binary trees of unbounded right height).
-
terms = 22; egf = 1/(3(1-z))(1/Sqrt[1-z^2] + (3z^3 - z^2 - 2z + 2)/((1-z)(1-z^2))) + O[z]^terms;
CoefficientList[egf, z] Range[0, terms-1]! (* Jean-François Alcover, Dec 13 2018 *)
A092081
Triangle of certain double factorials.
Original entry on oeis.org
1, 1, 2, 1, 3, 8, 1, 4, 15, 48, 1, 5, 24, 105, 384, 1, 6, 35, 192, 945, 3840, 1, 7, 48, 315, 1920, 10395, 46080, 1, 8, 63, 480, 3465, 23040, 135135, 645120, 1, 9, 80, 693, 5760, 45045, 322560, 2027025, 10321920, 1, 10, 99, 960, 9009, 80640, 675675, 5160960
Offset: 0
A316666
Number of simple relaxed compacted binary trees of right height at most one with no sequences on level 1 and no final sequences on level 0.
Original entry on oeis.org
1, 0, 1, 3, 15, 87, 597, 4701, 41787, 413691, 4512993, 53779833, 695000919, 9680369943, 144560191149, 2303928046437, 39031251610227, 700394126116851, 13270625547477177, 264748979672169681, 5547121478845459983, 121784530649198053263, 2795749225338111831429, 66981491857058929294653
Offset: 0
- G. C. Greubel, Table of n, a(n) for n = 0..448
- Antoine Genitrini, Bernhard Gittenberger, Manuel Kauers and Michael Wallner, Asymptotic Enumeration of Compacted Binary Trees, arXiv:1703.10031 [math.CO], 2017.
- Michael Wallner, A bijection of plane increasing trees with relaxed binary trees of right height at most one, arXiv:1706.07163 [math.CO], 2017.
-
m:=25; R:=PowerSeriesRing(Rationals(), m); b:=Coefficients(R!( (3*Exp(-x) + x-2)/(1-x)^2 )); [Factorial(n-1)*b[n]: n in [1..m]]; // G. C. Greubel, Dec 12 2018
-
aseq := n-> 3*round((n+2)*n!/exp(1))-(n+2)*n!: bseq := n-> (n+2)*n!- 2* round((n+2)*n!/exp(1)): s := (a,b,n)-> a*aseq(n) + b*bseq( n): seq(s(1,0,n),n = 0..20); # Gary Detlefs, Dec 11 2018
-
terms = 24;
CoefficientList[(3E^-z+z-2)/(1-z)^2 + O[z]^terms, z] Range[0, terms-1]! (* Jean-François Alcover, Sep 14 2018 *)
-
Vec(serlaplace((3*exp(-x + O(x^25)) + x - 2)/(1 - x)^2)) \\ Andrew Howroyd, Jul 10 2018
A102147
Second Eulerian transform of 1, 2, 3, 4, 5, ... (A000027).
Original entry on oeis.org
1, 1, 5, 35, 315, 3465, 45045, 675675, 11486475, 218243025, 4583103525, 105411381075, 2635284526875, 71152682225625, 2063427784543125, 63966261320836875, 2110886623587616875, 73881031825566590625
Offset: 1
- R. L. Graham, D. E. Knuth and O. Patashnik, Concrete Mathematics. Addison-Wesley, Reading, MA, 1990, p. 256.
Apparently equals
A051577(n-2), n > 1.
Comments