結果

問題 No.8030 ミラー・ラビン素数判定法のテスト
ユーザー kacho65535
提出日時 2020-11-15 13:43:16
言語 C++14
(gcc 13.3.0 + boost 1.89.0)
結果
WA  
実行時間 -
コード長 2,161 bytes
記録
記録タグの例:
初AC ショートコード 純ショートコード 純主流ショートコード 最速実行時間
コンパイル時間 1,662 ms
コンパイル使用メモリ 171,592 KB
実行使用メモリ 6,820 KB
最終ジャッジ日時 2024-11-18 18:48:28
合計ジャッジ時間 4,654 ms
ジャッジサーバーID
（参考情報） 		judge4 / judge1
このコードへのチャレンジ
（要ログイン）
ファイルパターン 結果
other AC * 4 WA * 6
権限があれば一括ダウンロードができます

ソースコード

diff #
raw source code 

#include <bits/stdc++.h>
using namespace std;
using ll = long long;
#define mod 1000000007ll
#define loop(i, n) for (int i = 0; i < n; i++)
#define flagcount(bit) __builtin_popcount(bit)
#define flag(x) (1ll << x)
#define flagadd(bit, x) bit |= flag(x)
#define flagpop(bit, x) bit &= ~flag(x)
#define flagon(bit, i) bit &flag(i)
#define flagoff(bit, i) !(bit & (1ll << i))
#define all(v) v.begin(), v.end()
#define putout(a) cout << a << '\n'
template <typename T>
bool chmax(T &a, const T &b)
{
    if (a < b)
    {
        a = b; // aをbで更新
        return true;
    }
    return false;
}
template <typename T>
bool chmin(T &a, const T &b)
{
    if (a > b)
    {
        a = b; // aをbで更新
        return true;
    }
    return false;
}
void AS(ll X, ll L, ll R) { assert(L <= X && X <= R); }
template <typename T>
bool primejudge(T n)
{
    if (n < 2)
        return false;
    else if (n == 2)
        return true;
    else if (n % 2 == 0)
        return false;
    double sqrtn = sqrt(n);
    for (T i = 3; i < sqrtn + 1; i++)
    {
        if (n % i == 0)
        {
            return false;
        }
        i++;
    }
    return true;
}
//modpow
long long modp(ll a, ll n, ll p)
{
    long long res = 1;
    while (n > 0)
    {
        if (n & 1)
            res = res * a % p;
        a = a * a % p;
        n >>= 1;
    }
    return res;
}
//素数ならtrue,合成数ならfalseを返す
bool MR(ll X)
{
    if (X < 10000)
        return primejudge(X);
    if (X % 2 == 0)
        return false;
    ll Z = X - 1;
    ll s = 0, t = 0;
    while (Z % 2 == 0)
    {
        s++;
        Z /= 2;
    }
    t = Z;
    vector<ll> A = {2, 3, 5, 7, 11, 13, 17, 19, 23, 29, 31, 37};
    for (auto a : A)
    {
        if (modp(a, t, X) == 1)
            return true;
        loop(i, s)
        {
            if (modp(a, modp(2, i, X) * t, X) == X - 1)
                return true;
        }
    }
    return false;
}
int main()
{
    ll Q;
    cin >> Q;
    loop(i, Q)
    {
        ll X;
        cin >> X;
        if (MR(X))
            cout << X << " " << 1 << endl;
        else
            cout << X << " " << 0 << endl;
    }
    return 0;
}
0