login

Year-end appeal: Please make a donation to the OEIS Foundation to support ongoing development and maintenance of the OEIS. We are now in our 61st year, we have over 378,000 sequences, and we’ve reached 11,000 citations (which often say “discovered thanks to the OEIS”).

A307811
Expansion of 1/AGM(1-45*x, sqrt((1-25*x)*(1-81*x))).
2
1, 49, 2601, 148225, 8970025, 570111129, 37678303881, 2567836387809, 179267329355625, 12754414737348025, 921185098227422161, 67340346156989933769, 4971327735657992896201, 369994703739586257235225, 27725052308247030792515625, 2089567204521186409129541025
OFFSET
0,2
COMMENTS
See A246923.
Also the squares of coefficients in g.f. 1/sqrt((1-5*x)*(1-9*x)).
LINKS
FORMULA
a(n) = A104454(n)^2 = (Sum_{k=0..n} 5^(n-k)*binomial(n,k)*binomial(2k,k))^2 = (Sum_{k=0..n} 9^(n-k)*(-1)^k*binomial(n,k)*binomial(2k,k))^2.
a(n) ~ 3^(4*n+2) / (4*Pi*n). - Vaclav Kotesovec, Sep 27 2019
MATHEMATICA
a[n_] := Sum[5^(n-k) * Binomial[n, k] * Binomial[2*k, k], {k, 0, n}]^2; Array[a, 16, 0] // Flatten (* Amiram Eldar, May 13 2021 *)
PROG
(PARI) N=20; x='x+O('x^N); Vec(1/agm(1-45*x, sqrt((1-25*x)*(1-81*x))))
(PARI) {a(n) = sum(k=0, n, 5^(n-k)*binomial(n, k)*binomial(2*k, k))^2}
(PARI) {a(n) = sum(k=0, n, 9^(n-k)*(-1)^k*binomial(n, k)*binomial(2*k, k))^2}
CROSSREFS
Cf. A104454.
(Sum_{k=0..n} c^(n-k)*e^k*binomial(n,k)*binomial(2k,k))^2 = (Sum_{k=0..n} d^(n-k)*(-e)^k*binomial(n,k)*binomial(2k,k))^2, where e = (d-c)/4: A002894 (c=0,d=4,e=1), A246467 (c=1,d=5,e=1), A246876 (c=2,d=6,e=1), A246906 (c=3,d=7,e=1), this sequence (c=5,d=9,e=1), A322240 (c=-3,d=5,e=2), A322243 (c=-1,d=7,e=2), A246923 (c=1,d=9,e=2), A248167 (c=3, d=11,e=2), A322247 (c=-1,d=11,e=3), A307810 (c=4,d=16,e=3), A322245 (c=-5,d=11,e=4), A322249 (c=-3,d=13,e=4).
Sequence in context: A069741 A203384 A099367 * A123841 A014773 A376092
KEYWORD
nonn
AUTHOR
Seiichi Manyama, Apr 30 2019
STATUS
approved