|
%I #13 Nov 27 2017 12:04:01
%S 1,4,27,236,2375,26090,304241,3704860,46622655,602035556,7937288062,
%T 106451074614,1448267147717,19944962832826,277565209168861,
%U 3898075200816892,55182857681572655,786731161113510584
%N a(n) = Sum_{k=0..n} binomial(n+k, k)*binomial(n+k+1, k+1)/(n+1).
%C The function B(x,r) = x*Sum_{n>=0} b(n,r)*x^n, where
%C b(n,r) = Sum_{k=0..n} binomial(n+k, k)*binomial(n+k+1, r*k+1)/(n+1), satisfies
%C B(x/(1+x) - x^r, r) = x for all positive integer r except at r=1;
%C B(x,1)/x is the generating function of this sequence.
%H G. C. Greubel, <a href="/A185655/b185655.txt">Table of n, a(n) for n = 0..830</a>
%F Recurrence: 2*(n+1)^2*(2*n + 1)*(3*n - 2)*(7*n - 2)*a(n) = (1365*n^5 - 607*n^4 - 821*n^3 + 411*n^2 + 80*n - 44)*a(n-1) - 4*(n-2)*(2*n - 1)^2*(3*n + 1)*(7*n + 5)*a(n-2). - _Vaclav Kotesovec_, Nov 27 2017
%F a(n) ~ 2^(4*n+3) / (3*Pi*n^2). - _Vaclav Kotesovec_, Nov 27 2017
%e G.f.: A(x) = 1 + 4*x + 27*x^2 + 236*x^3 + 2375*x^4 + 26090*x^5 +...
%e Let G(x*A(x)) = x, then the series reversion of x*A(x) begins:
%e G(x) = x - 4*x^2 + 5*x^3 - 16*x^4 - 12*x^5 - 218*x^6 - 1197*x^7 - 8974*x^8 - 65582*x^9 - 503614*x^10 - 3956461*x^11 - ...
%e Does G(x) satisfy a nice functional equation?
%t Table[Sum[Binomial[n + k, k]*Binomial[n + k + 1, k + 1]/(n + 1), {k, 0, n}], {n, 0, 50}] (* _G. C. Greubel_, Jul 09 2017 *)
%o (PARI) {a(n)=sum(k=0, n, binomial(n+k, k)*binomial(n+k+1, k+1))/(n+1)}
%K nonn
%O 0,2
%A _Paul D. Hanna_, Feb 15 2011
|
