Mersenne primes: Difference between revisions

m
Line 563:
 
lucasLehmer :: Int -> Bool
lucasLehmer p = go 4 [1..p - 2] == 0
where
go = foldr (\a_ b -> (b ^ 2 - 2) `mod` m) 4
m = 2 ^ p - 1
 
main = mapM_ (printf "M %d\n") $ take 20 mersenne
Anonymous user