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

 

Logo
Hints
(Greetings from The On-Line Encyclopedia of Integer Sequences!)
A346680 a(n) = Sum_{k=0..n} (-1)^(n-k) * binomial(4*k,k) / (3*k + 1). 5

%I #15 Aug 05 2021 13:12:24

%S 1,0,4,18,122,847,6237,47583,373149,2989111,24354777,201214021,

%T 1681719343,14193619647,120800146953,1035593096367,8934344395053,

%U 77510878324671,675799844685937,5918354494345863,52037647837001257,459200394617540288,4065477723321641932

%N a(n) = Sum_{k=0..n} (-1)^(n-k) * binomial(4*k,k) / (3*k + 1).

%H Seiichi Manyama, <a href="/A346680/b346680.txt">Table of n, a(n) for n = 0..1000</a>

%F G.f. A(x) satisfies: A(x) = 1 / (1 + x) + x * (1 + x)^3 * A(x)^4.

%F a(n) ~ 2^(8*n + 17/2) / (283 * sqrt(Pi) * n^(3/2) * 3^(3*n + 3/2)). - _Vaclav Kotesovec_, Jul 30 2021

%F D-finite with recurrence 3*n*(3*n-1)*(3*n+1)*a(n) -(n-1)*(229*n^2-155*n+24)*a(n-1) -8*(4*n-3)*(2*n-1)*(4*n-1)*a(n-2)=0. - _R. J. Mathar_, Aug 05 2021

%t Table[Sum[(-1)^(n - k) Binomial[4 k, k]/(3 k + 1), {k, 0, n}], {n, 0, 22}]

%t nmax = 22; A[_] = 0; Do[A[x_] = 1/(1 + x) + x (1 + x)^3 A[x]^4 + O[x]^(nmax + 1) // Normal, nmax + 1]; CoefficientList[A[x], x]

%o (PARI) a(n) = sum(k=0, n, (-1)^(n-k)*binomial(4*k, k)/(3*k + 1)); \\ _Michel Marcus_, Jul 29 2021

%Y Cf. A002293, A032357, A188678, A345367, A346664.

%Y Cf. A346681, A346682, A346683, A346684.

%K nonn

%O 0,3

%A _Ilya Gutkovskiy_, Jul 29 2021

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.

License Agreements, Terms of Use, Privacy Policy. .

Last modified March 29 00:26 EDT 2024. Contains 371264 sequences. (Running on oeis4.)