Nice primes: Difference between revisions

Content added Content deleted
Line 234: Line 234:
<lang julia>using Primes
<lang julia>using Primes


isnice(n, base=10) = isprime(n) && mod1(n, base - 1) in [2, 3, 5, 7, 11, 13, 17, 19])
function isnice(n, base=10)
pdroot(i) = (x = sum(digits(i, base=base)); x > 9 ? pdroot(x) : x in [2, 3, 5, 7])
return isprime(n) && pdroot(n)
end


filter_open_interval(500, 1000, isnice)
filter_open_interval(500, 1000, isnice)