fork download
copy 
  1. {-# OPTIONS_GHC -O2 #-}
  2. module Main where
  3.  
  4. import Data.Array
  5. import Data.List (unfoldr)
  6. import Control.Monad
  7.  
  8. primesFromTo a b = (if a<3 then [2] else []) 
  9.                      ++ [i | i <- [o,o+2..b], ar ! i]
  10.    where 
  11.     o  = max (if even a then a+1 else a) 3
  12.     r  = floor.sqrt.fromInteger $ b+1
  13.     ar = accumArray (\a b-> False) True (o,b) 
  14.           [(i,()) | p <- [3,5..r]
  15.                     , let q  = p*p 
  16.                           s  = 2*p 
  17.                           (n,x) = quotRem (o - q) s 
  18.                           q' = if  o <= q  then q
  19.                                else  q + (n + signum x)*s
  20.                     , i <- [q',q'+s..b] ]
  21. main = do
  22.   a <- getLine
  23.   b <- getLine
  24.   let ps = primesFromTo (read a) (read b)
  25.   forM_ (groups 4 ps) print
  26.   forM_ (groups 20 $ zipWith (-) (tail ps) ps) print
  27.  
  28. groups = unfoldr . listMaybe . splitAt
  29.  
  30. listMaybe _ [] = Nothing 
  31. listMaybe f xs = Just (f xs)
stdin
copy 
500111222000000 
500111222002350
compilation info 
[1 of 1] Compiling Main             ( prog.hs, prog.o )
Linking prog ...
stdout
copy 
[500111222000017,500111222000053,500111222000057,500111222000063]
[500111222000069,500111222000123,500111222000167,500111222000179]
[500111222000197,500111222000233,500111222000333,500111222000387]
[500111222000389,500111222000479,500111222000503,500111222000603]
[500111222000659,500111222000707,500111222000713,500111222000849]
[500111222000879,500111222000893,500111222000939,500111222001043]
[500111222001049,500111222001071,500111222001089,500111222001187]
[500111222001203,500111222001221,500111222001227,500111222001263]
[500111222001311,500111222001337,500111222001397,500111222001407]
[500111222001497,500111222001499,500111222001529,500111222001553]
[500111222001593,500111222001637,500111222001733,500111222001739]
[500111222001761,500111222001773,500111222001779,500111222001791]
[500111222001827,500111222001847,500111222001857,500111222001871]
[500111222001889,500111222001901,500111222001983,500111222002009]
[500111222002019,500111222002037,500111222002039,500111222002061]
[500111222002097,500111222002103,500111222002141,500111222002159]
[500111222002169,500111222002199,500111222002211,500111222002279]
[500111222002309,500111222002333]
[36,4,6,6,54,44,12,18,36,100,54,2,90,24,100,56,48,6,136,30]
[14,46,104,6,22,18,98,16,18,6,36,48,26,60,10,90,2,30,24,40]
[44,96,6,22,12,6,12,36,20,10,14,18,12,82,26,10,18,2,22,36]
[6,38,18,10,30,12,68,30,24]
https://ideone.com/HO2Hy
language:
Haskell (ghc 8.4.4)
created:
13 years ago
visibility:
public

Share or Embed source code

Discover > Sphere Engine API

The brand new service which powers Ideone!

Discover > IDE Widget

Widget for compiling and running the source code in a web browser!