cp's OEIS Frontend

This is a front-end for the Online Encyclopedia of Integer Sequences, made by Christian Perfect. The idea is to provide OEIS entries in non-ancient HTML, and then to think about how they're presented visually. The source code is on GitHub.

A047088 a(n) = A047080(2*n+1, n+2).

This page as a plain text file.
%I A047088 #16 Oct 31 2022 07:34:09
%S A047088 1,4,12,37,118,380,1229,3989,12987,42394,138709,454768,1493690,
%T A047088 4913969,16189534,53407853,176397299,583242159,1930349545,6394665589,
%U A047088 21201345460,70346920007,233581374587,776105485336,2580316142887,8583746045611,28570407158100
%N A047088 a(n) = A047080(2*n+1, n+2).
%H A047088 G. C. Greubel, <a href="/A047088/b047088.txt">Table of n, a(n) for n = 1..1000</a>
%F A047088 a(n+4) = ((16*n^5 + 324*n^4 + 2624*n^3 + 10509*n^2 + 20655*n + 15930)*a(n+3) - (8*n^5 + 148*n^4 + 1090*n^3 + 3953*n^2 + 7365*n + 5994)*a(n+2) + (4*n^4 + 84*n^3 + 701*n^2 + 2451*n + 2646)*a(n+1) - (n-3)*(n+6)*(2*n+7)*(2*n^2 + 23*n + 72)*a(n) )/((n+3)*(n+6)*(2*n+5)*(2*n^2 + 19*n + 51)). - _G. C. Greubel_, Oct 31 2022
%t A047088 A[n_, k_]:= Sum[(-1)^j*(n+k-3*j)!/(j!*(n-2*j)!*(k-2*j)!), {j,0,Floor[(n+k)/3]}] - Sum[(-1)^j*(n+k-3*j-2)!/(j!*(n-2*j-1)!*(k-2*j-1)!), {j,0,Floor[(n+k-2)/3]}];
%t A047088 Table[A[n-1, n+2], {n, 50}] (* _G. C. Greubel_, Oct 31 2022 *)
%o A047088 (Magma)
%o A047088 F:=Factorial;
%o A047088 p:= func< n,k | (&+[ (-1)^j*F(n+k-3*j)/(F(j)*F(n-2*j)*F(k-2*j)): j in [0..Min(Floor(n/2), Floor(k/2))]]) >;
%o A047088 q:= func< n,k | n eq 0 or k eq 0 select 0 else (&+[ (-1)^j*F(n+k-3*j-2)/(F(j)*F(n-2*j-1)*F(k-2*j-1)) : j in [0..Min(Floor((n-1)/2), Floor((k-1)/2))]]) >;
%o A047088 A:= func< n,k | p(n,k) - q(n,k) >;
%o A047088 [A(n-1,n+2): n in [1..50]]; // _G. C. Greubel_, Oct 31 2022
%o A047088 (SageMath)
%o A047088 f=factorial
%o A047088 def p(n,k): return sum( (-1)^j*f(n+k-3*j)/(f(j)*f(n-2*j)*f(k-2*j)) for j in range(1+min((n//2), (k//2))) )
%o A047088 def q(n,k): return sum( (-1)^j*f(n+k-3*j-2)/(f(j)*f(n-2*j-1)*f(k-2*j-1)) for j in range(1+min(((n-1)//2), ((k-1)//2))) )
%o A047088 def A(n,k): return p(n,k) - q(n,k)
%o A047088 [A(n-1,n+2) for n in range(1,50)] # _G. C. Greubel_, Oct 31 2022
%Y A047088 Cf. A047080, A047081, A047082, A047083, A047084, A047085, A047086, A047087.
%K A047088 nonn
%O A047088 1,2
%A A047088 _Clark Kimberling_
%E A047088 Corrected and extended by _Sean A. Irvine_, May 11 2021