login

Reminder: The OEIS is hiring a new managing editor, and the application deadline is January 26.

a(n) = floor( Gamma(n+1/7)/Gamma(1/7) ).
7

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

%S 1,0,0,0,1,4,23,143,1027,8370,76527,776207,8649165,105025575,

%T 1380336141,19521896861,295617295326,4772107767405,81807561726950,

%U 1484222905617529,28412267050392705,572304236300767347

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

%H G. C. Greubel, <a href="/A020079/b020079.txt">Table of n, a(n) for n = 0..445</a>

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

%p seq(floor(pochhammer(1/7,n)), n = 0..25); # _G. C. Greubel_, Nov 17 2019

%t Floor[Pochhammer[1/7, Range[0, 25]]] (* _G. C. Greubel_, Nov 17 2019 *)

%o (PARI) vector(26, n, my(x=1/7); gamma(n-1+x)\gamma(x) ) \\ _G. C. Greubel_, Nov 17 2019

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

%o (Sage) [floor(rising_factorial(1/7, n)) for n in (0..25)] # _G. C. Greubel_, Nov 17 2019

%Y Cf. A020074, A020075, A020076, A020077, A020078.

%K nonn

%O 0,6

%A _Simon Plouffe_