結果

問題 No.3030 ミラー・ラビン素数判定法のテスト
ユーザー glretoglreto
提出日時 2021-03-20 16:49:25
言語 C++14
(gcc 12.3.0 + boost 1.83.0)
結果
AC  
実行時間 109 ms / 9,973 ms
コード長 1,646 bytes
コンパイル時間 1,845 ms
コンパイル使用メモリ 167,188 KB
実行使用メモリ 6,944 KB
最終ジャッジ日時 2024-04-28 09:42:50
合計ジャッジ時間 2,612 ms
ジャッジサーバーID
（参考情報） 		judge5 / judge2
このコードへのチャレンジ
（要ログイン）

テストケース

テストケース表示
入力 結果 実行時間
実行使用メモリ
testcase_00 AC 2 ms
6,816 KB
testcase_01 AC 2 ms
6,940 KB
testcase_02 AC 2 ms
6,940 KB
testcase_03 AC 2 ms
6,944 KB
testcase_04 AC 70 ms
6,940 KB
testcase_05 AC 70 ms
6,940 KB
testcase_06 AC 42 ms
6,940 KB
testcase_07 AC 42 ms
6,940 KB
testcase_08 AC 42 ms
6,944 KB
testcase_09 AC 109 ms
6,944 KB
権限があれば一括ダウンロードができます

ソースコード

diff # 

#include <bits//stdc++.h>
using namespace std;
#define rep(i, n) for (int i = 0; i < (int)(n); i++)
#define req(i,n) for(int i = 1;i <=  n; i++)
#define rrep(i,n) for(ll i = n-1;i >= 0;i--)
#define ALL(obj) begin(obj), end(obj)
#define RALL(a) rbegin(a),rend(a)
typedef long long ll;
const ll mod = 1e9 + 7;
unsigned long long mod_mul(unsigned long long a, unsigned long long b, unsigned long long M) {
    long long ret = a * b - M * (unsigned long long)((long double)(a) * (long double)(b) / (long double)(M));
    return ret + M * (ret < 0) - M * (ret >= (ll)M);
}const unsigned long long numset[] = { 2, 325, 9375, 28178, 450775, 9780504, 1795265022ULL }; 
unsigned long long mod_pow(unsigned long long x, unsigned long long k, unsigned long long mod) {
    unsigned long long res = 1;
    while (k) {
        if (k & 1) res = mod_mul(res, x, mod);
        x = mod_mul(x, x, mod);
        k >>= 1;
    }return res;
}
// check if n is prime
bool check(unsigned long long n) {
    if (n < 2 || ((n % 6 != 1) && (n % 6 != 5))) return (n == 2) || (n == 3);
    unsigned long long d = n - 1, s = 0;
    while (d % 2 == 0) {
        d /= 2; s++;
    }for (unsigned long long a : numset) {
        if (a % n == 0) return 1;
        unsigned long long res = mod_pow(a, d, n);
        if (res == 1) continue;
        bool ok = 1;
        for (unsigned int r = 0; r < s; r++) {
            if (res == n - 1) {
                ok = 0;break;
            }res = mod_mul(res, res, n);
        }if (ok) return 0;
    }return 1;
}
int main() {
    int n; cin >> n;
    rep(i, n) {
        ll x; cin >> x;
        cout << x << " " << check(x) << endl;
    }
}
0