cp's OEIS Frontend

This is a front-end for the Online Encyclopedia of Integer Sequences, made by Christian Perfect. The idea is to provide OEIS entries in non-ancient HTML, and then to think about how they're presented visually. The source code is on GitHub.

Showing 1-6 of 6 results.

A360345 G.f. A(x) satisfies: [x^n] A(x)^(n+1) = [x^n] (1 + x*A(x)^(2*n+1))^(n+1) for n >= 0.

Original entry on oeis.org

1, 1, 5, 62, 1214, 31269, 973485, 34993597, 1412846469, 62926155294, 3053566438307, 160005640085764, 8992869671470675, 539298198547460797, 34364052537634696986, 2318526571023659653665, 165143229278977841236029, 12385688813185721332861730, 975844100444710104444582984
Offset: 0

Views

Author

Paul D. Hanna, Feb 05 2023

Keywords

Examples

			G.f.: A(x) = 1 + x + 5*x^2 + 62*x^3 + 1214*x^4 + 31269*x^5 + 973485*x^6 + 34993597*x^7 + 1412846469*x^8 + 62926155294*x^9 + ...
RELATED SERIES.
G.f. A(x) = B(x/A(x)) where B(x) = B(x*A(x)) begins:
B(x) = 1 + x + 6*x^2 + 78*x^3 + 1543*x^4 + 39810*x^5 + 1239252*x^6 + 44537587*x^7 + 1798314384*x^8 + ... + b(n)*x^n + ...
such that b(n) = [x^n] (1 + x*A(x)^(2*n+1))^(n+1) / (n+1),
as well as b(n) = [x^n] A(x)^(n+1) / (n+1),
so that b(n) begin:
[1/1, 2/2, 18/3, 312/4, 7715/5, 238860/6, 8674764/7, 356300696/8, ...].
ILLUSTRATION OF DEFINITION.
The table of coefficients of x^k in A(x)^(n+1) begins:
n=0: [1, 1,  5,  62,  1214,  31269,   973485,  34993597, ...];
n=1: [1, 2, 11, 134,  2577,  65586,  2025492,  72397390, ...];
n=2: [1, 3, 18, 217,  4104, 103212,  3161648, 112357788, ...];
n=3: [1, 4, 26, 312,  5811, 144428,  4387978, 155030276, ...];
n=4: [1, 5, 35, 420,  7715, 189536,  5710930, 200579975, ...];
n=5: [1, 6, 45, 542,  9834, 238860,  7137401, 249182232, ...];
n=6: [1, 7, 56, 679, 12187, 292747,  8674764, 301023241, ...];
n=7: [1, 8, 68, 832, 14794, 351568, 10330896, 356300696, ...]; ...
Compare to the table of coefficients in (1 + x*A(x)^(2*n+1))^(n+1):
n=0: [1, 1,   1,    5,    62,   1214,    31269,    973485, ...];
n=1: [1, 2,   7,   42,   479,   8750,   216258,   6562156, ...];
n=2: [1, 3,  18,  136,  1560,  26895,   633608,  18631701, ...];
n=3: [1, 4,  34,  312,  3767,  62888,  1412530,  40031684, ...];
n=4: [1, 5,  55,  595,  7715, 128041,  2763270,  75234930, ...];
n=5: [1, 6,  81, 1010, 14172, 238860,  5016947, 131313798, ...];
n=6: [1, 7, 112, 1582, 24059, 418166,  8674764, 219340759, ...];
n=7: [1, 8, 148, 2336, 38450, 696216, 14466592, 356300696, ...]; ...
to see that the main diagonals of the tables are the same.
		

Crossrefs

Programs

  • PARI
    {a(n) = my(A=[1]); for(m=1, n, A=concat(A, 0); A[m+1] = (Vec((1+x*Ser(A)^(2*m+1))^(m+1))[m+1] - Vec(Ser(A)^(m+1))[m+1])/(m+1) ); A[n+1]}
    for(n=0, 20, print1(a(n), ", "))

Formula

G.f. A(x) = Sum_{n>=0} a(n)*x^n satisfies:
(1) [x^n] A(x)^(n+1) = [x^n] (1 + x*A(x)^(2*n+1))^(n+1) for n>=0.
(2) A(x) = Sum_{n>=0} b(n) * x^n/A(x)^n, where b(n) = [x^n] (1 + x*A(x)^(2*n+1))^(n+1) / (n+1).
a(n) ~ c * d^n * n! * n^alpha, where d = 3.93464558322824528799..., alpha = 1.635402029299..., c = 0.0308525091280143... - Vaclav Kotesovec, Feb 06 2023

