結果
問題 | No.3030 ミラー・ラビン素数判定法のテスト |
ユーザー | kopricky |
提出日時 | 2018-11-08 05:42:40 |
言語 | C++11 (gcc 11.4.0) |
結果 |
AC
|
実行時間 | 396 ms / 9,973 ms |
コード長 | 2,936 bytes |
コンパイル時間 | 1,544 ms |
コンパイル使用メモリ | 159,984 KB |
実行使用メモリ | 5,248 KB |
最終ジャッジ日時 | 2024-11-16 23:09:56 |
合計ジャッジ時間 | 2,729 ms |
ジャッジサーバーID (参考情報) |
judge1 / judge2 |
(要ログイン)
テストケース
テストケース表示入力 | 結果 | 実行時間 実行使用メモリ |
---|---|---|
testcase_00 | AC | 2 ms
5,248 KB |
testcase_01 | AC | 1 ms
5,248 KB |
testcase_02 | AC | 1 ms
5,248 KB |
testcase_03 | AC | 1 ms
5,248 KB |
testcase_04 | AC | 201 ms
5,248 KB |
testcase_05 | AC | 185 ms
5,248 KB |
testcase_06 | AC | 52 ms
5,248 KB |
testcase_07 | AC | 52 ms
5,248 KB |
testcase_08 | AC | 52 ms
5,248 KB |
testcase_09 | AC | 396 ms
5,248 KB |
ソースコード
#include <bits/stdc++.h> #define INF 1000000005 #define MOD 1000000007 #define EPS 1e-10 #define rep(i,n) for(int i=0;i<(int)(n);++i) #define rrep(i,n) for(int i=(int)(n)-1;i>=0;--i) #define srep(i,s,t) for(int i=(int)(s);i<(int)(t);++i) #define each(a,b) for(auto& (a): (b)) #define along long(v) (v).begin(),(v).end() #define len(v) (int)(v).size() #define zip(v) sort(along long(v)),v.erase(unique(along long(v)),v.end()) #define cmx(x,y) x=max(x,y) #define cmn(x,y) x=min(x,y) #define fi first #define se second #define pb push_back #define show(x) cout<<#x<<" = "<<(x)<<endl #define spair(p) cout<<#p<<": "<<p.fi<<" "<<p.se<<endl #define sar(a,n) cout<<#a<<":";rep(pachico,n)cout<<" "<<a[pachico];cout<<endl #define svec(v) cout<<#v<<":";rep(pachico,v.size())cout<<" "<<v[pachico];cout<<endl #define svecp(v) cout<<#v<<":";each(pachico,v)cout<<" {"<<pachico.first<<":"<<pachico.second<<"}";cout<<endl #define sset(s) cout<<#s<<":";each(pachico,s)cout<<" "<<pachico;cout<<endl #define smap(m) cout<<#m<<":";each(pachico,m)cout<<" {"<<pachico.first<<":"<<pachico.second<<"}";cout<<endl using namespace std; typedef pair<int,int> P; typedef vector<int> vi; typedef vector<vi> vvi; typedef vector<long long> vl; typedef vector<vl> vvl; typedef vector<double> vd; typedef vector<P> vp; typedef vector<string> vs; const int MAX_N = 100005; //動作はint型のみ保証 //long long型の二乗とかが怪しい(__int128使ったらいける) // const vector<int> numset = {2,7,61}; // < 4,759,123,141 ≒ 4×10^9 // const vector<int> numset = {2,3,5,7,11,13,17}; // < 341,550,071,728,321 ≒ 3*10^14 const vector<int> numset = {2,3,5,7,11,13,17,19,23,29,31,37}; // 2^64までいける // long long mod_pow(long long x, long long k, long long mod){ // long long res = 1; // while(k){ // if(k & 1) res = res * x % mod; // x = x * x % mod; // k >>= 1; // } // return res; // } // long long型の素数判定の場合 long long mod_pow(__int128 x, long long k, long long mod){ __int128 res = 1; while(k){ if(k & 1) res = res * x % mod; x = x * x % mod; k >>= 1; } return res; } // check if n is prime bool check(long long n){ if(n < 2) return false; long long d = n - 1, s = 0; while(d % 2 == 0){ d /= 2; s++; } for(int a : numset){ if(a == n) return true; long long res = mod_pow(a, d, n); if(res != 1){ bool ok = true; for(int r = 0; r < s; r++){ if(mod_pow(res, 1LL << r, n) == n-1){ ok = false; break; } } if(ok) return false; } } return true; } int main() { cin.tie(0); ios::sync_with_stdio(false); int n; cin >> n; rep(i,n){ long long x; cin >> x; cout << x << " " << check(x) << "\n"; } return 0; }