結果
問題 | No.8030 ミラー・ラビン素数判定法のテスト |
ユーザー | とばり |
提出日時 | 2018-08-11 03:35:46 |
言語 | C++14 (gcc 13.3.0 + boost 1.87.0) |
結果 |
WA
|
実行時間 | - |
コード長 | 2,051 bytes |
コンパイル時間 | 808 ms |
コンパイル使用メモリ | 83,736 KB |
実行使用メモリ | 6,824 KB |
最終ジャッジ日時 | 2024-11-18 16:27:59 |
合計ジャッジ時間 | 1,607 ms |
ジャッジサーバーID (参考情報) |
judge4 / judge3 |
(要ログイン)
テストケース
テストケース表示入力 | 結果 | 実行時間 実行使用メモリ |
---|---|---|
testcase_00 | AC | 1 ms
6,820 KB |
testcase_01 | AC | 2 ms
6,820 KB |
testcase_02 | AC | 2 ms
6,816 KB |
testcase_03 | AC | 2 ms
6,816 KB |
testcase_04 | WA | - |
testcase_05 | WA | - |
testcase_06 | WA | - |
testcase_07 | WA | - |
testcase_08 | WA | - |
testcase_09 | WA | - |
ソースコード
#include <iostream>#include <random>using namespace std;using int64 = long long;namespace Random{static random_device rnd;static mt19937_64 mt(rnd());// n未満の非負整数をランダムに返すint64 randInt(int64 n){return mt() % n;}// [a, b)に属する整数をランダムに返すint64 randInt(int64 a, int64 b){int64 width = b - a;return randInt(width) + a;}};namespace Mod{const int64 MOD = 1e9 + 7;int64 powMod(int64 a, int64 n, int64 mod = MOD){int64 res = 1, p = a;while (n > 0){if (n & 1) res = (res * p) % mod;p = (p * p) % mod;n >>= 1;}return res;}};namespace Prime{enum Result { PROBABLY_PRIME, COMPOSITE };// 時間計算量:O(klog^3 n)// nが合成数なのに素数であるという判定をしてしまう確率はたかだか4^{-k}Result millerRabin(int64 n, int k){if (n <= 1) return COMPOSITE;if (n == 2) return PROBABLY_PRIME;if (n % 2 == 0) return COMPOSITE;int64 s = 0, d = n - 1;while (d % 2 == 0) { s++; d >>= 1; }for (int i = 0; i < k; i++){int64 a = Random::randInt(1, n);if (Mod::powMod(a, d, n) == 1) continue;bool all = true;for (int64 r = 0; r < s; r++){if (Mod::powMod(a, (1LL << r) * d, n) == n - 1){all = false;break;}}if (all) return COMPOSITE;}return PROBABLY_PRIME;}};int main(){int n;cin >> n;for (int i = 0; i < n; i++){int64 x;cin >> x;if (Prime::millerRabin(x, 30) == Prime::PROBABLY_PRIME){cout << x << " " << 1 << endl;}else{cout << x << " " << 0 << endl;}}return 0;}