A360346 G.f. A(x) satisfies: [x^n] A(x)^(n+1) = [x^n] (1 + x*A(x)^(2*n+2))^(n+1) for n >= 0.

Original entry on oeis.org

1, 1, 6, 82, 1724, 47223, 1555047, 58892186, 2496826094, 116434989450, 5900151126856, 322048641354617, 18810964989814291, 1169843128503194025, 77145176721564799777, 5376524285402806746719, 394887654026596322701724, 30489608056346314234108286, 2469347798211941105406473481
Offset: 0

Views

Author

Paul D. Hanna, Feb 05 2023

Keywords

Examples

			G.f.: A(x) = 1 + x + 6*x^2 + 82*x^3 + 1724*x^4 + 47223*x^5 + 1555047*x^6 + 58892186*x^7 + 2496826094*x^8 + 116434989450*x^9 + ...
RELATED SERIES.
G.f. A(x) = B(x/A(x)) where B(x) = B(x*A(x)) begins:
B(x) = 1 + x + 7*x^2 + 101*x^3 + 2161*x^4 + 59544*x^5 + 1965132*x^6 + 74504861*x^7 + 3161424763*x^8 + ... + b(n)*x^n + ...
such that b(n) = [x^n] (1 + x*A(x)^(2*n+2))^(n+1) / (n+1),
as well as b(n) = [x^n] A(x)^(n+1) / (n+1),
so that b(n) begin:
[1/1, 2/2, 21/3, 404/4, 10805/5, 357264/6, 13755924/7, 596038888/8, ...].
ILLUSTRATION OF DEFINITION.
The table of coefficients of x^k in A(x)^(n+1) begins:
n=0: [1, 1,  6,   82,  1724,  47223,  1555047,  58892186, ...];
n=1: [1, 2, 13,  176,  3648,  98878,  3231952, 121743878, ...];
n=2: [1, 3, 21,  283,  5790, 155319,  5039055, 188787837, ...];
n=3: [1, 4, 30,  404,  8169, 216924,  6985240, 260270488, ...];
n=4: [1, 5, 40,  540, 10805, 284096,  9079965, 336452690, ...];
n=5: [1, 6, 51,  692, 13719, 357264, 11333293, 417610542, ...];
n=6: [1, 7, 63,  861, 16933, 436884, 13755924, 504036226, ...];
n=7: [1, 8, 76, 1048, 20470, 523440, 16359228, 596038888, ...]; ...
Compare to the table of coefficients in (1 + x*A(x)^(2*n+2))^(n+1):
n=0: [1, 1,   2,   13,   176,   3648, 98878, 3231952, ...];
n=1: [1, 2,   9,   68,   884,  17386,   454318,  14493920, ...];
n=2: [1, 3,  21,  190,  2508,  47406,  1190949,  36928479, ...];
n=3: [1, 4,  38,  404,  5585, 103464,  2504568,  75227160, ...];
n=4: [1, 5,  60,  735, 10805, 200001,  4698210, 136509465, ...];
n=5: [1, 6,  87, 1208, 19011, 357264,  8227591, 231595008, ...];
n=6: [1, 7, 119, 1848, 31199, 602427, 13755924, 376756199, ...];
n=7: [1, 8, 156, 2680, 48518, 970712, 22218108, 596038888, ...]; ...
to see that the main diagonals of the tables are the same.
		

Crossrefs

Programs

  • PARI
    {a(n) = my(A=[1]); for(m=1, n, A=concat(A, 0); A[m+1] = (Vec((1+x*Ser(A)^(2*m+2))^(m+1))[m+1] - Vec(Ser(A)^(m+1))[m+1])/(m+1) ); A[n+1]}
    for(n=0, 20, print1(a(n), ", "))

Formula

G.f. A(x) = Sum_{n>=0} a(n)*x^n satisfies:
(1) [x^n] A(x)^(n+1) = [x^n] (1 + x*A(x)^(2*n+2))^(n+1) for n>=0.
(2) A(x) = Sum_{n>=0} b(n) * x^n/A(x)^n, where b(n) = [x^n] (1 + x*A(x)^(2*n+2))^(n+1) / (n+1).
a(n) ~ c * d^n * n! * n^alpha, where d = 3.93464558322824528799..., alpha = 2.153902930660..., c = 0.01676305987174... - Vaclav Kotesovec, Feb 06 2023

