OFFSET
2,2
LINKS
G. C. Greubel, Table of n, a(n) for n = 2..1000
FORMULA
a(n) = p(n)/p(n-1), where p(n) = ( floor(n*log(n)) / Product_{j=2..pi(floor(n*log(n)))} prime(j) )!.
MATHEMATICA
p[n_]:=Factorial[Floor[n*Log[n]]]/ Product[Prime[i], {i, 2, PrimePi[Floor[n*Log[n]]]}];
Table[p[n]/p[n-1], {n, 2, 50}]
PROG
(Magma)
m:=50;
b:= [ #PrimesUpTo(n): n in [1..2+Floor(2*m*Log(2*m))] ];
f:= func< n | Factorial( Floor(n*Log(n)) )/(&*[ NthPrime(j): j in [2..b[Floor(n*Log(n))]] ]) >;
A088301:= func< n | n le 3 select n-1 else f(n)/f(n-1) >;
[A088301(n): n in [2..m]]; // G. C. Greubel, Dec 18 2022
(SageMath)
def p(n): return factorial( floor(n*log(n)) )/product(nth_prime(j) for j in (2..prime_pi(floor(n*log(n)))))
def A088301(n): return p(n)/p(n-1)
[A088301(n) for n in range(2, 50)] # G. C. Greubel, Dec 18 2022
CROSSREFS
KEYWORD
nonn,less
AUTHOR
Roger L. Bagula, Nov 04 2003
EXTENSIONS
Edited by G. C. Greubel, Dec 18 2022
STATUS
approved