結果
問題 |
No.577 Prime Powerful Numbers
|
ユーザー |
|
提出日時 | 2025-04-27 20:45:46 |
言語 | C++17 (gcc 13.3.0 + boost 1.87.0) |
結果 |
WA
|
実行時間 | - |
コード長 | 1,740 bytes |
コンパイル時間 | 2,377 ms |
コンパイル使用メモリ | 199,560 KB |
実行使用メモリ | 40,072 KB |
最終ジャッジ日時 | 2025-04-27 20:45:53 |
合計ジャッジ時間 | 6,929 ms |
ジャッジサーバーID (参考情報) |
judge3 / judge1 |
(要ログイン)
ファイルパターン | 結果 |
---|---|
sample | AC * 1 |
other | AC * 1 WA * 9 |
ソースコード
#include <bits/stdc++.h> #define int long long 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; } } using namespace std; using namespace CZS_Office; int n,prime[6000005],prcnt; bool isp[10000005],isCZS[10000005],dp[10000005]; vector<int>ntt; inline void init(){ for(int i=2;i<=10000000;++i){ if(!isp[i])prime[++prcnt]=i; for(int j=1;j<=prcnt&&prime[j]*i<=10000000;++j){ isp[prime[j]*i]=1; if(i%prime[j]==0)break; } } for(int i=1;i<=prcnt;++i) for(int res=prime[i];res<=10000000;res*=prime[i])isCZS[res]=1; for(int i=1;i<=10000000;++i)if(isCZS[i])ntt.push_back(i); return; } unordered_map<int,bool>mp; inline int ksm(int d,int z){ int res=1; while(z){ if(z&1){ if(1ull*d*res>1e18)return -1; res*=d; } d*=d;z>>=1; } return 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 1; if(res>x)R=mid-1; else L=mid+1; } return 0; } inline bool chk(int x){ for(int i=2;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<=61;++j,i*=2){ if(chk(x-i))return 1; } return 0; } for(auto i:ntt)if(isCZS[x-i])return 1; return 0; } signed main(){ // freopen("number.in","r",stdin); // freopen("number.out","w",stdout); init();int Q=read(); while(Q--){ int n=read(); if(liliang(n))cout<<"Yes\n"; else cout<<"No\n"; } return 0; }