A090362
G.f. satisfies A^6 = BINOMIAL(A)^5 and also equals A090358^5.
Original entry on oeis.org
1, 5, 40, 460, 7220, 148276, 3831760, 120333680, 4460572870, 190679906990, 9230084185456, 498734395394840, 29740372199558420, 1939241402832412180, 137222625361036807760, 10470376552560151801616, 856818090423771231257245
Offset: 0
-
nmax = 16; sol = {a[0] -> 1};
Do[A[x_] = Sum[a[k] x^k, {k, 0, n}] /. sol; eq = CoefficientList[A[x]^6 - A[x/(1 - x)]^5/(1 - x)^5 + O[x]^(n + 1), x] == 0 /. sol; sol = sol ~Join~ Solve[eq][[1]], {n, 1, nmax}];
sol /. Rule -> Set;
a /@ Range[0, nmax] (* Jean-François Alcover, Nov 02 2019 *)
-
{a(n)=local(A); if(n<0,0,A=1+x+x*O(x^n); for(k=1,n,B=subst(A,x,x/(1-x))/(1-x)+x*O(x^n); A=A-A^6+B^5);polcoeff(A,n,x))}
A090359
Self-convolution equals the binomial transform of A090358: A^2 = BINOMIAL(A090358), where A090358^6 = BINOMIAL(A090358^5).
Original entry on oeis.org
1, 1, 4, 40, 640, 13816, 374636, 12229364, 466769330, 20391705290, 1003264704212, 54885373562372, 3304609250020008, 217139910688424400, 15461303963210314980, 1185856988993966140380, 97466557932008735970465
Offset: 0
-
{a(n)=local(A); if(n<1,0,A=1+x+x*O(x^n); for(k=1,n,B=subst(A^5,x,x/(1-x))/(1-x)+x*O(x^n); A=A-A^6+B);B=subst(A,x,x/(1-x))/(1-x)+x*O(x^n); polcoeff(B^(1/2),n,x))}
Original entry on oeis.org
1, 1, 5, 53, 855, 18471, 500651, 16334075, 623157898, 27214793210, 1338641476882, 73220527990322, 4408033139790830, 289619155572769806, 20620846574162157750, 1581511116725744872022, 129981009079318502310677
Offset: 0
-
{a(n)=local(A); if(n<1,0,A=1+x+x*O(x^n); for(k=1,n,B=subst(A^5,x,x/(1-x))/(1-x)+x*O(x^n); A=A-A^6+B);B=subst(A^2,x,x/(1-x))/(1-x)+x*O(x^n); polcoeff(B^(1/3),n,x))}
Original entry on oeis.org
1, 2, 12, 130, 2075, 44196, 1183492, 38254712, 1449394470, 62974889140, 3085705347128, 168283580164356, 10107641429213386, 662869582493716400, 47124829795282593000, 3609674673146922124600, 296355186635275737517175
Offset: 0
-
{a(n)=local(A); if(n<1,0,A=1+x+x*O(x^n); for(k=1,n,B=subst(A^5,x,x/(1-x))/(1-x)+x*O(x^n); A=A-A^6+B);B=subst(A^3,x,x/(1-x))/(1-x)+x*O(x^n); polcoeff(B^(1/2),n,x))}
A090351
G.f. satisfies A^3 = BINOMIAL(A^2).
Original entry on oeis.org
1, 1, 3, 15, 108, 1032, 12388, 179572, 3052986, 59555338, 1310677726, 32114051862, 866766965308, 25547102523604, 816335926158372, 28107705687291892, 1037367351120788551, 40852168787823027351, 1709792654612819858341
Offset: 0
A^3 = BINOMIAL(A090352), since A090352=A^2.
-
m:=40;
f:= func< n,x | Exp((&+[(&+[2^(j-1)*Factorial(j)* StirlingSecond(k,j)*x^k/k: j in [1..k]]): k in [1..n+2]])) >;
R:=PowerSeriesRing(Rationals(), m+1); // A090351
Coefficients(R!( f(m,x) )); // G. C. Greubel, Jun 08 2023
-
nmax = 18; sol = {a[0] -> 1};
Do[A[x_] = Sum[a[k] x^k, {k, 0, n}] /. sol; eq = CoefficientList[A[x]^3 - A[x/(1 - x)]^2/(1 - x) + O[x]^(n + 1), x] == 0 /. sol; sol = sol ~Join~ Solve[eq][[1]], {n, 1, nmax}];
sol /. Rule -> Set;
a /@ Range[0, nmax] (* Jean-François Alcover, Nov 02 2019 *)
With[{m=40}, CoefficientList[Series[Exp[Sum[Sum[2^(j-1)*j!* StirlingS2[k,j], {j,k}]*x^k/k, {k,m+1}]], {x,0,m}], x]] (* G. C. Greubel, Jun 08 2023 *)
-
{a(n) = my(A); if(n<0,0,A = 1+x +x*O(x^n); for(k=1,n, B = subst(A^2,x,x/(1-x))/(1-x) +x*O(x^n); A = A - A^3 + B); polcoef(A,n,x))}
for(n=0,25,print1(a(n),", "))
-
m=50
def f(n, x): return exp(sum(sum(2^(j-1)*factorial(j)* stirling_number2(k,j)*x^k/k for j in range(1,k+1)) for k in range(1,n+2)))
def A090351_list(prec):
P. = PowerSeriesRing(QQ, prec)
return P( f(m,x) ).list()
A090351_list(m-9) # G. C. Greubel, Jun 08 2023
A090353
G.f. satisfies A^4 = BINOMIAL(A^3).
Original entry on oeis.org
1, 1, 4, 28, 286, 3886, 66260, 1361972, 32784353, 904412593, 28124223808, 973106096392, 37073604836768, 1541948625066176, 69513081435903392, 3376138396206853792, 175739519606046355540, 9760024269508314079444
Offset: 0
A^4 = BINOMIAL(A090355), since A090355=A^3. Also, BINOMIAL(A) = A090354^2.
-
m:=40;
f:= func< n,x | Exp((&+[(&+[3^(j-1)*Factorial(j)* StirlingSecond(k,j)*x^k/k: j in [1..k]]): k in [1..n+2]])) >;
R:=PowerSeriesRing(Rationals(), m+1); // A090353
Coefficients(R!( f(m,x) )); // G. C. Greubel, Jun 09 2023
-
nmax = 17; sol = {a[0] -> 1};
Do[A[x_] = Sum[a[k] x^k, {k, 0, n}] /. sol; eq = CoefficientList[A[x]^4 - A[x/(1 - x)]^3/(1 - x) + O[x]^(n + 1), x] == 0 /. sol; sol = sol ~Join~ Solve[eq][[1]], {n, 1, nmax}];
sol /. Rule -> Set;
a /@ Range[0, nmax] (* Jean-François Alcover, Nov 02 2019 *)
With[{m=40}, CoefficientList[Series[Exp[Sum[Sum[3^(j-1)*j!* StirlingS2[k,j], {j,k}]*x^k/k, {k,m+1}]], {x,0,m}], x]] (* G. C. Greubel, Jun 09 2023 *)
-
{a(n) = my(A); if(n<0,0,A=1+x +x*O(x^n); for(k=1,n, B = subst(A^3,x,x/(1-x))/(1-x)+x*O(x^n); A = A - A^4 + B); polcoef(A,n,x))}
for(n=0,20,print1(a(n),", "))
-
m=50
def f(n, x): return exp(sum(sum(3^(j-1)*factorial(j)* stirling_number2(k,j)*x^k/k for j in range(1,k+1)) for k in range(1,n+2)))
def A090353_list(prec):
P. = PowerSeriesRing(QQ, prec)
return P( f(m,x) ).list()
A090353_list(m-9) # G. C. Greubel, Jun 09 2023
A090356
G.f. A(x) satisfies A(x)^5 = BINOMIAL(A(x)^4); that is, the binomial transform of the coefficients in A(x)^4 yields the coefficients in A(x)^5.
Original entry on oeis.org
1, 1, 5, 45, 595, 10475, 231255, 6148495, 191276600, 6815243040, 273601200136, 12217471594856, 600580173151560, 32224787998758280, 1873909224391774760, 117388347849375956328, 7880739469498103077588, 564440024187816634143380
Offset: 0
G.f.: A(x) = 1 + x + 5*x^2 + 45*x^3 + 595*x^4 + 10475*x^5 + 231255*x^6 + ...
The coefficients in A(x)^4 are given by A090357 and begin
A(x)^4: [1, 4, 26, 244, 3131, 52600, 1111940, ..., A090357(n), ...].
The binomial transform of A090357 yields the coefficients of A(x)^5:
A(x)^5: [1, 5, 35, 335, 4280, 70976, 1479800, ...]
as shown by
1 = 1*1,
5 = 1*1 + 1*4,
35 = 1*1 + 2*4 + 1*26,
335 = 1*1 + 3*4 + 3*26 + 1*244,
4280 = 1*1 + 4*4 + 6*26 + 4*244 + 1*3131, ...
-
m:=40;
f:= func< n,x | Exp((&+[(&+[4^(j-1)*Factorial(j)* StirlingSecond(k,j)*x^k/k: j in [1..k]]): k in [1..n+2]])) >;
R:=PowerSeriesRing(Rationals(), m+1); // A090356
Coefficients(R!( f(m,x) )); // G. C. Greubel, Jun 09 2023
-
nmax = 17; sol = {a[0] -> 1};
Do[A[x_] = Sum[a[k] x^k, {k, 0, n}] /. sol; eq = CoefficientList[A[x]^5 - A[x/(1 - x)]^4/(1 - x) + O[x]^(n + 1), x] == 0 /. sol; sol = sol ~Join~ Solve[eq][[1]], {n, 1, nmax}];
sol /. Rule -> Set;
a /@ Range[0, nmax] (* Jean-François Alcover, Nov 02 2019 *)
With[{m = 40}, CoefficientList[Series[Exp[Sum[Sum[4^(j-1)*j!* StirlingS2[k,j], {j,k}]*x^k/k, {k,m+1}]], {x,0,m}], x]] (* G. C. Greubel, Jun 09 2023 *)
-
{a(n)=local(A); if(n<1,0,A=1+x+x*O(x^n); for(k=1,n,B=subst(A^4,x,x/(1-x))/(1-x)+x*O(x^n); A=A-A^5+B);polcoeff(A,n,x))}
-
m=40
def f(n, x): return exp(sum(sum(4^(j-1)*factorial(j)* stirling_number2(k,j)*x^k/k for j in range(1,k+1)) for k in range(1,n+2)))
def A090356_list(prec):
P. = PowerSeriesRing(QQ, prec)
return P( f(m,x) ).list()
A090356_list(m) # G. C. Greubel, Jun 09 2023
A258377
O.g.f. satisfies A^2(z) = 1/(1 - z)*( BINOMIAL(BINOMIAL(A(z))) ).
Original entry on oeis.org
1, 3, 13, 79, 649, 6955, 93813, 1539991, 29884881, 669628819, 17005862301, 482399018527, 15108642099673, 517599894435643, 19247498583665029, 771922934908235751, 33206411983713679009, 1525025984109289947171, 74466779211331635306029, 3852255519421356879419631
Offset: 0
Cf.
A019538,
A080253,
A084784,
A084785,
A090351,
A090352,
A090353,
A090355,
A090356,
A090357,
A090358,
A090362,
A145901,
A258378 (N = 2),
A258379 (N = 3),
A258380 (N = 4),
A258381 (N = 5).
-
#A258377
with(combinat):
#recursively define row polynomials R(n,x) of A145901
R := proc (n, x) option remember; if n = 0 then 1 else 1 + x*add(binomial(n, i)*2^(n-i)*R(i,x), i = 0..n-1) end if; end proc:
#define a family of sequences depending on an integer parameter k
a := proc (n, k) option remember; if n = 0 then 1 else 1/n*add(R(i+1,k)*a(n-1-i,k), i = 0..n-1) end if; end proc:
# display the case k = 1
seq(a(n,1), n = 0..19);
-
R[n_, x_] := R[n, x] = If[n==0, 1, 1+x*Sum[Binomial[n, i]*2^(n-i)*R[i, x], {i, 0, n-1}]];
a[n_, k_] := a[n, k] = If[n==0, 1, 1/n*Sum[R[i+1, k]*a[n-1-i, k], {i, 0, n-1}]];
a[n_] := a[n, 1];
a /@ Range[0, 19] (* Jean-François Alcover, Oct 02 2019 *)
A384305
Expansion of Product_{k>=1} 1/(1 - k*x)^((5/6)^k).
Original entry on oeis.org
1, 30, 615, 11260, 205695, 4013406, 88035585, 2255192280, 68859250020, 2506898720040, 107238427737876, 5281094776037040, 293625956135692020, 18139856902224931080, 1229886945212115522060, 90641666662687182976896, 7206758883035555464430370, 614391718014749017022916060
Offset: 0
-
terms = 20; A[] = 1; Do[A[x] = -5*A[x] + 6*A[x/(1-x)]^(5/6) / (1-x)^5 + O[x]^j // Normal, {j, 1, terms}]; CoefficientList[A[x], x] (* Vaclav Kotesovec, May 31 2025 *)
-
my(N=20, x='x+O('x^N)); Vec(exp(6*sum(k=1, N, sum(j=0, k, 5^j*j!*stirling(k, j, 2))*x^k/k)))
Showing 1-9 of 9 results.
Comments