A020085 a(n) = floor( Gamma(n + 1/5)/Gamma(1/5) ).

1, 0, 0, 0, 1, 7, 36, 228, 1647, 13507, 124268, 1267542, 14196477, 173197024, 2286200718, 32464050199, 493453563029, 7993947721071, 137495900802423, 2502425394604113, 48046567576398986, 970540665043259525

Offset: 0

Author

Simon Plouffe

Keywords

• nonn

Links

• G. C. Greubel, Table of n, a(n) for n = 0..450

Crossrefs

Cf. A020040, A020130.

Cf. A020082, A020083, A020084.

Programs

Magma

[Truncate(Gamma(n + 1/5)/Gamma(1/5)): n in [0..30]]; // G. C. Greubel, Nov 19 2018

Maple

Digits:= 64: x:=1/5: f:= proc(n,x) trunc(GAMMA(n+x)/GAMMA(x)); end;
seq(floor(pochhammer(1/5,n)), n = 0..25); # G. C. Greubel, Nov 17 2019

Mathematica

Floor[Pochhammer[1/5, Range[0, 25]]] (* G. C. Greubel, Nov 19 2018 *)

PARI

vector(26, n, my(x=1/5); gamma(n-1+x)\gamma(x) ) \\ G. C. Greubel, Nov 19 2018

Sage

[int(gamma(n + 1/5)/gamma(1/5))  for n in range(30)] # G. C. Greubel, Nov 19 2018