結果

問題 No.8030 ミラー・ラビン素数判定法のテスト
ユーザー pogyomo
提出日時 2024-04-03 05:06:13
言語 C++11（廃止可能性あり）
(gcc 13.3.0)
結果
AC  
実行時間 269 ms / 9,973 ms
コード長 1,230 bytes
コンパイル時間 726 ms
コンパイル使用メモリ 61,392 KB
実行使用メモリ 6,820 KB
最終ジャッジ日時 2024-09-30 23:37:41
合計ジャッジ時間 2,065 ms
ジャッジサーバーID
（参考情報） 		judge1 / judge2
このコードへのチャレンジ
（要ログイン）
ファイルパターン 結果
other AC * 10
権限があれば一括ダウンロードができます

ソースコード

diff # 

#include <iostream>
#include <vector>
using namespace std;

template<class T> T pow_mod(T A, T N, T M) {
    T res = 1 % M;
    A %= M;
    while (N) {
        if (N & 1) res = (res * A) % M;
        A = (A * A) % M;
        N >>= 1;
    }
    return res;
}

bool MillerRabin(long long N, vector<long long> A) {
    long long s = 0, d = N - 1;
    while (d % 2 == 0) {
        ++s;
        d >>= 1;
    }
    for (auto a : A) {
        if (N <= a) return true;
        long long t, x = pow_mod<__int128_t>(a, d, N);
        if (x != 1) {
            for (t = 0; t < s; ++t) {
                if (x == N - 1) break;
                x = __int128_t(x) * x % N;
            }
            if (t == s) return false;
        }
    }
    return true;
}

bool is_prime(long long N) {
    if (N <= 1) return false;
    if (N == 2) return true;
    if (N % 2 == 0) return false;
    if (N < 4759123141LL)
        return MillerRabin(N, {2, 7, 61});
    else
        return MillerRabin(N, {2, 325, 9375, 28178, 450775, 9780504, 1795265022});
}

int main() {
	int n;
	cin >> n;
	for (int i = 0; i < n; i++) {
		long long x;
		cin >> x;
		if (is_prime(x)) {
			cout << x << " " << 1 << endl;
		} else {
			cout << x << " " << 0 << endl;
		}
	}
}
0