A020089 Integer part of Gamma(n + 1/3)/Gamma(1/3).

1, 0, 0, 1, 3, 14, 79, 505, 3710, 30920, 288591, 2982109, 33797240, 416832633, 5557768440, 79661347641, 1221473997164, 19950741953681, 345812860530487, 6339902443058941, 122571447232472865, 2492286093726948268

Offset

0,5

Links

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

Maple

Digits:= 64: x:=1/3: f:= proc(n, x) trunc(GAMMA(n+x)/GAMMA(x)); end;
# alternative
f:= n -> trunc(simplify(GAMMA(n+1/3)/GAMMA(1/3))): map(f, [$0..30]); # Robert Israel, Nov 20 2018

Mathematica

Table[IntegerPart[Pochhammer[1/3, n]], {n, 0, 30}] (* G. C. Greubel, Nov 19 2018 *)

Prog

(PARI) default(realprecision, 100); for(n=0, 30, print1(truncate( gamma(n+1/3)/gamma(1/3)), ", ")) \\ G. C. Greubel, Nov 19 2018
(Magma) [Truncate(Gamma(n + 1/3)/Gamma(1/3)): n in [0..30]]; // G. C. Greubel, Nov 19 2018
(Sage) [int(gamma(n + 1/3)/gamma(1/3)) for n in range(30)] # G. C. Greubel, Nov 19 2018

Crossrefs

Cf. A020044, A020134.

Sequence in context: A059276 A003169 A086621 * A218677 A353079 A305128

Adjacent sequences: A020086 A020087 A020088 * A020090 A020091 A020092

Keyword

nonn

Author

Simon Plouffe

Status

approved