A360347 G.f. A(x) satisfies: [x^n] A(x)^(n+1) = [x^n] (1 + x*A(x)^(2*n+3))^(n+1) for n >= 0.

Original entry on oeis.org

1, 1, 7, 105, 2366, 68776, 2390230, 95166058, 4228436480, 206090296497, 10887958126763, 618187895371965, 37479711430699245, 2414492049517400164, 164626564026042206780, 11841796830661101527618, 896184803460067359995232, 71189783172592806474895908
Offset: 0

Views

Author

Paul D. Hanna, Feb 05 2023

Keywords

Comments

Sequences with g.f. A(x,k) such that [x^n] A(x,k)^(n+1) = [x^n] (1 + x*A(x,k)^(2*n+k))^(n+1) have a rate of growth: a(n) ~ c(k) * d^n * n! * n^alpha(k), where d = 3.93464558322824528799... (independent on k) and alpha(k) = 1.1169011279372... + k*0.518500901361... - Vaclav Kotesovec, Feb 06 2023

Examples

			G.f.: A(x) = 1 + x + 7*x^2 + 105*x^3 + 2366*x^4 + 68776*x^5 + 2390230*x^6 + 95166058*x^7 + 4228436480*x^8 + 206090296497*x^9 + ...
RELATED SERIES.
G.f. A(x) = B(x/A(x)) where B(x) = B(x*A(x)) begins:
B(x) = 1 + x + 8*x^2 + 127*x^3 + 2927*x^4 + 85892*x^5 + 2998264*x^6 + 119665415*x^7 + 5325877575*x^8 + ... + b(n)*x^n + ...
such that b(n) = [x^n] (1 + x*A(x)^(2*n+3))^(n+1) / (n+1),
as well as b(n) = [x^n] A(x)^(n+1) / (n+1),
so that b(n) begin:
[1/1, 2/2, 24/3, 508/4, 14635/5, 515352/6, 20987848/7, 957323320/8, ...].
ILLUSTRATION OF DEFINITION.
The table of coefficients of x^k in A(x)^(n+1) begins:
n=0: [1, 1,  7,  105,  2366,  68776,  2390230,  95166058, ...];
n=1: [1, 2, 15,  224,  4991, 143754,  4962161, 196572300, ...];
n=2: [1, 3, 24,  358,  7896, 225396,  7727644, 304572936, ...];
n=3: [1, 4, 34,  508, 11103, 314192, 10699244, 419541832, ...];
n=4: [1, 5, 45,  675, 14635, 410661, 13890275, 541873525, ...];
n=5: [1, 6, 57,  860, 18516, 515352, 17314836, 671984280, ...];
n=6: [1, 7, 70, 1064, 22771, 628845, 20987848, 810313190, ...];
n=7: [1, 8, 84, 1288, 27426, 751752, 24925092, 957323320, ...]; ...
Compare to the table of coefficients in (1 + x*A(x)^(2*n+3))^(n+1):
n=0: [1, 1,   3,   24,   358,   7896,   225396,   7727644, ...];
n=1: [1, 2,  11,  100,  1465,  31070,   859367,  28808972, ...];
n=2: [1, 3,  24,  253,  3780,  77994,  2089024,  68277867, ...];
n=3: [1, 4,  42,  508,  7915, 161316,  4196916, 133476480, ...];
n=4: [1, 5,  65,  890, 14635, 298981,  7602705, 235213110, ...];
n=5: [1, 6,  93, 1424, 24858, 515352, 12914214, 389369448, ...];
n=6: [1, 7, 126, 2135, 39655, 842331, 20987848, 619044602, ...];
n=7: [1, 8, 164, 3048, 60250, 1320480, 32998388, 957323320, ...]; ...
to see that the main diagonals of the tables are the same.
		

Crossrefs

Programs

  • PARI
    {a(n) = my(A=[1]); for(m=1, n, A=concat(A, 0); A[m+1] = (Vec((1+x*Ser(A)^(2*m+3))^(m+1))[m+1] - Vec(Ser(A)^(m+1))[m+1])/(m+1) ); A[n+1]}
    for(n=0, 20, print1(a(n), ", "))

Formula

