結果
問題 | No.8030 ミラー・ラビン素数判定法のテスト |
ユーザー | fumofumofuni |
提出日時 | 2021-06-10 18:34:11 |
言語 | C++17 (gcc 12.3.0 + boost 1.83.0) |
結果 |
WA
|
実行時間 | - |
コード長 | 1,781 bytes |
コンパイル時間 | 2,005 ms |
コンパイル使用メモリ | 203,020 KB |
実行使用メモリ | 5,248 KB |
最終ジャッジ日時 | 2024-11-30 04:40:25 |
合計ジャッジ時間 | 4,797 ms |
ジャッジサーバーID (参考情報) |
judge5 / judge3 |
(要ログイン)
テストケース
テストケース表示入力 | 結果 | 実行時間 実行使用メモリ |
---|---|---|
testcase_00 | WA | - |
testcase_01 | AC | 2 ms
5,248 KB |
testcase_02 | AC | 2 ms
5,248 KB |
testcase_03 | AC | 2 ms
5,248 KB |
testcase_04 | AC | 239 ms
5,248 KB |
testcase_05 | AC | 241 ms
5,248 KB |
testcase_06 | AC | 246 ms
5,248 KB |
testcase_07 | AC | 252 ms
5,248 KB |
testcase_08 | AC | 248 ms
5,248 KB |
testcase_09 | AC | 243 ms
5,248 KB |
ソースコード
#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)break; 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; } }