結果
問題 | No.3030 ミラー・ラビン素数判定法のテスト |
ユーザー | hamray |
提出日時 | 2020-01-05 16:03:16 |
言語 | C++11 (gcc 11.4.0) |
結果 |
AC
|
実行時間 | 501 ms / 9,973 ms |
コード長 | 1,747 bytes |
コンパイル時間 | 562 ms |
コンパイル使用メモリ | 62,456 KB |
実行使用メモリ | 5,376 KB |
最終ジャッジ日時 | 2024-04-28 09:27:24 |
合計ジャッジ時間 | 2,365 ms |
ジャッジサーバーID (参考情報) |
judge1 / judge3 |
(要ログイン)
テストケース
テストケース表示入力 | 結果 | 実行時間 実行使用メモリ |
---|---|---|
testcase_00 | AC | 2 ms
5,248 KB |
testcase_01 | AC | 2 ms
5,376 KB |
testcase_02 | AC | 1 ms
5,376 KB |
testcase_03 | AC | 2 ms
5,376 KB |
testcase_04 | AC | 282 ms
5,376 KB |
testcase_05 | AC | 249 ms
5,376 KB |
testcase_06 | AC | 79 ms
5,376 KB |
testcase_07 | AC | 83 ms
5,376 KB |
testcase_08 | AC | 82 ms
5,376 KB |
testcase_09 | AC | 501 ms
5,376 KB |
コンパイルメッセージ
main.cpp: In function ‘int main()’: main.cpp:102:10: warning: ignoring return value of ‘int scanf(const char*, ...)’ declared with attribute ‘warn_unused_result’ [-Wunused-result] 102 | scanf("%d", &n); | ~~~~~^~~~~~~~~~ main.cpp:107:14: warning: ignoring return value of ‘int scanf(const char*, ...)’ declared with attribute ‘warn_unused_result’ [-Wunused-result] 107 | scanf("%lld", &Q); | ~~~~~^~~~~~~~~~~~
ソースコード
#include <iostream> #include <cstdio> #include <stdlib.h> #include <cstdint> #include <vector> #include <fstream> #include <iostream> using std::vector; using Int32 = __int32_t; using Int64 = __int64_t; using Int128 = __int128_t; using ll = long long int; Int128 powmod(Int128 x, Int128 n, ll mod) { Int128 res = 1; while (n > 0) { if (n & 1) { res = (res * x) % mod; } x = x * x % mod; n >>= 1; } return res; } Int128 mulmod(Int128 a, Int128 b, ll mod) { Int128 x = 0, y = a % mod; while (b > 0) { if (b % 2 == 1) { x = (x + y) % mod; } y = (y * 2) % mod; b /= 2; } return x % mod; } ll gcd(ll a, ll b) { if (a < b) std::swap(a, b); while (b > 0) { ll tmp = a; a = b; b = tmp % b; } return a; } inline ll funct(ll x, ll c) { return x * x + c; } bool MillerRabinTest(Int128 n) { constexpr Int128 witness[] = {2, 3, 5, 7, 11, 13, 17, 19, 23, 29, 31, 37, -1}; if (n == 2) return true; if (n <= 1 || n % 2 == 0) return false; Int128 d = n - 1; int s = 0; while (!(d & 1)) { ++s; d >>= 1; } for (int i = 0; witness[i] != -1 && witness[i] < n; i++) { Int128 y = powmod(witness[i], d, n); if (y == 1) continue; for (int j = 0; j < s - 1 && y != n - 1; j++) y = mulmod(y, y, n); if (y != n - 1) return false; } return true; } int main() { int n; scanf("%d", &n); ll Q; while (n--) { scanf("%lld", &Q); printf("%lld %d\n", Q, MillerRabinTest(Q)); } }