G.f. A(x) = Sum_{n>=0} a(n)*x^n satisfies:
(1) [x^n] A(x)^(n+1) = [x^n] (1 + x*A(x)^(2*n+3))^(n+1) for n>=0.
(2) A(x) = Sum_{n>=0} b(n) * x^n/A(x)^n, where b(n) = [x^n] (1 + x*A(x)^(2*n+3))^(n+1) / (n+1).
a(n) ~ c * d^n * n! * n^alpha, where d = 3.93464558322824528799..., alpha = 2.672403832022..., c = 0.0085435225111... - Vaclav Kotesovec, Feb 06 2023

A360342 G.f. A(x) satisfies: [x^n] A(x)^(n+1) = [x^n] (1 + x*A(x)^(2*n-2))^(n+1) for n >= 0.

Original entry on oeis.org

1, 1, 2, 20, 316, 6686, 173379, 5255624, 180911070, 6938866748, 292678301988, 13446616806957, 668017569348751, 35678261176871802, 2038906890461704040, 124171127134721710130, 8030684434410398312840, 549848454475826567644385, 39744302449387229743134043
Offset: 0

Views

Author

Paul D. Hanna, Feb 05 2023

Keywords

Examples

			G.f.: A(x) = 1 + x + 2*x^2 + 20*x^3 + 316*x^4 + 6686*x^5 + 173379*x^6 + 5255624*x^7 + 180911070*x^8 + 6938866748*x^9 + ...
RELATED SERIES.
G.f. A(x) = B(x/A(x)) where B(x) = B(x*A(x)) begins:
B(x) = 1 + x + 3*x^2 + 27*x^3 + 417*x^4 + 8727*x^5 + 225018*x^6 + 6800714*x^7 + 233778499*x^8 + ... + b(n)*x^n + ...
such that b(n) = [x^n] (1 + x*A(x)^(2*n-2))^(n+1) / (n+1),
as well as b(n) = [x^n] A(x)^(n+1) / (n+1),
so that b(n) begin:
[1/1, 2/2, 9/3, 108/4, 2085/5, 52362/6, 1575126/7, 54405712/8, ...].
ILLUSTRATION OF DEFINITION.
The table of coefficients of x^k in A(x)^(n+1) begins:
n=0: [1, 1,  2,  20,  316,  6686,  173379,  5255624, ...];
n=1: [1, 2,  5,  44,  676, 14084,  361794, 10897390, ...];
n=2: [1, 3,  9,  73, 1086, 22266,  566441, 16950588, ...];
n=3: [1, 4, 14, 108, 1553, 31312,  788620, 23442284, ...];
n=4: [1, 5, 20, 150, 2085, 41311, 1029745, 30401460, ...];
n=5: [1, 6, 27, 200, 2691, 52362, 1291355, 37859166, ...];
n=6: [1, 7, 35, 259, 3381, 64575, 1575126, 45848685, ...];
n=7: [1, 8, 44, 328, 4166, 78072, 1882884, 54405712, ...]; ...
Compare to the table of coefficients in (1 + x*A(x)^(2*n-2))^(n+1):
n=0: [1, 1,  -2,   -1,   -32,   -519,  -11490,  -305967, ...];
n=1: [1, 2,   1,    0,     0,      0,       0,        0, ...];
n=2: [1, 3,   9,   28,   180,   2379,   47111,  1182009, ...];
n=3: [1, 4,  22,  108,   745,   8556,  153292,  3658316, ...];
n=4: [1, 5,  40,  265,  2085,  22706,  366450,  8157230, ...];
n=5: [1, 6,  63,  524,  4743,  52362,  781973, 16041192, ...];
n=6: [1, 7,  91,  910,  9415, 109536, 1575126, 29886445, ...];
n=7: [1, 8, 124, 1448, 16950, 211840, 3042820, 54405712, ...]; ...
to see that the main diagonals of the tables are the same.
		

Crossrefs

Programs

  • PARI
    {a(n) = my(A=[1]); for(m=1, n, A=concat(A, 0); A[m+1] = (Vec((1+x*Ser(A)^(2*m-2))^(m+1))[m+1] - Vec(Ser(A)^(m+1))[m+1])/(m+1) ); A[n+1]}
    for(n=0, 20, print1(a(n), ", "))

Formula

G.f. A(x) = Sum_{n>=0} a(n)*x^n satisfies:
(1) [x^n] A(x)^(n+1) = [x^n] (1 + x*A(x)^(2*n-2))^(n+1) for n>=0.
(2) A(x) = Sum_{n>=0} b(n) * x^n/A(x)^n, where b(n) = [x^n] (1 + x*A(x)^(2*n-2))^(n+1) / (n+1).
a(n) ~ c * d^n * n! * n^alpha, where d = 3.93464558322824528799..., alpha = 0.0798993252137..., c = 0.118957192149397... - Vaclav Kotesovec, Feb 06 2023

