OFFSET
1,1
LINKS
David A. Corneth, Table of n, a(n) for n = 1..10000
PROG
(PARI) is(n) = omega(n) == 8
(PARI) A246655(lim)=my(v=List(primes([2, lim\=1]))); for(e=2, logint(lim, 2), forprime(p=2, sqrtnint(lim, e), listput(v, p^e))); Set(v)
list(lim, pr=8)=if(pr==1, return(A246655(lim))); my(v=List(), pr1=pr-1, mx=prod(i=1, pr1, prime(i))); forprime(p=prime(pr), lim\mx, my(u=list(lim\p, pr1)); for(i=1, #u, listput(v, p*u[i]))); Set(v) \\ Charles R Greathouse IV, Feb 03 2023
CROSSREFS
KEYWORD
nonn
AUTHOR
David A. Corneth, Oct 10 2021
STATUS
approved