A301305
G.f. L(x,y) satisfies: L(x,y) = x * (1 + y*x*L'(x,y)) / (1 + x*L'(x,y)) where L'(x,y) = d/dx L(x,y), as a triangle read by rows.
Original entry on oeis.org
1, -1, 1, 3, -5, 2, -14, 33, -25, 6, 85, -261, 292, -140, 24, -621, 2363, -3516, 2546, -892, 120, 5236, -23872, 44537, -43405, 23228, -6444, 720, -49680, 264860, -596396, 733983, -532095, 226644, -52356, 5040, 521721, -3193029, 8448004, -12605668, 11586756, -6707208, 2383248, -474144, 40320, -5994155, 41506739, -126480376, 222424796, -248535142, 182793154, -88379152, 27046632, -4745376, 362880, 74701055, -578419961, 1998774636, -4056699966, 5351696394, -4791391134, 2945757656, -1226765624, 330797184, -52079040, 3628800
Offset: 1
G.f. L(x,y) = Sum_{n>=1, k=0..n-1} T(n,k) * x^n*y^k begins:
L(x,y) = x + (-1 + y)*x^2 + (3 - 5*y + 2*y^2)*x^3 + (-14 + 33*y - 25*y^2 + 6*y^3)*x^4 + (85 - 261*y + 292*y^2 - 140*y^3 + 24*y^4)*x^5 + (-621 + 2363*y - 3516*y^2 + 2546*y^3 - 892*y^4 + 120*y^5)*x^6 + (5236 - 23872*y + 44537*y^2 - 43405*y^3 + 23228*y^4 - 6444*y^5 + 720*y^6)*x^7 + (-49680 + 264860*y - 596396*y^2 + 733983*y^3 - 532095*y^4 + 226644*y^5 - 52356*y^6 + 5040*y^7)*x^8 + ...
where L = L(x,y) satisfies:
L = x*(1 + y*x*L') / (1 + x*L').
TRIANGLE.
This triangle of coefficients T(n,k) in L(x,y) begins:
1;
-1, 1;
3, -5, 2;
-14, 33, -25, 6;
85, -261, 292, -140, 24;
-621, 2363, -3516, 2546, -892, 120;
5236, -23872, 44537, -43405, 23228, -6444, 720;
-49680, 264860, -596396, 733983, -532095, 226644, -52356, 5040;
521721, -3193029, 8448004, -12605668, 11586756, -6707208, 2383248, -474144, 40320;
-5994155, 41506739, -126480376, 222424796, -248535142, 182793154, -88379152, 27046632, -4745376, 362880; ...
LIMITS.
In this triangle, the largest real root of the n-th row polynomial in y converges to the constant t = 2.845344903202547217277843362090557097661... (A301389).
RELATED SERIES.
exp(L(x,y)) = 1 + x + (-1 + 2*y)*x^2/2! + (13 - 24*y + 12*y^2)*x^3/3! + (-263 + 660*y - 540*y^2 + 144*y^3)*x^4/4! + (8381 - 26800*y + 31380*y^2 - 15840*y^3 + 2880*y^4)*x^5/5! + (-379409 + 1485870*y - 2280180*y^2 + 1706520*y^3 - 619200*y^4 + 86400*y^5)*x^6/6! + (22915369 - 106759128*y + 203726880*y^2 - 203269920*y^3 + 111449520*y^4 - 31691520*y^5 + 3628800*y^6)*x^7/7! + ...
satisfies: [x^n] exp(-n*L(x,y)) = ((y-1)*(n-1) - 1) * [x^(n-1)] exp(-n*L(x,y)) for n>=1.
-
{T(n,k) = my(L=x); for(i=0,n, L = x*(1 + y*x*L')/(1 + x*L' +x*O(x^n)) ); polcoeff(polcoeff(L,n,x),k,y)}
for(n=1,12, for(k=0,n-1, print1(T(n,k),", "));print(""))
A301385
O.g.f. A(x) satisfies: A(x) = x * (1 + 3*x*A'(x)) / (1 + x*A'(x)).
Original entry on oeis.org
1, 2, 6, 22, 94, 474, 2974, 24630, 271710, 3799570, 63378806, 1208997078, 25736584670, 602485683530, 15356903176110, 423032451327510, 12518043710674878, 395909541133928226, 13325077980379707238, 475466006418129789206, 17926802213221278261726, 712095926927360739006522, 29722097317161256669118142, 1300445348644716445771904502
Offset: 1
G.f.: A(x) = x + 2*x^2 + 6*x^3 + 22*x^4 + 94*x^5 + 474*x^6 + 2974*x^7 + 24630*x^8 + 271710*x^9 + 3799570*x^10 + ...
where
A(x) = x*(1 + 3*x*A'(x)) / (1 + x*A'(x)).
RELATED SERIES.
A'(x) = 1 + 4*x + 18*x^2 + 88*x^3 + 470*x^4 + 2844*x^5 + 20818*x^6 + 197040*x^7 + 2445390*x^8 + 37995700*x^9 + ...
exp(A(x)) = 1 + x + 5*x^2/2! + 49*x^3/3! + 745*x^4/4! + 16001*x^5/5! + 472621*x^6/6! + 19659025*x^7/7! + 1211940689*x^8/8! + ... + A301386*x^n/n! + ...
-
Rest[CoefficientList[AsymptoticDSolveValue[{A[x] == x*(1 + 3*x*A'[x])/(1 + x*A'[x]), A[1] == 1}, A[x], {x, 0, 20}], x]] (* Requires Mathematica version 11.3 or later *) (* Vaclav Kotesovec, Mar 20 2018 *)
-
{a(n) = my(A=x); for(i=0,n, A = x*(1 + 3*x*A')/(1 +x*A' +x*O(x^n)) ); polcoeff(A,n)}
for(n=1,30,print1(a(n),", "))
A301387
E.g.f. A(x) satisfies: [x^n] A(x)^n = (2*n - 1) * [x^(n-1)] A(x)^n for n>=1.
Original entry on oeis.org
1, 1, 5, 73, 2185, 108881, 8012941, 809101945, 106751544593, 17777715999265, 3641538959077141, 899203905438721961, 263360370074829490585, 90266849931865728247153, 35797787561010642901996445, 16266200052343102435435467481, 8396767484613064218299219696161, 4887270793881840799764234723049025, 3185975842406919756361250644217366053
Offset: 0
E.g.f.: A(x) = 1 + x + 5*x^2/2! + 73*x^3/3! + 2185*x^4/4! + 108881*x^5/5! + 8012941*x^6/6! + 809101945*x^7/7! + 106751544593*x^8/8! + 17777715999265*x^9/9! + ...
The table of coefficients in A(x)^n begins:
n=1: [(1), (1), 5/2, 73/6, 2185/24, 108881/120, 8012941/720, ...];
n=2: [1, (2), (6), 88/3, 638/3, 10288/5, 1110424/45, ...];
n=3: [1, 3, (21/2), (105/2), 2979/8, 140241/40, 3288981/80, ...];
n=4: [1, 4, 16, (248/3), (1736/3), 79768/15, 2744776/45, ...];
n=5: [1, 5, 45/2, 725/6, (20185/24), (60555/8), 12237185/144, ...];
n=6: [1, 6, 30, 168, 1170, (51744/5), (569184/5), ...]; ...
in which the coefficients in parenthesis are related by
1 = 1*(1); 6 = 3*(2); 105/2 = 5*(21/2); 1736/3 = 7*(248/3); 60555/8 = 9*(20185/24); 569184/5 = 11*(51744/5); ...
illustrating: [x^n] A(x)^n = (2*n - 1) * [x^(n-1)] A(x)^n.
LOGARITHMIC PROPERTY.
The logarithm of the e.g.f. is an integer power series in x satisfying
log(A(x)) = x * (1 + x*A'(x)/A(x)) / (1 - x*A'(x)/A(x));
explicitly,
log(A(x)) = x + 2*x^2 + 6*x^3 + 22*x^4 + 94*x^5 + 474*x^6 + 2974*x^7 + 24630*x^8 + 271710*x^9 + 3799570*x^10 + ... + A301388(n)*x^n + ...
-
{a(n) = my(A=1); for(i=1, n, A = exp( x*(A + x*A')/(A - x*A' +x*O(x^n)) ) ); n!*polcoeff(A, n)}
for(n=0, 25, print1(a(n), ", "))
Showing 1-3 of 3 results.
Comments