login
Binomial transform of A084624.
2

%I #6 Mar 24 2023 17:18:43

%S 1,3,8,21,55,143,366,919,2265,5491,13125,31000,72485,168042,386709,

%T 884161,2009742,4543830,10222264,22891099,51041560,113359224,

%U 250839510,553173006,1216070081,2665518207,5826533103,12703217438,27628250142

%N Binomial transform of A084624.

%H G. C. Greubel, <a href="/A084625/b084625.txt">Table of n, a(n) for n = 0..1000</a>

%F a(n) = Sum_{k=0..n} C(n, k)*floor(C(k+5, 5)/C(k+2, 2)).

%F a(n) = 2^n + Sum_{k=0..n} binomial(n,k)*floor(k*(k^2 +12*k +47)/60). - _G. C. Greubel_, Mar 24 2023

%t a[n_]:= a[n]= 2^n +Sum[Binomial[n,j]*Floor[j*(j^2+12*j+47)/60], {j,0, n}];

%t Table[a[n], {n,0,50}] (* _G. C. Greubel_, Mar 24 2023 *)

%o (Magma)

%o A084625:= func< n | (&+[Binomial(n,j)*Floor(Binomial(j+5,3)/10): j in [0..n]]) >;

%o [A084625(n): n in [0..50]]; // _G. C. Greubel_, Mar 24 2023

%o (SageMath)

%o def A084625(n): return sum(binomial(n,j)*(binomial(j+5,3)//10) for j in range(n+1))

%o [A084625(n) for n in range(51)] # _G. C. Greubel_, Mar 24 2023

%Y Cf. A084624, A084626, A084627, A084628, A084630, A084631.

%K easy,nonn

%O 0,2

%A _Paul Barry_, Jun 01 2003