結果
問題 | No.577 Prime Powerful Numbers |
ユーザー |
|
提出日時 | 2025-04-27 22:00:40 |
言語 | C++17 (gcc 13.3.0 + boost 1.87.0) |
結果 |
WA
|
実行時間 | - |
コード長 | 1,858 bytes |
コンパイル時間 | 2,126 ms |
コンパイル使用メモリ | 197,316 KB |
実行使用メモリ | 7,844 KB |
最終ジャッジ日時 | 2025-04-27 22:00:50 |
合計ジャッジ時間 | 9,458 ms |
ジャッジサーバーID (参考情報) |
judge1 / judge5 |
(要ログイン)
ファイルパターン | 結果 |
---|---|
sample | AC * 1 |
other | AC * 7 WA * 1 TLE * 2 |
ソースコード
#include <bits/stdc++.h> #define int __int128 namespace CZS_Office{ inline int read(){ int x=0,f=1; char ch=getchar(); while(ch<'0'||ch>'9'){ if(ch=='-')f=-1; ch=getchar(); } while(ch>='0'&&ch<='9')x=x*10+ch-'0',ch=getchar(); return x*f; } inline void write(int x){ if(x<0){ putchar('-'); x=-x; } if(x>9)write(x/10); putchar(x%10+'0'); return; } inline int ksm(int d,int z){ int res=1; while(z){ if(z&1){ if(d*res>1e18)return 2e18; res*=d; } if(d*d>1e18)d=2e18; d*=d;z>>=1; } return res; } inline int shr(int d,int z,int p){ int res=1; while(z){ if(z&1){ (res*=d)%=p; } (d*=d)%=p;z>>=1; } return res; } } using namespace std; using namespace CZS_Office; const __int128 ONE = 1;mt19937 rnd(time(0)); bool millerRabin(int n) { if (n < 3 || n % 2 == 0) return n == 2; int a = n - 1 , b = 0 , tt = 11; while (a % 2 == 0) a /= 2 , ++ b; for(int i = 1 ; i <= tt ; i ++) { int x=rnd()%(n-2)+2,v=shr(x,a,n); if (v == 1) continue;int j; for(j = 1 ; j <= b ; j ++) { if (v == n - 1) break; v = (int)v * v % n; } if (j > b) return 0; } return 1; } inline bool Isp(int res){ return millerRabin(res); } inline int czs(int x,int z){ int L=1,R=x; while(L<=R){ int mid=(L+R)>>1; int res=ksm(mid,z); if(res==x)return Isp(mid); if(res>x)R=mid-1; else L=mid+1; } return 0; } inline bool chk(int x){ for(int i=1;i<=60;++i) if(czs(x,i))return 1; return 0; } inline bool liliang(int x){ if(x%2==0){ if(x>=4)return 1; else return 0; } else{ for(int i=2,j=1;i<=x&&j<60;++j,i*=2){ if(chk(x-i))return 1; } return 0; } return 0; } signed main(){ // freopen("number.in","r",stdin); // freopen("number.out","w",stdout); int Q=read(); while(Q--){ int n=read(); if(liliang(n))cout<<"Yes\n"; else cout<<"No\n"; } return 0; }