The OEIS is supported by the many generous donors to the OEIS Foundation.

 Hints (Greetings from The On-Line Encyclopedia of Integer Sequences!)
 A357793 a(n) = coefficient of x^n in A(x) = Sum_{n>=0} x^n*F(x)^n * (1 - x^n*F(x)^n)^n, where F(x) = 1 + x*F(x)^3 is a g.f. of A001764. 2

%I #15 Mar 14 2023 04:56:55

%S 1,1,1,4,14,64,314,1633,8826,49107,279349,1617290,9498099,56445918,

%T 338817460,2051182532,12509647159,76785827812,474000090118,

%U 2940761033970,18327028477625,114677403429121,720191795608082,4537925593859911,28679991910774479,181761824439041725

%N a(n) = coefficient of x^n in A(x) = Sum_{n>=0} x^n*F(x)^n * (1 - x^n*F(x)^n)^n, where F(x) = 1 + x*F(x)^3 is a g.f. of A001764.

%C Related identity: 0 = Sum_{n=-oo..+oo} x^n * (1 - x^n*F(x)^n)^n / (1 - x*F(x)^2)^n, where F(x) = 1 + x*F(x)^3.

%H Paul D. Hanna, <a href="/A357793/b357793.txt">Table of n, a(n) for n = 0..400</a>

%F Given F(x) = 1 + x*F(x)^3, g.f. A(x) = Sum_{n>=0} a(n)*x^n may be defined by the following.

%F (1) A(x) = Sum_{n>=0} x^n * F(x)^n * (1 - x^n*F(x)^n)^n.

%F (2) A(x) = Sum_{n>=1} (-1)^(n-1) * (x*F(x))^(n*(n-1)) / (1 - x^n*F(x)^n)^n.

%F (3) A(x) = Sum_{n>=0} x^n * (1 - x^n*F(x)^n)^n / (1 - x*F(x)^2)^n.

%F (4) A(x) = Sum_{n>=1} (-1)^(n-1) * x^(n*(n-1)) * F(x)^(n^2) * (1 - x*F(x)^2)^n / (1 - x^n*F(x)^n)^n.

%F a(n) ~ c * 3^(3*n) / (n^(3/2) * 2^(2*n)), where c = 0.0403028056146458801802487899052088995113692232406693619.... - _Vaclav Kotesovec_, Mar 14 2023

%e G.f.: A(x) = 1 + x + x^2 + 4*x^3 + 14*x^4 + 64*x^5 + 314*x^6 + 1633*x^7 + 8826*x^8 + 49107*x^9 + 279349*x^10 + 1617290*x^11 + 9498099*x^12 + ...

%e where

%e F(x) = 1 + x*F(x)*(1 - x*F(x)) + x^2*F(x)^2*(1 - x^2*F(x)^2) + x^3*F(x)^3*(1 - x^3*F(x)^3) + x^4*F(x)^4*(1 - x^4*F(x)^4) + ... + x^n * F(x)^n * (1 - x^n*F(x)^n)^n + ...

%e also,

%e F(x) = 1/(1 - x*F(x)) - (x*F(x))^2/(1 - x^2*F(x)^2)^2 + (x*F(x))^6/(1 - x^3*F(x)^3)^3 - (x*F(x))^12/(1 - x^4*F(x)^4)^4 + (x*F(x))^20/(1 - x^5*F(x)^4)^5 +- ... + (-1)^(n-1) * (x*F(x))^(n*(n-1)) / (1 - x^n*F(x)^n)^n + ...

%e Where F(x) = 1 + x*F(x)^3 begins

%e F(x) = 1 + x + 3*x^2 + 12*x^3 + 55*x^4 + 273*x^5 + 1428*x^6 + 7752*x^7 + 43263*x^8 + 246675*x^9 + 1430715*x^10 + ... + A001764(n)*x^n + ...

%e SPECIFIC VALUES.

%e The radius of convergence of the power series A(x) equals 4/27.

%e The power series A(x) converges at x = 4/27 to

%e A(4/27) = 1.2311920996301390036800654138630946234233891541082821783156...

%e which equals the following sums:

%e (1) A(4/27) = Sum_{n>=0} 2^n * (9^n - 2^n)^n / 9^(n*(n+1)),

%e (2) A(4/27) = Sum_{n>=1} (-1)^(n-1) * 2^(n*(n-1)) * 9^n / (9^n - 2^n)^n.

%o (PARI) {a(n) = my(A=1, F = (serreverse(x/(1+x + O(x^(n+2)))^3)/x)^(1/3));

%o A = sum(m=0,n, x^m * F^m * (1 - x^m*F^m)^m); polcoeff(A,n)}

%o for(n=0,30, print1(a(n),", "))

%o (PARI) {a(n) = my(A=1, F = (serreverse(x/(1+x + O(x^(n+3)))^3)/x)^(1/3));

%o A = sum(m=1,n+1, (-1)^(m-1) * (x*F)^(m*(m-1)) / (1 - x^m*F^m)^m); polcoeff(A,n)}

%o for(n=0,30, print1(a(n),", "))

%Y Cf. A357792, A001764.

%K nonn

%O 0,4

%A _Paul D. Hanna_, Dec 20 2022

Lookup | Welcome | Wiki | Register | Music | Plot 2 | Demos | Index | Browse | More | WebCam
Contribute new seq. or comment | Format | Style Sheet | Transforms | Superseeker | Recents
The OEIS Community | Maintained by The OEIS Foundation Inc.

Last modified September 15 03:41 EDT 2024. Contains 375931 sequences. (Running on oeis4.)