OFFSET
1,1
COMMENTS
EXAMPLE
PROG
(PARI) list(lim) = my(v=primes(primepi(lim)), t); forprime(p=2, sqrt(lim), t=p; while((t=t^2)<=lim, v=concat(v, t))); vecsort(v); \\ A050376
lista(nn) = my(vfd=list(nn), res=List()); for (n=1, nn, my(vless = select(x->(x<n), vfd)); if (#vless, my(vmore = select(x->(x>n), vfd)); if (#vmore, my(p = vecprod(vless)); if (numdiv(p*n) > numdiv(p*vmore[1]), listput(res, n)); ); ); ); res; \\ Michel Marcus, Jan 10 2022
CROSSREFS
KEYWORD
nonn
AUTHOR
J. Lowell, Jan 09 2022
EXTENSIONS
More terms from Jinyuan Wang, Jan 09 2022
STATUS
approved