OFFSET
1,2
PROG
(PARI) a(n) = b=prod(i=1, n, prime(i)); k={my(k=1); while(omega(b-k)<(n-1), k++); k; }; b-k
(PARI) a(n) = if(n==1, return(1)); my(A=vecprod(primes(n-1)), B=vecprod(primes(n))-1); (f(m, p, j) = my(r=0); forprime(q=p, sqrtnint(B\m, j), my(v=m*q); while(v <= B, if(j==1, if(v>=A && v > r, r = v), if(v*(q+1) <= B, r = max(r, f(v, q+1, j-1)))); v *= q)); r); f(1, 2, n-1); \\ Daniel Suteu, Sep 22 2023
CROSSREFS
KEYWORD
nonn
AUTHOR
Gionata Neri, Sep 29 2017
EXTENSIONS
a(12)-a(19) from Daniel Suteu, Sep 22 2023
STATUS
approved