結果

問題 No.8030 ミラー・ラビン素数判定法のテスト
ユーザー Dr_Geometry
提出日時 2020-07-25 16:42:55
言語 C++14
(gcc 13.3.0 + boost 1.87.0)
結果
WA  
実行時間 -
コード長 1,095 bytes
コンパイル時間 1,753 ms
コンパイル使用メモリ 170,204 KB
実行使用メモリ 6,824 KB
最終ジャッジ日時 2024-11-18 18:01:18
合計ジャッジ時間 2,693 ms
ジャッジサーバーID
（参考情報） 		judge4 / judge2
このコードへのチャレンジ
（要ログイン）
ファイルパターン 結果
other AC * 4 WA * 6
権限があれば一括ダウンロードができます

ソースコード

diff # 

#include <bits/stdc++.h>
using namespace std;
typedef long long ll;
int challenge = 100;

ll powmod(ll a, ll x, ll p){
    ll res=1;
    while(x>0){
        if(x%2 == 1){
            res = (res * a) % p;
            x--;
        }
        a = (a*a)%p;
        x >>= 1;
    }
    return res;
}

bool is_prime(ll n){
    if(n == 2) return(true);
    if(n == 1 || (n&1) == 0) return (false);
    
    mt19937 mt;
    ll s=1;
    ll d=n;
    do{
        d >>= 1;
    }while((d&1) == 0);

    for(int i=0;i<challenge;i++){
        ll a = 1 + mt() % (n-1);
        ll t = d;
        ll y = powmod(a,t,n);
        while(t != n-1 && y != 1 && y != n-1){
            y = powmod(y,2,n);
            t <<= 1;
        }
        if(y != n-1 && (t&1) == 0) return false;
    }
    
    return(true);
}

int main() {
    int n;
    cin >> n;
    vector<ll> x(n);
    for(int i=0;i<n;i++){
        cin >> x[i];
    }

    for(int i=0;i<n;i++){
        cout << x[i] << " ";
        if(is_prime(x[i])){
            cout << 1 << endl;
        }else{
            cout << 0 << endl;
        }
    }

    return(0);
}
0