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”).
%I #16 Aug 31 2014 01:49:23
%S 1,4,24,168,1286,10440,88112,764368,6766278,60828024,553529808,
%T 5086837680,47127896444,439608960656,4124536224864,38891699480992,
%U 368326082421446,3501654020899800,33403335855108368,319612386771594608,3066480362268978804,29493401582426082032,284301304326376855200
%N G.f.: Sum_{n>=0} x^n / (1-3*x)^(2*n+1) * [Sum_{k=0..n} C(n,k)^2 * 2^k * x^k] * [Sum_{k=0..n} C(n,k)^2 * 3^k * x^k].
%C a(n) == 2 (mod 4) iff n = 2^k for k>=2, and a(n) == 0 (mod 4) elsewhere except at a(0)=1 (conjecture).
%H Paul D. Hanna, <a href="/A246423/b246423.txt">Table of n, a(n) for n = 0..300</a>
%F G.f.: Sum_{n>=0} x^n / (1-x)^(2*n+1) * [Sum_{k=0..n} C(n,k)^2 * x^k] * [Sum_{k=0..n} C(n,k)^2 * 3^(n-k) * 2^k * x^k].
%F G.f.: Sum_{n>=0} x^n * Sum_{k=0..n} C(n,k)^2 * Sum_{j=0..k} C(k,j)^2 * 3^(k-j) * 2^j * x^j.
%F G.f.: Sum_{n>=0} x^n * Sum_{k=0..n} C(n,k)^2 * 3^(n-k) * Sum_{j=0..k} C(k,j)^2 * 2^j * x^j.
%F a(n) = Sum_{k=0..[n/2]} 2^k * Sum_{j=0..n-2*k} C(n-k, k+j)^2 * C(k+j, j)^2 * 3^j.
%F Recurrence: (n-5)*(n-4)*(n-2)*n^2*a(n) = 4*(n-5)*(n-4)*(4*n^3 - 12*n^2 + 10*n - 3)*a(n-1) - 16*(n-5)*(n-4)*(n-1)*(2*n-3)^2*a(n-2) + 8*(n-5)*(n-2)*(4*n^3 - 28*n^2 + 54*n - 27)*a(n-3) + 24*(n-3)*(5*n^4 - 60*n^3 + 248*n^2 - 408*n + 216)*a(n-4) + 16*(n-4)*(n-1)*(4*n^3 - 44*n^2 + 150*n - 153)*a(n-5) - 64*(n-5)*(n-2)*(n-1)*(2*n-9)^2*a(n-6) + 32*(n-2)*(n-1)*(4*n^3 - 60*n^2 + 298*n - 489)*a(n-7) - 16*(n-6)^2*(n-4)*(n-2)*(n-1)*a(n-8). - _Vaclav Kotesovec_, Aug 26 2014
%F a(n) ~ c * d^n / n, where d = 10.094399065494857710014687346... is the root of the equation 16 - 128*d + 256*d^2 - 64*d^3 - 120*d^4 - 32*d^5 + 64*d^6 - 16*d^7 + d^8 = 0, and c = 0.5132545324612697424702223429844481717... . - _Vaclav Kotesovec_, Aug 26 2014
%e G.f.: A(x) = 1 + 4*x + 24*x^2 + 168*x^3 + 1286*x^4 + 10440*x^5 +...
%e where the g.f. is given by the binomial series identity:
%e A(x) = 1/(1-3*x) + x/(1-3*x)^3 * (1 + 2*x) * (1 + 3*x)
%e + x^2/(1-3*x)^5 * (1 + 2^2*2*x + 4*x^2) * (1 + 2^2*3*x + 9*x^2)
%e + x^3/(1-3*x)^7 * (1 + 3^2*2*x + 3^2*4*x^2 + 8*x^3) * (1 + 3^2*3*x + 3^2*9*x^2 + 27*x^3)
%e + x^4/(1-3*x)^9 * (1 + 4^2*2*x + 6^2*4*x^2 + 4^2*8*x^3 + 16*x^4) * (1 + 4^2*3*x + 6^2*9*x^2 + 4^2*27*x^3 + 81*x^4)
%e + x^5/(1-3*x)^11 * (1 + 5^2*2*x + 10^2*4*x^2 + 10^2*8*x^3 + 5^2*16*x^4 + 32*x^5) * (1 + 5^2*3*x + 10^2*9*x^2 + 10^2*27*x^3 + 5^2*81*x^4 + 243*x^5) +...
%e equals the series
%e A(x) = 1/(1-x) + x/(1-x)^3 * (1 + x) * (3+2*x)
%e + x^2/(1-x)^5 * (1 + 2^2*x + x^2) * (9+2^2*3*2*x+4*x^2)
%e + x^3/(1-x)^7 * (1 + 3^2*x + 3^2*x^2 + x^3) * (27+3^2*9*2*x+3^2*4*3*x^2+8*x^3)
%e + x^4/(1-x)^9 * (1 + 4^2*x + 6^2*x^2 + 4^2*x^3 + x^4) * (81+4^2*27*2*x+6^2*9*4*x^2+4^2*3*8*x^3+16*x^4)
%e + x^5/(1-x)^11 * (1 + 5^2*x + 10^2*x^2 + 10^2*x^3 + 5^2*x^4 + x^5) * (243+5^2*81*2*x+10^2*27*4*x^2+10^2*9*16*x^3+5^2*3*18*x^4+32*x^5) +...
%e We can also express the g.f. by another binomial series identity:
%e A(x) = 1 + x*(3 + (1+2*x)) + x^2*(9 + 2^2*3*(1+2*x) + (1+2^2*2*x+4*x^2))
%e + x^3*(27 + 3^2*9*(1+2*x) + 3^2*3*(1+2^2*2*x+4*x^2) + (1+3^2*2*x+3^2*4*x^2+8*x^3))
%e + x^4*(81 + 4^2*27*(1+2*x) + 6^2*9*(1+2^2*2*x+4*x^2) + 4^2*3*(1+3^2*2*x+3^2*4*x^2+8*x^3) + (1+4^2*2*x+6^2*4*x^2+4^2*8*x^3+16*x^4))
%e + x^5*(243 + 5^2*81*(1+2*x) + 10^2*27*(1+2^2*2*x+4*x^2) + 10^2*9*(1+3^2*2*x+3^2*4*x^2+8*x^3) + 5^2*3*(1+4^2*2*x+6^2*4*x^2+4^2*8*x^3+16*x^4) + (1+5^2*2*x+10^2*4*x^2+10^2*8*x^3+5^2*16*x^4+32*x^5)) +...
%e equals the series
%e A(x) = 1 + x*(1 + (3+2*x)) + x^2*(1 + 2^2*(3+2*x) + (9+2^2*3*2*x+4*x^2))
%e + x^3*(1 + 3^2*(3+2*x) + 3^2*(9+2^2*3*2*x+4*x^2) + (27+3^2*9*2*x+3^2*4*3*x^2+8*x^3))
%e + x^4*(1 + 4^2*(3+2*x) + 6^2*(9+2^2*3*2*x+4*x^2) + 4^2*(27+3^2*9*2*x+3^2*4*3*x^2+8*x^3) + (81+4^2*27*2*x+6^2*9*4*x^2+4^2*3*8*x^3+16*x^4))
%e + x^5*(1 + 5^2*(3+2*x) + 10^2*(9+2^2*3*2*x+4*x^2) + 10^2*(27+3^2*9*2*x+3^2*4*3*x^2+8*x^3) + 5^2*(81+4^2*27*2*x+6^2*9*4*x^2+4^2*3*8*x^3+16*x^4) + (243+5^2*81*2*x+10^2*27*4*x^2+10^2*9*16*x^3+5^2*3*18*x^4+32*x^5)) +...
%o (PARI) /* By definition: */
%o {a(n)=local(A=1); A=sum(m=0, n, x^m/(1-3*x)^(2*m+1) * sum(k=0, m, binomial(m, k)^2 * 2^k * x^k) * sum(k=0, m, binomial(m, k)^2 * 3^k * x^k) +x*O(x^n)); polcoeff(A, n)}
%o for(n=0, 25, print1(a(n), ", "))
%o (PARI) /* By a binomial identity: */
%o {a(n)=local(A=1); A=sum(m=0, n, x^m/(1-x)^(2*m+1) * sum(k=0, m, binomial(m, k)^2 * 3^(m-k) * 2^k * x^k) * sum(k=0, m, binomial(m, k)^2 * x^k) +x*O(x^n)); polcoeff(A, n)}
%o for(n=0, 25, print1(a(n), ", "))
%o (PARI) /* By a binomial identity: */
%o {a(n)=polcoeff(sum(m=0, n, x^m*sum(k=0, m, binomial(m, k)^2 * 3^(m-k) * sum(j=0, k, binomial(k, j)^2 * 2^j * x^j)+x*O(x^n))), n)}
%o for(n=0, 25, print1(a(n), ", "))
%o (PARI) /* By a binomial identity: */
%o {a(n)=polcoeff(sum(m=0, n, x^m*sum(k=0, m, binomial(m, k)^2 * sum(j=0, k, binomial(k, j)^2 * 3^(k-j) * 2^j * x^j)+x*O(x^n))), n)}
%o for(n=0, 25, print1(a(n), ", "))
%o (PARI) /* Formula for a(n): */
%o {a(n)=sum(k=0, n\2, sum(j=0, n-2*k, 2^k * binomial(n-k, k+j)^2 * binomial(k+j, j)^2 * 3^j))}
%o for(n=0, 25, print1(a(n), ", "))
%Y Cf. A246056 (dual), A243948, A245929, A227845, A245925.
%K nonn
%O 0,2
%A _Paul D. Hanna_, Aug 25 2014