A291820
G.f. A(x,y) satisfies: A( x - x*y*A(x,y), y) = x + x*(1-y)*A(x,y), where the coefficients T(n,k) of x^n*y^k form a triangle read by rows n>=1, for k=0..n-1.
Original entry on oeis.org
1, 1, 0, 1, 2, 0, 1, 7, 5, 0, 1, 16, 38, 14, 0, 1, 30, 157, 189, 42, 0, 1, 50, 477, 1245, 904, 132, 0, 1, 77, 1197, 5616, 8791, 4242, 429, 0, 1, 112, 2632, 19881, 55566, 57854, 19723, 1430, 0, 1, 156, 5250, 59327, 265204, 491947, 363880, 91366, 4862, 0, 1, 210, 9714, 155783, 1035442, 3062271, 4039551, 2220933, 423124, 16796, 0, 1, 275, 16929, 370205, 3472513, 15217674, 31979723, 31463341, 13285415, 1963169, 58786, 0, 1, 352, 28094, 811877, 10331673, 63678254, 197983540, 310618856, 235959185, 78419541, 9138416, 208012, 0
Offset: 1
G.f.: A(x,y) = x + x^2 + (2*y + 1)*x^3 + (5*y^2 + 7*y + 1)*x^4 +
(14*y^3 + 38*y^2 + 16*y + 1)*x^5 +
(42*y^4 + 189*y^3 + 157*y^2 + 30*y + 1)*x^6 +
(132*y^5 + 904*y^4 + 1245*y^3 + 477*y^2 + 50*y + 1)*x^7 +
(429*y^6 + 4242*y^5 + 8791*y^4 + 5616*y^3 + 1197*y^2 + 77*y + 1)*x^8 +
(1430*y^7 + 19723*y^6 + 57854*y^5 + 55566*y^4 + 19881*y^3 + 2632*y^2 + 112*y + 1)*x^9 +
(4862*y^8 + 91366*y^7 + 363880*y^6 + 491947*y^5 + 265204*y^4 + 59327*y^3 + 5250*y^2 + 156*y + 1)*x^10 +
(16796*y^9 + 423124*y^8 + 2220933*y^7 + 4039551*y^6 + 3062271*y^5 + 1035442*y^4 + 155783*y^3 + 9714*y^2 + 210*y + 1)*x^11 +
(58786*y^10 + 1963169*y^9 + 13285415*y^8 + 31463341*y^7 + 31979723*y^6 + 15217674*y^5 + 3472513*y^4 + 370205*y^3 + 16929*y^2 + 275*y + 1)*x^12 +...
such that
A( x - x*y*A(x,y), y) = x + x*(1-y)*A(x,y).
Also,
A(x,y) = x + Z*A(Z, y) where Z = y*A(x,y) + (1-y)*x.
...
This triangle of coefficients T(n,k) of x^n*y^k in g.f. A(x,y) begins:
1;
1, 0;
1, 2, 0;
1, 7, 5, 0;
1, 16, 38, 14, 0;
1, 30, 157, 189, 42, 0;
1, 50, 477, 1245, 904, 132, 0;
1, 77, 1197, 5616, 8791, 4242, 429, 0;
1, 112, 2632, 19881, 55566, 57854, 19723, 1430, 0;
1, 156, 5250, 59327, 265204, 491947, 363880, 91366, 4862, 0;
1, 210, 9714, 155783, 1035442, 3062271, 4039551, 2220933, 423124, 16796, 0;
1, 275, 16929, 370205, 3472513, 15217674, 31979723, 31463341, 13285415, 1963169, 58786, 0;
1, 352, 28094, 811877, 10331673, 63678254, 197983540, 310618856, 235959185, 78419541, 9138416, 208012, 0;
1, 442, 44759, 1666522, 27896583, 232505790, 1014785477, 2355151627, 2859824058, 1721756609, 458956233, 42718416, 742900, 0; ...
RELATED SEQUENCES.
Given T(n,k) is the coefficient of x^n*y^k in g.f. A(x,y),
if b(n) = Sum_{k=0..n-1} T(n,k) * p^k * q^(n-k-1)
then B(x) = Sum_{n>=1} b(n)*x^n satisfies
(E1) B(x - p*x*B(x)) = x + (q-p)*x*B(x)
(E2) B(x) = x + Z*B(Z) where Z = p*B(x) + (q-p)*x.
...
G.f.s of columns of this triangle begin:
C.0: 1/(1-x)
C.1: (2 - x)/(1-x)^4
C.2: (5 + 3*x - 4*x^2 + x^3)/(1-x)^7
C.3: (14 + 49*x - 15*x^2 - 9*x^3 + 6*x^4 - x^5)/(1-x)^10
C.4: (42 + 358*x + 315*x^2 - 217*x^3 + 30*x^4 + 18*x^5 - 8*x^6 + x^7)/(1-x)^13
C.5: (132 + 2130*x + 5822*x^2 + 1403*x^3 - 1681*x^4 + 602*x^5 - 50*x^6 - 30*x^7 + 10*x^8 - x^9)/(1-x)^16
C.6: (429 + 11572*x + 62502*x^2 + 82763*x^3 + 2951*x^4 - 9760*x^5 + 5395*x^6 - 1329*x^7 + 75*x^8 + 45*x^9 - 12*x^10 + x^11)/(1-x)^19
C.7: (1430 + 59906*x + 541211*x^2 + 1506161*x^3 + 1217687*x^4 + 16416*x^5 - 35746*x^6 + 36682*x^7 - 13502*x^8 + 2550*x^9 - 105*x^10 - 63*x^11 + 14*x^12 - x^13)/(1-x)^22
C.8: (4862 + 301574*x + 4165915*x^2 + 19578410*x^3 + 34788033*x^4 + 20899306*x^5 + 1681742*x^6 + 174039*x^7 + 195964*x^8 - 103084*x^9 + 28953*x^10 - 4444*x^11 + 140*x^12 + 84*x^13 - 16*x^14 + x^15)/(1-x)^25
...
Thus A(x, y*(1-x)^3)*(1-x) = x + 2*y*x^3 + (5*y^2 - y)*x^4 + (14*y^3 + 3*y^2)*x^5 + (42*y^4 + 49*y^3 - 4*y^2)*x^6 + (132*y^5 + 358*y^4 - 15*y^3 + y^2)*x^7 +...
-
nmax = 13; A[x_] = x;
Do[A[x_] = x + (y A[x] + (1-y) x) A[y A[x] + (1-y) x] + x O[x]^n // Normal // Expand // Collect[#, x]&, {n, nmax}];
T[n_, k_] := SeriesCoefficient[A[x], {x, 0, n}, {y, 0, k}];
Table[T[n, k], {n, 1, nmax}, {k, 0, n-1}] // Flatten (* Jean-François Alcover, Oct 20 2019 *)
-
{T(n, k) = my(A=x); for(i=1, n, A = x + subst(x*A, x, y*A + (1-y)*x +x*O(x^n)) ); polcoeff(polcoeff(A, n, x), k, y)}
for(n=1, 12, for(k=0, n-1, print1(T(n, k), ", ")); print(""))
A291814
G.f. A(x) satisfies: A(x - 3*x*A(x)) = x - 2*x*A(x).
Original entry on oeis.org
1, 1, 7, 67, 769, 10009, 143359, 2218255, 36625657, 639659737, 11741022235, 225390779647, 4508109360985, 93665093491381, 2016669357747667, 44905700922069463, 1032419000661778213, 24472819932819733957, 597384952530618840715, 15000294032677574361955, 387082666821619977435277, 10256260095368150955828565, 278811213889895147327704519, 7770474960716476086765483619
Offset: 1
G.f.: A(x) = x + x^2 + 7*x^3 + 67*x^4 + 769*x^5 + 10009*x^6 + 143359*x^7 + 2218255*x^8 + 36625657*x^9 + 639659737*x^10 + 11741022235*x^11 + 225390779647*x^12 +...
such that A(x - 3*x*A(x)) = x - 2*x*A(x).
RELATED SERIES.
A(x - 3*x*A(x)) = x - 2*x^2 - 2*x^3 - 14*x^4 - 134*x^5 - 1538*x^6 - 20018*x^7 +...
which equals x - 2*x*A(x).
Series_Reversion( x - 3*x*A(x) ) = x + 3*x^2 + 21*x^3 + 201*x^4 + 2307*x^5 + 30027*x^6 + 430077*x^7 + 6654765*x^8 +...
which equals 3*A(x) - 2*x.
A( 3*A(x) - 2*x ) = x + 4*x^2 + 34*x^3 + 382*x^4 + 5038*x^5 + 74134*x^6 + 1184650*x^7 + 20224990*x^8 + 364994554*x^9 + 6911857450*x^10 + 136622440786*x^11 + 2807805653098*x^12 +...
which equals (A(x) - x) / (3*A(x) - 2*x).
-
{a(n) = my(A=x); for(i=1, n, A = (1/3)*serreverse( x - 3*x*A +x*O(x^n) ) + 2*x/3 ); polcoeff(A, n)}
for(n=1, 30, print1(a(n), ", "))
A291815
G.f. A(x) satisfies: A(x - 4*x*A(x)) = x - 3*x*A(x).
Original entry on oeis.org
1, 1, 9, 109, 1569, 25481, 454105, 8730373, 178996865, 3881556561, 88477557289, 2109927671453, 52443846331297, 1354646602217945, 36275862587452281, 1005099719255707829, 28765965099599741953, 849204340574458575777, 25827102287376124267593, 808349897942417046805197, 26011340193853765710238241, 859773626049480606121078057, 29168437337569276216572259097
Offset: 1
G.f.: A(x) = x + x^2 + 9*x^3 + 109*x^4 + 1569*x^5 + 25481*x^6 + 454105*x^7 + 8730373*x^8 + 178996865*x^9 + 3881556561*x^10 + 88477557289*x^11 + 2109927671453*x^12 +...
such that A(x - 4*x*A(x)) = x - 3*x*A(x).
RELATED SERIES.
A(x - 4*x*A(x)) = x - 3*x^2 - 3*x^3 - 27*x^4 - 327*x^5 - 4707*x^6 - 76443*x^7 +...
which equals x - 3*x*A(x).
Series_Reversion( x - 4*x*A(x) ) = x + 4*x^2 + 36*x^3 + 436*x^4 + 6276*x^5 + 101924*x^6 + 1816420*x^7 + 34921492*x^8 +...
which equals 4*A(x) - 3*x.
A( 4*A(x) - 3*x ) = x + 5*x^2 + 53*x^3 + 741*x^4 + 12153*x^5 + 222405*x^6 + 4421501*x^7 + 93949493*x^8 + 2110952881*x^9 + 49786323589*x^10 + 1225967873349*x^11 + 31395927333829*x^12 +...
which equals (A(x) - x) / (4*A(x) - 3*x).
-
{a(n) = my(A=x); for(i=1, n, A = (1/4)*serreverse( x - 4*x*A +x*O(x^n) ) + 3*x/4 ); polcoeff(A, n)}
for(n=1, 30, print1(a(n), ", "))
A291816
G.f. A(x) satisfies: A(x - 3*x*A(x)) = x - x*A(x).
Original entry on oeis.org
1, 2, 16, 182, 2524, 39992, 699520, 13231034, 266985280, 5694001172, 127481465536, 2981125793144, 72532301230672, 1830526849868000, 47802726801684544, 1289123410465365782, 35841130838977837348, 1025903099063974343984, 30195807234087904770952, 912951678159786641659796, 28327442752528049524839856, 901289532361030971832330544, 29382621186595702051011638128
Offset: 1
G.f.: A(x) = x + 2*x^2 + 16*x^3 + 182*x^4 + 2524*x^5 + 39992*x^6 + 699520*x^7 + 13231034*x^8 + 266985280*x^9 + 5694001172*x^10 +...
such that A(x - 3*x*A(x)) = x - x*A(x).
RELATED SERIES.
A(x - 3*x*A(x)) = x - x^2 - 2*x^3 - 16*x^4 - 182*x^5 - 2524*x^6 - 39992*x^7 - 699520*x^8 +...
which equals x - x*A(x).
Series_Reversion( x - 3*x*A(x) ) = x + 3*x^2 + 24*x^3 + 273*x^4 + 3786*x^5 + 59988*x^6 + 1049280*x^7 + 19846551*x^8 +...
which equals (3/2)*A(x) - x/2.
A( (3*A(x) - x)/2 ) = x + 5*x^2 + 52*x^3 + 713*x^4 + 11458*x^5 + 205160*x^6 + 3984304*x^7 + 82576109*x^8 + 1807215616*x^9 + 41461917398*x^10 +...
which equals (A(x) - x) / (3*A(x) - x).
-
{a(n) = my(A=x); for(i=1, n, A = (2/3)*serreverse( x - 3*x*A +x*O(x^n) ) + x/3 ); polcoeff(A, n)}
for(n=1, 30, print1(a(n), ", "))
A291817
G.f. A(x) satisfies: A(x - 3*x*A(x)) = x + x*A(x).
Original entry on oeis.org
1, 4, 40, 580, 10312, 209752, 4707952, 114128308, 2946787192, 80268150808, 2290811949904, 68149759850680, 2105074730357968, 67308231895605520, 2222306773263886624, 75615701295449074084, 2647156154616207962920, 95215874465318554556776, 3514739264129342710455184, 133012394550946993742673304, 5156112210399927390763631056, 204570581814658024128260509360
Offset: 1
G.f.: A(x) = x + 4*x^2 + 40*x^3 + 580*x^4 + 10312*x^5 + 209752*x^6 + 4707952*x^7 + 114128308*x^8 + 2946787192*x^9 + 80268150808*x^10 +...
such that A(x - 3*x*A(x)) = x + x*A(x).
RELATED SERIES.
A(x - 3*x*A(x)) = x + x^2 + 4*x^3 + 40*x^4 + 580*x^5 + 10312*x^6 + 209752*x^7 + 4707952*x^8 +...
which equals x + x*A(x).
Series_Reversion( x - 3*x*A(x) ) = x + 3*x^2 + 30*x^3 + 435*x^4 + 7734*x^5 + 157314*x^6 + 3530964*x^7 + 85596231*x^8 +...
which equals (3/4)*A(x) + x/4.
A( (3*A(x) + x)/4 ) = x + 7*x^2 + 94*x^3 + 1651*x^4 + 33886*x^5 + 773458*x^6 + 19117780*x^7 + 503529979*x^8 + 13983485770*x^9 + 406470316978*x^10 +...
which equals (A(x) - x) / (3*A(x) + x).
-
{a(n) = my(A=x); for(i=1, n, A = (4/3)*serreverse( x - 3*x*A +x*O(x^n) ) - x/3 ); polcoeff(A, n)}
for(n=1, 30, print1(a(n), ", "))
A291818
G.f. A(x) satisfies: A(x - 4*x*A(x)) = x - x*A(x).
Original entry on oeis.org
1, 3, 33, 519, 9969, 218907, 5307201, 139123215, 3889995297, 114928234611, 3563543673825, 115375173490839, 3885328265571345, 135675583665864843, 4900856792035006593, 182756242210436579871, 7023982500750575903553, 277842871320960134512611, 11297961688442941015761825, 471773677417286920645721895, 20211597594930636918024401457, 887652829316087359743197592315
Offset: 1
G.f.: A(x) = x + 3*x^2 + 33*x^3 + 519*x^4 + 9969*x^5 + 218907*x^6 + 5307201*x^7 + 139123215*x^8 + 3889995297*x^9 + 114928234611*x^10 +...
such that A(x - 4*x*A(x)) = x - x*A(x).
RELATED SERIES.
A(x - 4*x*A(x)) = x - x^2 - 3*x^3 - 33*x^4 - 519*x^5 - 9969*x^6 - 218907*x^7 - 5307201*x^8 +...
which equals x - x*A(x).
Series_Reversion( x - 4*x*A(x) ) = x + 4*x^2 + 44*x^3 + 692*x^4 + 13292*x^5 + 291876*x^6 + 7076268*x^7 + 185497620*x^8 +...
which equals (4/3)*A(x) - x/3.
A( (4*A(x) - x)/3 ) = x + 7*x^2 + 101*x^3 + 1919*x^4 + 42713*x^5 + 1058967*x^6 + 28469325*x^7 + 816617535*x^8 + 24729787889*x^9 + 784895219495*x^10 +...
which equals (A(x) - x) / (4*A(x) - x).
-
{a(n) = my(A=x); for(i=1, n, A = (3/4)*serreverse( x - 4*x*A +x*O(x^n) ) + x/4 ); polcoeff(A, n)}
for(n=1, 30, print1(a(n), ", "))
A291819
G.f. A(x) satisfies: A(x - x*A(x)) = x + 3*x*A(x).
Original entry on oeis.org
1, 4, 24, 196, 1944, 21944, 272080, 3627412, 51288200, 761782104, 11805102064, 189901153112, 3158767322992, 54165347282960, 955189096759776, 17289056525343716, 320678326091307448, 6087009196570756488, 118109764108446889008, 2340448760238788518488, 47324471620802426563376, 975739573623235107473968, 20500725692629852174532192, 438679922664144046444438488, 9555430871381022848971028208
Offset: 1
G.f.: A(x) = x + 4*x^2 + 24*x^3 + 196*x^4 + 1944*x^5 + 21944*x^6 + 272080*x^7 + 3627412*x^8 + 51288200*x^9 + 761782104*x^10 +...
such that A(x - x*A(x)) = x + 3*x*A(x).
RELATED SERIES.
A(x - x*A(x)) = x + 3*x^2 + 12*x^3 + 72*x^4 + 588*x^5 + 5832*x^6 + 65832*x^7 + 816240*x^8 +...
which equals x + 3*x*A(x).
Series_Reversion( x - x*A(x) ) = x + x^2 + 6*x^3 + 49*x^4 + 486*x^5 + 5486*x^6 + 68020*x^7 + 906853*x^8 +...
which equals (1/4)*A(x) + 3*x/4.
A( (A(x) + 3*x)/4 ) = x + 5*x^2 + 38*x^3 + 369*x^4 + 4158*x^5 + 51870*x^6 + 698036*x^7 + 9974297*x^8 + 149755186*x^9 + 2345335606*x^10 +...
which equals (A(x) - x) / (A(x) + 3*x).
-
{a(n) = my(A=x); for(i=1, n, A = 4*serreverse( x - x*A +x*O(x^n) ) - 3*x ); polcoeff(A, n)}
for(n=1, 30, print1(a(n), ", "))
Showing 1-7 of 7 results.
Comments