結果
問題 | No.577 Prime Powerful Numbers |
ユーザー |
|
提出日時 | 2025-04-27 22:09:12 |
言語 | C++17 (gcc 13.3.0 + boost 1.87.0) |
結果 |
WA
|
実行時間 | - |
コード長 | 1,666 bytes |
コンパイル時間 | 2,250 ms |
コンパイル使用メモリ | 195,664 KB |
実行使用メモリ | 7,844 KB |
最終ジャッジ日時 | 2025-04-27 22:09:18 |
合計ジャッジ時間 | 4,828 ms |
ジャッジサーバーID (参考情報) |
judge5 / judge3 |
(要ログイン)
ファイルパターン | 結果 |
---|---|
sample | AC * 1 |
other | AC * 9 WA * 1 |
ソースコード
#include <bits/stdc++.h> #define int unsigned __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 maxx){ int res=1; while(z){ if(res>maxx)return res; if(z&1)res*=d; 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 Isp(int n) { if (n<3||n%2==0)return n==2; int a=n-1,b=0,tt=114; 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=v*v%n;} if(j>b)return 0; } return 1; } 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,x); 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;i<=x;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; }