login
a(n) = floor( Gamma(n+9/10)/Gamma(9/10) ).
4

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

%S 1,0,1,4,19,94,559,3857,30477,271252,2685395,29270806,348322596,

%T 4493361500,62457724853,930620100318,14796859595058,250066927156481,

%U 4476197996101023,84600142126309335,1683542828313555774

%N a(n) = floor( Gamma(n+9/10)/Gamma(9/10) ).

%H G. C. Greubel, <a href="/A020060/b020060.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(9/10,n)), n = 0..20); # _G. C. Greubel_, Nov 13 2019

%t Floor[Pochhammer[9/10, Range[0, 20]]] (* _G. C. Greubel_, Nov 13 2019 *)

%o (PARI) vector(21, n, my(x=9/10); gamma(n-1+x)\gamma(x) ) \\ _G. C. Greubel_, Nov 13 2019

%o (Magma) [Floor(Gamma(n+9/10)/Gamma(9/10)): n in [0..20]]; // _G. C. Greubel_, Nov 13 2019

%o (Sage) [floor(rising_factorial(9/10, n)) for n in (0..20)] # _G. C. Greubel_, Nov 13 2019

%Y Cf. A020015, A020105.

%Y Cf. A020061, A020062, A020063.

%K nonn

%O 0,4

%A _Simon Plouffe_