結果
問題 | 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
ソースコード
importControl.MonadimportControl.Monad.STimportData.BitsimportqualifiedData.VectorasVimportqualifiedData.Vector.MutableasVMimportqualifiedData.Vector.UnboxedasVUimportqualifiedData.Vector.Unboxed.MutableasVUMmain::IO()main = readLn >>= zeta >>= printzeta::Int->IODoublezeta n = dolet primes = getPrimeVector 10000xm <- V.unsafeThaw $ V.fromList [1.0]forM_ [0 .. (VU.length primes - 1)] $ \i -> dolet p = fromIntegral $ primes VU.! ix = 1.0 / (1.0 - p ** (- (fromIntegral n)))VM.modify xm (\a -> a * x) 0(V.! 0) <$> V.unsafeFreeze xmgetPrimeVector::Int->VUVectorIntgetPrimeVector top = VU.filter (/= -1) . VU.imap (\i check -> if i == 0 then 2 else if check then i * 2 + 1 else -1) $! runST $ dolet m = (top - 1) `shiftR` 1r = floor . sqrt . fromIntegral $ (top + 1)sieve <- VU.unsafeThaw $ VU.replicate (m + 1) TrueforM_ [1 .. r `unsafeShiftR` 1] $ \i -> doisPrime <- VUM.unsafeRead sieve iwhen isPrime $ doforM_ [2 * i * (i + 1), 2 * i * (i + 2) + 1 .. m] $ \j -> doVUM.unsafeWrite sieve j FalseVU.unsafeFreeze sieve