login
A323688
G.f.: Sum_{n>=0} 3^n * ((1+x)^n + i)^n / (4 + 3*i*(1+x)^n)^(n+1), where i^2 = -1.
2
1, 39, 4914, 1036674, 306362988, 116429448762, 54085355493876, 29694009536806140, 18811303988392369116, 13506257834309548172844, 10838311025576842625186286, 9612978973113439745046522234, 9338259937723855837235460233514, 9860260989532664468069402304530610, 11244403367868838452318619557009699582, 13772731580491669420660957058197302913122, 18033025522536379279790808472409496887914884
OFFSET
0,2
COMMENTS
It is remarkable that the generating function results in a power series in x with only real coefficients.
FORMULA
G.f.: Sum_{n>=0} 3^n * ((1+x)^n + i)^n / (4 + 3*i*(1+x)^n)^(n+1).
G.f.: Sum_{n>=0} 3^n * ((1+x)^n - i)^n / (4 - 3*i*(1+x)^n)^(n+1).
G.f.: Sum_{n>=0} 3^n * ((1+x)^n + i)^n * (4 - 3*i*(1+x)^n)^(n+1) / (16 + 9*(1+x)^(2*n))^(n+1).
G.f.: Sum_{n>=0} 3^n * ((1+x)^n - i)^n * (4 + 3*i*(1+x)^n)^(n+1) / (16 + 9*(1+x)^(2*n))^(n+1).
EXAMPLE
G.f.: A(x) = 1 + 39*x + 4914*x^2 + 1036674*x^3 + 306362988*x^4 + 116429448762*x^5 + 54085355493876*x^6 + 29694009536806140*x^7 + 18811303988392369116*x^8 + ...
such that
A(x) = 1/(4+3*i) + 3*((1+x) + i)/(4 + 3*i*(1+x))^2 + 3^2*((1+x)^2 + i)^2/(4 + 3*i*(1+x)^2)^3 + 3^3*((1+x)^3 + i)^3/(4 + 3*i*(1+x)^3)^4 + 3^4*((1+x)^4 + i)^4/(4 + 3*i*(1+x)^4)^5 + 3^5*((1+x)^5 + i)^5/(4 + 3*i*(1+x)^5)^6 + 3^6*((1+x)^6 + i)^6/(4 + 3*i*(1+x)^6)^7 + ...
also
A(x) = 1/(4-3*i) + 3*((1+x) - i)/(4 - 3*i*(1+x))^2 + 3^2*((1+x)^2 - i)^2/(4 - 3*i*(1+x)^2)^3 + 3^3*((1+x)^3 - i)^3/(4 - 3*i*(1+x)^3)^4 + 3^4*((1+x)^4 - i)^4/(4 - 3*i*(1+x)^4)^5 + 3^5*((1+x)^5 - i)^5/(4 - 3*i*(1+x)^5)^6 + 3^6*((1+x)^6 - i)^6/(4 - 3*i*(1+x)^6)^7 + ...
RELATED INFINITE SERIES.
At x = -1/4, the g.f. as a power series in x diverges, but the related series converges:
S = Sum_{n>=0} 3^n * ((3/4)^n + i)^n / (4 + 3*i*(3/4)^n)^(n+1).
Equivalently,
S = Sum_{n>=0} 12^n * (3^n + 4^n*i)^n / (4^(n+1) + 3^(n+1)*i)^(n+1) ;
written explicitly,
S = 1/(4+3*i) + 12*(3+4*i)/(4^2+3^2*i)^2 + 12^2*(3^2+4^2*i)^2/(4^3+3^3*i)^3
+ 12^3*(3^3+4^3*i)^3/(4^4+3^4*i)^4 + 12^4*(3^4+4^4*i)^4/(4^5+3^5*i)^5
+ 12^5*(3^5+4^5*i)^5/(4^6+3^6*i)^6 + 12^6*(3^6+4^6*i)^6/(4^7+3^7*i)^7 + ...
which equals the real number
S = 0.33582962080420548450017336765358292462921076160047749762227134638...
PROG
(PARI) {a(n) = my(A = sum(m=0, n*60 + 600, 3^m*((1+x +x*O(x^n))^m + I)^m/(4 + 3*I*(1+x +x*O(x^n))^m)^(m+1)*1. )); round(polcoeff(A, n))}
for(n=0, 20, print1(a(n), ", "))
(PARI) {a(n) = my(A = sum(m=0, n*60 + 600, 3^m*((1+x +x*O(x^n))^m - I)^m/(4 - 3*I*(1+x +x*O(x^n))^m)^(m+1)*1. )); round(polcoeff(A, n))}
for(n=0, 20, print1(a(n), ", "))
CROSSREFS
Sequence in context: A176073 A145619 A027490 * A266453 A279116 A266307
KEYWORD
nonn
AUTHOR
Paul D. Hanna, Feb 12 2019
STATUS
approved