A020049 - OEIS

%I #13 Sep 08 2022 08:44:44

%S 1,0,0,0,0,2,12,73,518,4192,38083,384010,4256112,51428023,672849973,

%T 9475970454,142929221023,2298778304796,39270796040272,710146895061597,

%U 13551969914092152,272168729108017392,5738224038694033364

%N a(n) = floor(Gamma(n+1/12)/Gamma(1/12)).

%H G. C. Greubel, <a href="/A020049/b020049.txt">Table of n, a(n) for n = 0..450</a>

%p Digits := 64:f := proc(n,x) trunc(GAMMA(n+x)/GAMMA(x)); end;

%p seq(floor(pochhammer(1/12,n)), n = 0..25); # _G. C. Greubel_, Dec 01 2019

%t Floor[Pochhammer[1/12, Range[0, 25]]] (* _G. C. Greubel_, Dec 01 2019 *)

%o (PARI) x=1/12; vector(26, n, gamma(n-1+x)\gamma(x) ) \\ _G. C. Greubel_, Dec 01 2019

%o (Magma) [Floor(Gamma(n+1/12)/Gamma(1/12)): n in [0..25]]; // _G. C. Greubel_, Dec 01 2019

%o (Sage) [floor(rising_factorial(1/12, n)) for n in (0..25)] # _G. C. Greubel_, Dec 01 2019

%Y Cf. A020004, A020094.

%K nonn

%O 0,6

%A _Simon Plouffe_