結果
問題 |
No.8030 ミラー・ラビン素数判定法のテスト
|
ユーザー |
![]() |
提出日時 | 2019-10-16 08:31:31 |
言語 | C++14 (gcc 13.3.0 + boost 1.87.0) |
結果 |
WA
|
実行時間 | - |
コード長 | 931 bytes |
コンパイル時間 | 1,610 ms |
コンパイル使用メモリ | 169,388 KB |
実行使用メモリ | 10,148 KB |
最終ジャッジ日時 | 2024-11-18 17:22:14 |
合計ジャッジ時間 | 67,994 ms |
ジャッジサーバーID (参考情報) |
judge5 / judge3 |
(要ログイン)
ファイルパターン | 結果 |
---|---|
other | AC * 1 WA * 3 TLE * 6 |
ソースコード
#include "bits/stdc++.h" using namespace std; #define int long long random_device rd; mt19937 mt(rd()); int modpro(int x, int n, int md) { int res = 0; for (; n; n >>= 1) { if (n & 1)res =(res+x)%md; x = (x + x) % md; } return res; } int modpow(int x, int n, int md) { int res = 1; for (; n; n >>= 1) { if (n & 1)res = modpro(res, x, md); x = modpro(x, x, md); } return res; } int gcd(int a, int b) { if (!b)return a; return gcd(b, a%b); } bool isprime(int p) { if (p == 2 || p == 3 || p == 5)return true; if (p == 1 || p == 4)return false; for(int z=0;z<1;++z) { int b = mt() % (p - 2) + 2; if (gcd(b, p) > 1) { return false; } if (modpow(b, p - 1, p) > 1) { return false; } } return true; } signed main() { ios::sync_with_stdio(false); cin.tie(nullptr); int n; cin >> n; int p; for (int i = 0; i < n; i++) { cin >> p; cout << p; if (isprime(p))puts(" 1"); else puts(" 0"); } }