%I #16 Jan 08 2022 17:14:45
%S 1,0,0,0,1,6,29,177,1269,10366,95020,966032,10787363,131246245,
%T 1728075563,24481070470,371296235455,6002622473192,103045019123133,
%U 1871984514070245,35879703186346359,723574014257984910
%N Nearest integer to Gamma(n + 1/6)/Gamma(1/6).
%C Gamma(n + 1/6)/Gamma(1/6) = 1, 1/6, 7/36, 91/216, 1729/1296, 43225/7776, 1339975/46656, 49579075/279936, 2131900225/1679616, ... - _R. J. Mathar_, Sep 04 2016
%p Digits := 64:f := proc(n,x) round(GAMMA(n+x)/GAMMA(x)); end;
%t With[{c=Gamma[1/6]},Table[Round[Gamma[n+1/6]/c],{n,0,30}]] (* _Harvey P. Dale_, Jan 08 2022 *)
%o (PARI) a(n) = round(gamma(n+1/6)/gamma(1/6)); \\ _Michel Marcus_, Jun 14 2018
%Y Cf. A008542, A020081, A020126.
%K nonn
%O 0,6
%A _Simon Plouffe_
|