結果

問題 No.577 Prime Powerful Numbers
コンテスト
ユーザー vjudge1
提出日時 2025-04-25 21:48:54
言語 C++23
(gcc 13.3.0 + boost 1.87.0)
結果
WA  
実行時間 -
コード長 1,146 bytes
コンパイル時間 3,080 ms
コンパイル使用メモリ 280,960 KB
実行使用メモリ 23,076 KB
最終ジャッジ日時 2025-04-25 21:49:09
合計ジャッジ時間 4,865 ms
ジャッジサーバーID
（参考情報） 		judge1 / judge2
このコードへのチャレンジ
（要ログイン）
ファイルパターン 結果
sample AC * 1
other AC * 3 WA * 7
権限があれば一括ダウンロードができます

ソースコード

diff # 

#include<bits/stdc++.h>
#define int long long
using namespace std;
int fastmul(int x,int y,int mod)
{
	int z=(double)x*y/mod;
	return ((x*y-z*mod)%mod+mod)%mod;
}
int fastpow(int x,int y,int mod)
{
	int ans=1;
	for(;y;y>>=1,x=fastmul(x,x,mod))
	{
		if(y&1)
		{
			ans=fastmul(ans,x,mod);
		}
	}
	return ans;
}
bool miller_rabin(int n)
{
	for(int i=1;i<=50;i++)
	{
		int x=rand()*rand()%n+1;
		if(fastpow(x,n-1,n)!=1)return 0;
	}
	return 1;
}
vector<int>prm;
set<int>s;
int q,n,f[1000005];
signed main()
{
	srand(time(0));
	for(int i=2,w;i<=1000000;i++)
	{
		if(!f[i])
		{
			prm.push_back(i);w=i;
			for(;;)
			{
				s.insert(w);
				if(w>1000000000000000000ll/i)break;
				w*=i;
			}
		}
		for(int j=0;j<prm.size()&&i*prm[j]<=1000000;j++)
		{
			f[i*prm[j]]=true;
			if(i%prm[j]==0)break;
		}
	}
	cin>>q;
	while(q--)
	{
		cin>>n;int f=0;
		if(n==2){cout<<"No\n";continue;}
		if(n%2==0){cout<<"Yes\n";continue;}
		for(int i=2,sq;i<n-1;i*=2)
		{
			sq=sqrt(n-i);
			while(sq*sq<n-i)sq++;
			while(sq*sq>n-i)sq--;
			if(s.count(n-i)||miller_rabin(n-i)||(sq*sq==n-i&&miller_rabin(sq))){f=1;break;}
		}
		if(f)cout<<"Yes\n";
		else cout<<"No\n";
	}
}
0