結果
問題 | No.8030 ミラー・ラビン素数判定法のテスト |
ユーザー |
|
提出日時 | 2021-06-10 18:48:58 |
言語 | C++17 (gcc 13.3.0 + boost 1.87.0) |
結果 |
AC
|
実行時間 | 302 ms / 9,973 ms |
コード長 | 1,786 bytes |
コンパイル時間 | 2,254 ms |
コンパイル使用メモリ | 196,896 KB |
最終ジャッジ日時 | 2025-01-22 05:16:09 |
ジャッジサーバーID (参考情報) |
judge2 / judge3 |
(要ログイン)
ファイルパターン | 結果 |
---|---|
other | AC * 10 |
ソースコード
#include<bits/stdc++.h> using namespace std; #define rep(i,n) for(ll i=0;i<n;i++) #define repl(i,l,r) for(ll i=(l);i<(r);i++) #define per(i,n) for(ll i=(n)-1;i>=0;i--) #define perl(i,r,l) for(ll i=r-1;i>=l;i--) #define fi first #define se second #define pb push_back #define ins insert #define pqueue(x) priority_queue<x,vector<x>,greater<x>> #define all(x) (x).begin(),(x).end() #define CST(x) cout<<fixed<<setprecision(x) #define rev(x) reverse(x); using ll=long long; using vl=vector<ll>; using vvl=vector<vector<ll>>; using pl=pair<ll,ll>; using vpl=vector<pl>; using vvpl=vector<vpl>; const ll MOD=1000000007; const ll MOD9=998244353; const int inf=1e9+10; const ll INF=4e18; const ll dy[8]={1,0,-1,0,1,1,-1,-1}; const ll dx[8]={0,1,0,-1,1,-1,1,-1}; template <typename T> inline bool chmax(T &a, T b) { return ((a < b) ? (a = b, true) : (false)); } template <typename T> inline bool chmin(T &a, T b) { return ((a > b) ? (a = b, true) : (false)); } struct fast_factorial{ using i128=__int128_t; vector<ll> wit={2, 325, 9375, 28178, 450775, 9780504, 1795265022}; ll modpow(ll a,ll b,ll m){ ll ret=1,now=a; while(b){ if(b&1)ret=i128(ret)*now%m; now=i128(now)*now%m; b>>=1; } return ret; } bool miller_rabin(ll p){ if(p==2)return true; if(p==1||p%2==0)return false; bool ans=true; ll s=0,d=p-1; while(d%2==0){ d/=2;s++; } for(auto a:wit){ if(a%p==0)continue; bool ret=true; ll x=modpow(a,d,p); if(x==1)ret=false; rep(i,s){ if(x==p-1)ret=false; x=i128(x)*x%p; } if(ret)ans=false; } return ans; } }; int main(){ ll q;cin >> q; fast_factorial f; while(q--){ ll m;cin >> m; cout << m <<" " << f.miller_rabin(m) <<endl; } }