Square-free integers: Difference between revisions
Content added Content deleted
Line 888: | Line 888: | ||
const maxrootprime = Int64(floor(sqrt(1000000000145))) |
const maxrootprime = Int64(floor(sqrt(1000000000145))) |
||
sqprimes = map(x -> x * x, primes(2, maxrootprime)) |
const sqprimes = map(x -> x * x, primes(2, maxrootprime)) |
||
possdivisorsfor(n) = vcat(filter(x -> x <= n / 2, sqprimes), n in sqprimes ? n : []) |
possdivisorsfor(n) = vcat(filter(x -> x <= n / 2, sqprimes), n in sqprimes ? n : []) |
||
issquarefree(n) = all(x -> floor(n / x) != n / x, possdivisorsfor(n)) |
issquarefree(n) = all(x -> floor(n / x) != n / x, possdivisorsfor(n)) |