{-# LANGUAGE BangPatterns #-} import Control.Monad import Control.Monad.Fix import Data.IORef import qualified Data.Vector.Fusion.Stream.Monadic as VFSM main :: IO () main = do l <- readLn :: IO Int ans <- newIORef (0 :: Int) rep l $ \m -> repStep (m `mod` 2 + 1) m 2 $ \n -> do when (2 * m * (m + m) <= l) $ do when (gcd m n == 1) $ modifyIORef' ans succ print =<< readIORef ans stream :: Monad m => Int -> Int -> VFSM.Stream m Int stream !l !r = VFSM.Stream step l where step x | x * x < r = return $ VFSM.Yield x (x + 1) | otherwise = return $ VFSM.Done {-# INLINE [0] step #-} {-# INLINE [1] stream #-} rep :: Monad m => Int -> (Int -> m ()) -> m () rep n = flip VFSM.mapM_ (stream 2 n) {-# INLINE rep #-} streamStep :: Monad m => Int -> Int -> Int -> VFSM.Stream m Int streamStep !l !r !d = VFSM.Stream step l where step x | x < r = return $ VFSM.Yield x (x + d) | otherwise = return $ VFSM.Done {-# INLINE [0] step #-} {-# INLINE [1] streamStep #-} repStep :: Monad m => Int -> Int -> Int -> (Int -> m ()) -> m () repStep l r d = flip VFSM.mapM_ (streamStep l r d) {-# INLINE repStep #-}