Primality by Wilson's theorem: Difference between revisions

m
Line 195:
end
 
wilsonprimesbetween(n, m) = [i for i in n:m if iswilsonprime(i)]
function nwilsonprimesabove(N, bottom=1)
wprimes = Int[]
for i in bottom:typemax(Int)
if iswilsonprime(i)
push!(wprimes, i)
if length(wprimes) >= N
break
end
end
end
return wprimes
end
 
println("First 120 Wilson primes: ", nwilsonprimesabovewilsonprimesbetween(1201, 1000)[1:120])
println("\nThe first 40 Wilson primes above 7900 are: ", nwilsonprimesabovewilsonprimesbetween(407900, 79009000)[1:40])
 
println("\nThe first 40 Wilson primes above 7900 are: ", nwilsonprimesabove(40, 7900))
</lang>{{out}}
<pre>
4,108

edits