結果

問題 No.8030 ミラー・ラビン素数判定法のテスト
ユーザー konishi
提出日時 2023-11-07 09:35:35
言語 C++23
(gcc 13.3.0 + boost 1.87.0)
結果
WA  
実行時間 -
コード長 2,029 bytes
コンパイル時間 2,678 ms
コンパイル使用メモリ 245,668 KB
実行使用メモリ 6,944 KB
最終ジャッジ日時 2024-09-25 23:16:53
合計ジャッジ時間 3,719 ms
ジャッジサーバーID
（参考情報） 		judge3 / judge2
このコードへのチャレンジ
（要ログイン）
ファイルパターン 結果
other AC * 4 WA * 6
権限があれば一括ダウンロードができます

ソースコード

diff #
raw source code 

#include <bits/stdc++.h>
//#include<atcoder/all>
//using namespace atcoder;
using namespace std;
typedef long long ll;
typedef unsigned long long ull;
typedef long double ld;

#define pii pair<int,int>
#define pll pair<ll,ll>

#define rep(i,num,n) for(int i=num;i<(int)(n);i++) //for_loop
#define REP(i,n) for(int i=0;i<(int)(n);i++)
#define rrep(i,num,n) for(int i=num-1;i>=(int)(n);i--) //reverse_for>
#define in(x,a,b) (a<=x && x<b)//範囲
#define reo return 0
#define all(v) v.begin(),v.end()

#define INFI 1010000000
#define INFL 4000000000000000000
//#define ten(x) ((int)1e##x)
//#define tenll(x) ((ll)1e##x)
#define PI 3.14159265358979
//#define mint modint1000000007
#define mint modint998244353

bool chmin(ll &a,ll b){if(a>b){a=b;return true;}return false;}
bool chmax(ll &a,ll b){if(a<b){a=b;return true;}return false;}
bool chmin(int &a,int b){if(a>b){a=b;return true;}return false;}
bool chmax(int &a,int b){if(a<b){a=b;return true;}return false;}

int dx[]={0,0,1,-1};int dy[]={1,-1,0,0};
//int dx[]={1,1,1,0,0,-1,-1,-1},dy[]={1,0,-1,1,-1,1,0,-1};
//----------------------------
ll modpow(ll a,ll b,ll n){
    ll t=1;
    while(b){
        if(b&1){
            t=(t*a)%n;
        }
        a=(a*a)%n;
        b>>=1;
    }
    return t;
}
ll modmul(ll a,ll b,ll n){
    return (a*b)%n;
    
}
const uint64_t bases[] = {2,325,9375,28178,450775,9780504,1795265022};
bool miller_rabin(ll n){
    if(n==2)return true;
    if(n<2 ||(n&1)==0)return false;
    ll n1=n-1,d=n-1;
    ll s=0;
    for(;(d&1)==0;d>>=1)s++;//s=dの二進桁数
    for(auto &base:bases){
        ll a=base;
        if(a>=n)
        {
            a%=n;
            if(a==0)continue;
        }
    
    ll t=modpow(a,d,n);
    if(t==1)continue;
    for(ll j=1;t!=n1;j++){
        if(j>=s)return false;
    t=modmul(t,t,n);

    }
    }
    return true;
}

int  main(){
    int n;cin>>n;
    rep(i,0,n){
        ll a;cin>>a;cout<<a<<" ";
        //ll a=i;cout<<a<<" ";
        if(miller_rabin(a))cout<<1<<endl;
        else cout<<0<<endl;
    }
}
0