A360343 G.f. A(x) satisfies: [x^n] A(x)^(n+1) = [x^n] (1 + x*A(x)^(2*n-1))^(n+1) for n >= 0.

Original entry on oeis.org

1, 1, 3, 31, 526, 11907, 328980, 10580531, 384937042, 15549217485, 688430225102, 33096289502982, 1715499922758709, 95339852384471586, 5655337634718941111, 356683962066445400017, 23840465113068534382248, 1683771696557415075462436, 125327912444852044066759399
Offset: 0

Views

Author

Paul D. Hanna, Feb 05 2023

Keywords

Examples

			G.f.: A(x) = 1 + x + 3*x^2 + 31*x^3 + 526*x^4 + 11907*x^5 + 328980*x^6 + 10580531*x^7 + 384937042*x^8 + 15549217485*x^9 + ...
RELATED SERIES.
G.f. A(x) = B(x/A(x)) where B(x) = B(x*A(x)) begins:
B(x) = 1 + x + 4*x^2 + 41*x^3 + 687*x^4 + 15433*x^5 + 424524*x^6 + 13620842*x^7 + 495005025*x^8 + ... + b(n)*x^n + ...
such that b(n) = [x^n] (1 + x*A(x)^(2*n-1))^(n+1) / (n+1),
as well as b(n) = [x^n] A(x)^(n+1) / (n+1),
so that {b(n)} begins:
[1/1, 2/2, 12/3, 164/4, 3435/5, 92598/6, 2971668/7, 108966736/8, ...].
ILLUSTRATION OF DEFINITION.
The table of coefficients of x^k in A(x)^(n+1) begins:
  n=0: [1, 1,  3,  31,  526,  11907,  328980,  10580531, ...];
  n=1: [1, 2,  7,  68, 1123,  25052,  685891,  21923076, ...];
  n=2: [1, 3, 12, 112, 1800,  39555, 1072896,  34076544, ...];
  n=3: [1, 4, 18, 164, 2567,  55548, 1492336,  47093172, ...];
  n=4: [1, 5, 25, 225, 3435,  73176, 1946745,  61028770, ...];
  n=5: [1, 6, 33, 296, 4416,  92598, 2438866,  75942984, ...];
  n=6: [1, 7, 42, 378, 5523, 113988, 2971668,  91899578, ...];
  n=7: [1, 8, 52, 472, 6770, 137536, 3548364, 108966736, ...]; ...
Compare to the table of coefficients in (1 + x*A(x)^(2*n-1))^(n+1):
  n=0: [1, 1,  -1,   -2,   -26,   -463,  -10778,   -303048, ...];
  n=1: [1, 2,   3,    8,    69,   1120,   24937,    683012, ...];
  n=2: [1, 3,  12,   55,   444,   6351,  132492,   3504654, ...];
  n=3: [1, 4,  26,  164,  1411,  18560,  357624,   9024812, ...];
  n=4: [1, 5,  45,  360,  3435,  43926,  785715,  18700710, ...];
  n=5: [1, 6,  69,  668,  7134,  92598, 1570420,  35086104, ...];
  n=6: [1, 7,  98, 1113, 13279, 179816, 2971668,  62645353, ...];
  n=7: [1, 8, 132, 1720, 22794, 327032, 5403036, 108966736, ...]; ...
to see that the main diagonals of the tables are the same.
		

Crossrefs

Programs

  • PARI
    {a(n) = my(A=[1]); for(m=1, n, A=concat(A, 0); A[m+1] = (Vec((1+x*Ser(A)^(2*m-1))^(m+1))[m+1] - Vec(Ser(A)^(m+1))[m+1])/(m+1) ); A[n+1]}
    for(n=0, 20, print1(a(n), ", "))

Formula

G.f. A(x) = Sum_{n>=0} a(n)*x^n satisfies:
(1) [x^n] A(x)^(n+1) = [x^n] (1 + x*A(x)^(2*n-1))^(n+1) for n>=0.
(2) A(x) = Sum_{n>=0} b(n) * x^n/A(x)^n, where b(n) = [x^n] (1 + x*A(x)^(2*n-1))^(n+1) / (n+1).
a(n) ~ c * d^n * n! * n^alpha, where d = 3.93464558322824528799..., alpha = 0.5984002265754..., c = 0.08321697608093... - Vaclav Kotesovec, Feb 06 2023

