A020082 a(n) = floor( Gamma(n + 4/5)/Gamma(4/5) ).

1, 0, 1, 4, 15, 73, 426, 2900, 22624, 199094, 1951127, 21072175, 248651668, 3182741351, 43921830647, 650043093588, 10270680878693, 172547438762055, 3071344409964594, 57741274907334370, 1143277243165220533

Offset: 0

Author

Simon Plouffe

Keywords

• nonn

Links

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

Crossrefs

Cf. A020037, A020127.

Cf. A020083, A020084, A020085.

Programs

Magma

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

Maple

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

Mathematica

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

PARI

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

Sage

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