結果
問題 | No.8030 ミラー・ラビン素数判定法のテスト |
ユーザー | amylase_pepsin |
提出日時 | 2020-10-15 10:01:23 |
言語 | C++14 (gcc 12.3.0 + boost 1.83.0) |
結果 |
WA
|
実行時間 | - |
コード長 | 2,233 bytes |
コンパイル時間 | 1,568 ms |
コンパイル使用メモリ | 171,160 KB |
実行使用メモリ | 6,820 KB |
最終ジャッジ日時 | 2024-11-18 18:30:16 |
合計ジャッジ時間 | 2,419 ms |
ジャッジサーバーID (参考情報) |
judge3 / judge4 |
(要ログイン)
テストケース
テストケース表示入力 | 結果 | 実行時間 実行使用メモリ |
---|---|---|
testcase_00 | AC | 2 ms
6,820 KB |
testcase_01 | AC | 2 ms
6,816 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 "bits/stdc++.h" #include <vector> namespace amylase { using li = long long; li _powmod(const li x, const li n, const li mod) { if (n == 0) { return 1; } li sq = _powmod(x, n / 2, mod); if (n & 1) { return sq * sq % mod * x % mod; } else { return sq * sq % mod; } } bool _miller_rabin(const long long x, const std::vector<long long>& witnesses) { if (x <= 1) { return false; } if ((x & 1) == 0) { return x == 2; } long long d = x - 1; long long s = 0; while ((d & 1) == 0) { d >>= 1; s++; } for (const auto& a : witnesses) { if (a % x <= 1) { continue; } bool is_composite = true; is_composite &= _powmod(a, d, x) != 1; long long dd = d; for (int i = 0; i < s; ++i) { is_composite &= _powmod(a, dd, x) != x - 1; dd <<= 1; } if (is_composite) { return false; } } return true; } bool is_prime(const long long x) { return _miller_rabin(x, {2LL, 325LL, 9375LL, 28178LL, 450775LL, 9780504LL, 1795265022LL}); } std::vector<long long> factor(long long x) { std::vector<long long> factors; long long p = 2; while (p * p <= x) { while (x % p == 0) { factors.emplace_back(p); x /= p; } p++; } if (x > 1) { factors.emplace_back(x); } return factors; } /** * solves ax + by = gcd(a, b) * @param a * @param b * @param x ref to variable to store root x * @param y ref to variable to store root x * @return gcd(a, b) */ long long extgcd(long long a, long long b, long long& x, long long& y) { long long d = a; if (b != 0) { d = extgcd(b, a % b, y, x); y -= (a / b) * x; } else { x = 1; y = 0; } return d; } } // namespace amylase using namespace std; typedef long long li; int main() { cin.tie(0); ios::sync_with_stdio(false); li n; cin >> n; for (int i = 0; i < n; ++i) { li x; cin >> x; cout << x << " "; li ispr = amylase::is_prime(x); cout << ispr << '\n'; } return 0; }