login
A360339
a(n) = coefficient of x^n*y^(2*n+1)/n! in log( Sum_{n>=0} (n + y)^(3*n) * x^n/n! ).
3
1, 6, 99, 2832, 117405, 6423408, 438143391, 35869775616, 3430351996569, 375544727136000, 46333978359977979, 6362713275564589056, 962689133095843525749, 159139760744994666835968, 28539360163037720058960375, 5518961894002049077780611072, 1144859158421455331276272257201
OFFSET
1,2
LINKS
FORMULA
E.g.f. A(x) = Sum_{n>=1} a(n)*x^n/n! may be defined as follows.
(1) A(x) = Limit_{N->oo} (1/N) * log( Sum_{n>=0} (N + n)^(3*n) * (x/N^2)^n/n! ).
(2) a(n) = [x^n*y^(2*n+1)/n!] log( Sum_{n>=0} (n + y)^(3*n) * x^n/n! ).
a(n) ~ c * d^n * n! / n^(5/2), where d = (3/2) * (3 + sqrt(6)) * exp(3 - sqrt(6)) = 14.175247991325192557234088913125084764719990898660219459... and c = sqrt((3 - sqrt(6))/Pi)/4 = 0.1046520596183180437324097699670683850916674939335504... - Vaclav Kotesovec, Feb 12 2023, updated Mar 20 2024
EXAMPLE
E.g.f.: A(x) = x + 6*x^2/2! + 99*x^3/3! + 2832*x^4/4! + 117405*x^5/5! + 6423408*x^6/6! + 438143391*x^7/7! + 35869775616*x^8/8! + ... + a(n)*x^n/n! + ...
where a(n) equals the coefficient of y^(2*n+1)*x^n/n! in the series given by
log( Sum_{n>=0} (n + y)^(3*n) * x^n/n! ) = (y^3 + 3*y^2 + 3*y + 1)*x + (6*y^5 + 45*y^4 + 140*y^3 + 225*y^2 + 186*y + 63)*x^2/2! + (99*y^7 + 1305*y^6 + 7722*y^5 + 26514*y^4 + 56844*y^3 + 75780*y^2 + 57915*y + 19493)*x^3/3! + (2832*y^9 + 56214*y^8 + 521784*y^7 + 2965716*y^6 + 11339280*y^5 + 30131946*y^4 + 55424512*y^3 + 67771380*y^2 + 49792368*y + 16686958)*x^4/4! + (117405*y^11 + 3214647*y^10 + 42201705*y^9 + 349928235*y^8 + 2030468625*y^7 + 8627152275*y^6 + 27284511927*y^5 + 63980788365*y^4 + 108602299435*y^3 + 126629082945*y^2 + 90978438315*y + 30421607649)*x^5/5! + ...
Exponentiation yields the e.g.f. of A266482:
exp(A(x)) = 1 + x + 7*x^2/2! + 118*x^3/3! + 3373*x^4/4! + 139096*x^5/5! + 7565779*x^6/6! + 513277024*x^7/7! + 41820455065*x^8/8! + 3982842285184*x^9/9! + 434457816912991*x^10/10! + ... + A266483(n)*x^n/n! + ...
which equals
lim_{N->oo} [ Sum_{n>=0} (N + n)^(3*n) * (x/N^2)^n/n! ]^(1/N).
RELATED SEQUENCES.
a(n) is divisible by n where a(n)/n begins:
[1, 3, 33, 708, 23481, 1070568, 62591913, 4483721952, ...].
PROG
(PARI) /* Using logarithmic formula */
{a(n) = n! * polcoeff( polcoeff( log( sum(m=0, n+1, (m + y)^(3*m) *x^m/m! ) +x*O(x^n) ), n, x), 2*n+1, y)}
for(n=1, 20, print1(a(n), ", "))
CROSSREFS
KEYWORD
nonn
AUTHOR
Paul D. Hanna, Feb 10 2023
STATUS
approved