結果

問題 No.1235 ζ関数
ユーザー こまるこまる
提出日時 2020-10-02 13:12:42
言語 Haskell
(9.10.1)
結果
AC  
実行時間 2 ms / 2,000 ms
コード長 1,256 bytes
コンパイル時間 9,173 ms
コンパイル使用メモリ 201,216 KB
実行使用メモリ 5,376 KB
最終ジャッジ日時 2024-07-08 07:07:15
合計ジャッジ時間 9,966 ms
ジャッジサーバーID
(参考情報)
judge3 / judge4
このコードへのチャレンジ
(要ログイン)

テストケース

テストケース表示
入力 結果 実行時間
実行使用メモリ
testcase_00 AC 2 ms
5,248 KB
testcase_01 AC 1 ms
5,376 KB
testcase_02 AC 2 ms
5,376 KB
testcase_03 AC 1 ms
5,376 KB
testcase_04 AC 2 ms
5,376 KB
testcase_05 AC 2 ms
5,376 KB
testcase_06 AC 2 ms
5,376 KB
testcase_07 AC 1 ms
5,376 KB
testcase_08 AC 2 ms
5,376 KB
testcase_09 AC 2 ms
5,376 KB
testcase_10 AC 1 ms
5,376 KB
testcase_11 AC 1 ms
5,376 KB
testcase_12 AC 1 ms
5,376 KB
testcase_13 AC 2 ms
5,376 KB
testcase_14 AC 2 ms
5,376 KB
testcase_15 AC 1 ms
5,376 KB
testcase_16 AC 2 ms
5,376 KB
testcase_17 AC 1 ms
5,376 KB
testcase_18 AC 1 ms
5,376 KB
testcase_19 AC 1 ms
5,376 KB
testcase_20 AC 2 ms
5,376 KB
権限があれば一括ダウンロードができます
コンパイルメッセージ
Loaded package environment from /home/judge/.ghc/x86_64-linux-9.8.2/environments/default
[1 of 2] Compiling Main             ( Main.hs, Main.o )
[2 of 2] Linking a.out

ソースコード

diff #
プレゼンテーションモードにする

import Control.Monad
import Control.Monad.ST
import Data.Bits
import qualified Data.Vector as V
import qualified Data.Vector.Mutable as VM
import qualified Data.Vector.Unboxed as VU
import qualified Data.Vector.Unboxed.Mutable as VUM
main :: IO ()
main = readLn >>= zeta >>= print
zeta :: Int -> IO Double
zeta n = do
let primes = getPrimeVector 10000
xm <- V.unsafeThaw $ V.fromList [1.0]
forM_ [0 .. (VU.length primes - 1)] $ \i -> do
let p = fromIntegral $ primes VU.! i
x = 1.0 / (1.0 - p ** (- (fromIntegral n)))
VM.modify xm (\a -> a * x) 0
(V.! 0) <$> V.unsafeFreeze xm
getPrimeVector :: Int -> VU.Vector Int
getPrimeVector top = VU.filter (/= -1) . VU.imap (\i check -> if i == 0 then 2 else if check then i * 2 + 1 else -1) $! runST $ do
let m = (top - 1) `shiftR` 1
r = floor . sqrt . fromIntegral $ (top + 1)
sieve <- VU.unsafeThaw $ VU.replicate (m + 1) True
forM_ [1 .. r `unsafeShiftR` 1] $ \i -> do
isPrime <- VUM.unsafeRead sieve i
when isPrime $ do
forM_ [2 * i * (i + 1), 2 * i * (i + 2) + 1 .. m] $ \j -> do
VUM.unsafeWrite sieve j False
VU.unsafeFreeze sieve
הההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההה
XXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXX
0