結果
問題 | No.2751 429-like Number |
ユーザー | hiro1729 |
提出日時 | 2024-05-10 21:25:38 |
言語 | C++23 (gcc 12.3.0 + boost 1.83.0) |
結果 |
AC
|
実行時間 | 934 ms / 4,000 ms |
コード長 | 2,597 bytes |
コンパイル時間 | 2,978 ms |
コンパイル使用メモリ | 262,304 KB |
実行使用メモリ | 6,944 KB |
最終ジャッジ日時 | 2024-05-10 21:26:00 |
合計ジャッジ時間 | 18,786 ms |
ジャッジサーバーID (参考情報) |
judge3 / judge2 |
(要ログイン)
テストケース
テストケース表示入力 | 結果 | 実行時間 実行使用メモリ |
---|---|---|
testcase_00 | AC | 1 ms
6,812 KB |
testcase_01 | AC | 2 ms
6,944 KB |
testcase_02 | AC | 2 ms
6,944 KB |
testcase_03 | AC | 2 ms
6,940 KB |
testcase_04 | AC | 2 ms
6,940 KB |
testcase_05 | AC | 2 ms
6,944 KB |
testcase_06 | AC | 3 ms
6,944 KB |
testcase_07 | AC | 365 ms
6,940 KB |
testcase_08 | AC | 458 ms
6,944 KB |
testcase_09 | AC | 349 ms
6,940 KB |
testcase_10 | AC | 877 ms
6,940 KB |
testcase_11 | AC | 863 ms
6,940 KB |
testcase_12 | AC | 934 ms
6,940 KB |
testcase_13 | AC | 429 ms
6,944 KB |
testcase_14 | AC | 615 ms
6,944 KB |
testcase_15 | AC | 24 ms
6,940 KB |
testcase_16 | AC | 632 ms
6,944 KB |
testcase_17 | AC | 627 ms
6,944 KB |
testcase_18 | AC | 896 ms
6,944 KB |
testcase_19 | AC | 901 ms
6,940 KB |
testcase_20 | AC | 906 ms
6,944 KB |
testcase_21 | AC | 864 ms
6,940 KB |
testcase_22 | AC | 880 ms
6,940 KB |
testcase_23 | AC | 878 ms
6,944 KB |
testcase_24 | AC | 865 ms
6,944 KB |
testcase_25 | AC | 877 ms
6,940 KB |
testcase_26 | AC | 865 ms
6,940 KB |
testcase_27 | AC | 855 ms
6,940 KB |
コンパイルメッセージ
main.cpp: In function 'll _findprime(ll)': main.cpp:29:276: warning: control reaches end of non-void function [-Wreturn-type] 29 | ll _findprime(ll n){if(n%2==0)return 2;for(int c=1;c<n;c++){auto f=[&](ll a)->ll {return((bll)a*a+c)%n;};ll x=0,y=0,g=1;while(g==1){x=f(x);y=f(f(y));g=gcd(abs(x-y),n);}if(g==n)continue;if(MillerRabin(g))return g;else if(MillerRabin(n/g))return n/g;else return _findprime(g);}} | ^
ソースコード
#include <bits/stdc++.h> using namespace std; using ll = long long; using bll = __int128_t; // O(√n) verified bool is_prime(ll n){if(n<2)return false;if(n==2)return true;if(n%2==0)return false;for(ll i=3;i*i<=n;i+=2)if(n%i==0)return false;return true;} // O(√n) verified vector<ll>divisors(ll n){vector<ll>r;int i;for(i=1;i*i<n;i++)if(n%i==0){r.push_back(i);r.push_back(n/i);}if(i*i==n)r.push_back(i);sort(r.begin(),r.end());return r;} // O(√n) verified vector<pair<ll,ll>>prime_factorize(ll n){vector<pair<ll,ll>>r;int c=0;while(n%2==0){n/=2;c++;}r.push_back({2,c});for(ll p=3;p*p<=n;p+=2){if(n%p!=0)continue;int e=0;while(n%p==0){e++;n/=p;}r.push_back({p,e});}if(n!=1)r.push_back({n,1});return r;} // エラトステネスの篩 // O(n*log(log(n))) verified vector<bool>EraSieve(ll n){vector<bool>p(n+1,true);p[0]=p[1]=false;for(int i=2;i<=n;i++){if(!p[i])continue;for(int q=i*2;q<=n;q+=i){p[q]=false;}}return p;} // エラトステネスの区間篩 verified // a以上b未満のcに対してp[c - a] vector<bool>EraSecSieve(ll a,ll b){ll sq=(ll)ceil(sqrt(b));vector<bool>ps(sq,true);ps[0]=ps[1]=false;for(ll i=2;i<sq;i++){if(!ps[i])continue;for(ll j=i*i;j<sq;j+=i)ps[j]=false;}vector<bool>p(b-a,true);for(ll i=2;i<sq;i++){if(!ps[i])continue;ll k=max(i,(a+i-1)/i);for(ll j=k*i;j<b;j+=i)p[j-a]=false;}return p;} // PollardRhoまでverified template<class T>T modpow(T a,T n,T m){T r=1;while(n){if(n&1)r=r*a%m;a=a*a%m;n>>=1;}return r;} // 失敗率 1/(4^30) bool MillerRabin(ll n) {if(n==2)return true;if(n==1||(n>2&&n%2==0))return false;ll s=0,t=n-1;while(t%2==0){s++;t>>=1;}random_device gen;mt19937_64 mtrand(gen());for(int i=0;i<30;i++){ll a=(mtrand())%(n-1)+1;if(a%n==0)continue;bool ok=false;if(modpow<bll>(a,t,n)==1)ok=true;else for(int i=0;i<s;i++)if(modpow<bll>(a,modpow<bll>(2,i,1ull<<63)*t,n)==n-1){ok=true;break;}if(!ok)return false;}return true;} ll _findprime(ll n){if(n%2==0)return 2;for(int c=1;c<n;c++){auto f=[&](ll a)->ll {return((bll)a*a+c)%n;};ll x=0,y=0,g=1;while(g==1){x=f(x);y=f(f(y));g=gcd(abs(x-y),n);}if(g==n)continue;if(MillerRabin(g))return g;else if(MillerRabin(n/g))return n/g;else return _findprime(g);}} vector<pair<ll,ll>>PollardRho(ll n){vector<pair<ll,ll>>r;while(!MillerRabin(n)&&n>1){ll p=_findprime(n),s=0;while(n%p==0){n/=p;s++;}r.push_back({p,s});}if(n>1)r.push_back({n,1});return r;} int main() { int q; ll a; cin >> q; while (cin >> a) { auto p = PollardRho(a); int c = 0; for (auto [i, j]: p) c += j; if (c == 3) { cout << "Yes\n"; } else { cout << "No\n"; } } }