A360336 G.f. A(x) satisfies: [x^n] A(x)^(n+1) = [x^n] (1 + x*A(x)^(3*n))^(n+1) for n >= 0.

Original entry on oeis.org

1, 1, 6, 99, 2608, 90800, 3835458, 187727106, 10356030404, 632391914502, 42217751766193, 3053486035335835, 237640678130730437, 19794116975373467259, 1756875217029906875379, 165552614838271944281933, 16509692094523556884973416, 1737510282985845400007263814
Offset: 0

Views

Author

Paul D. Hanna, Feb 06 2023

Keywords

Examples

			G.f.: A(x) = 1 + x + 6*x^2 + 99*x^3 + 2608*x^4 + 90800*x^5 + 3835458*x^6 + 187727106*x^7 + 10356030404*x^8 + 632391914502*x^9 + ...
RELATED SERIES.
G.f. A(x) = B(x/A(x)) where B(x) = B(x*A(x)) begins:
B(x) = 1 + x + 7*x^2 + 118*x^3 + 3113*x^4 + 108221*x^5 + 4564720*x^6 + 223208259*x^7 + 12307249017*x^8 + ... + b(n)*x^n + ...
such that b(n) = [x^n] (1 + x*A(x)^(3*n))^(n+1) / (n+1),
as well as b(n) = [x^n] A(x)^(n+1) / (n+1),
so that b(n) begin:
[1/1, 2/2, 21/3, 472/4, 15565/5, 649326/6, 31953040/7, 1785666072/8, ...].
ILLUSTRATION OF DEFINITION.
The table of coefficients of x^k in A(x)^(n+1) begins:
n=0: [1, 1,  6,   99,  2608,  90800,  3835458,  187727106, ...];
n=1: [1, 2, 13,  210,  5450, 188004,  7893613,  384731112, ...];
n=2: [1, 3, 21,  334,  8544, 292017, 12186069,  591418401, ...];
n=3: [1, 4, 30,  472, 11909, 403268, 16725042,  808213780, ...];
n=4: [1, 5, 40,  625, 15565, 522211, 21523390, 1035561335, ...];
n=5: [1, 6, 51,  794, 19533, 649326, 26594644, 1273925322, ...];
n=6: [1, 7, 63,  980, 23835, 785120, 31953040, 1523791095, ...];
n=7: [1, 8, 76, 1184, 28494, 930128, 37613552, 1785666072, ...]; ...
Compare to the table of coefficients in (1 + x*A(x)^(3*n))^(n+1):
n=0: [1, 1,   0,    0,     0,       0,        0,          0, ...];
n=1: [1, 2,   7,   48,   719,   17882,   603567,   25021464, ...];
n=2: [1, 3,  21,  190,  2814,   65460,  2105997,   84726534, ...];
n=3: [1, 4,  42,  472,  7303,  162828,  4982706,  193437168, ...];
n=4: [1, 5,  70,  940, 15565,  341796, 10002300,  373126910, ...];
n=5: [1, 6, 105, 1640, 29340,  649326, 18377374,  658075230, ...];
n=6: [1, 7, 147, 2618, 50729, 1150968, 31953040, 1101647800, ...];
n=7: [1, 8, 196, 3920, 82194, 1934296, 53433184, 1785666072, ...]; ...
to see that the main diagonals of the tables are the same.
		

Crossrefs

Programs

  • PARI
    {a(n) = my(A=[1]); for(m=1, n, A=concat(A, 0); A[m+1] = (Vec((1+x*Ser(A)^(3*m))^(m+1))[m+1] - Vec(Ser(A)^(m+1))[m+1])/(m+1) ); A[n+1]}
    for(n=0, 20, print1(a(n), ", "))

Formula

G.f. A(x) = Sum_{n>=0} a(n)*x^n satisfies:
(1) [x^n] A(x)^(n+1) = [x^n] (1 + x*A(x)^(3*n))^(n+1) for n>=0.
(2) A(x) = Sum_{n>=0} b(n) * x^n/A(x)^n, where b(n) = [x^n] (1 + x*A(x)^(3*n))^(n+1) / (n+1).
a(n) ~ c * d^n * n! * n^alpha, where d = 5.7189630165873859806..., alpha = 1.0541176773983..., c = 0.03951220887392... - Vaclav Kotesovec, Feb 06 2023
Showing 1-6 of 6 results.