結果
問題 | No.577 Prime Powerful Numbers |
ユーザー | ciel |
提出日時 | 2017-10-14 19:54:50 |
言語 | C (gcc 12.3.0) |
結果 |
WA
|
実行時間 | - |
コード長 | 1,066 bytes |
コンパイル時間 | 164 ms |
コンパイル使用メモリ | 32,896 KB |
実行使用メモリ | 6,824 KB |
最終ジャッジ日時 | 2024-11-17 18:48:07 |
合計ジャッジ時間 | 6,867 ms |
ジャッジサーバーID (参考情報) |
judge2 / judge4 |
(要ログイン)
テストケース
テストケース表示入力 | 結果 | 実行時間 実行使用メモリ |
---|---|---|
testcase_00 | WA | - |
testcase_01 | WA | - |
testcase_02 | AC | 4 ms
6,816 KB |
testcase_03 | AC | 459 ms
6,816 KB |
testcase_04 | AC | 31 ms
6,816 KB |
testcase_05 | TLE | - |
testcase_06 | AC | 179 ms
6,816 KB |
testcase_07 | TLE | - |
testcase_08 | AC | 390 ms
6,816 KB |
testcase_09 | AC | 389 ms
6,816 KB |
testcase_10 | AC | 1 ms
6,816 KB |
ソースコード
#include <stdio.h> #include <stdlib.h> #include <stdbool.h> #include <math.h> typedef unsigned long long ull; ull mul(ull x,ull y,ull m){ ull z=0; for(;y;y>>=1){ if(y&1)z=(z+x)%m; x=(x+x)%m; } return z; } ull powmod(ull x,ull y,ull m){ ull z=1; for(;y;y>>=1){ if(y&1)z=mul(z,x,m); x=mul(x,x,m); } return z; } ull gcd(ull x,ull y){return y?gcd(y,x%y):x;} bool miller_rabin(ull n){ if(n==2)return true; if(n==1||n%2==0)return false; ull d=n-1,s=0,a=1; for(;d%2==0;d/=2)s+=1; for(int k=5;k--;){ //todo for(a++;gcd(a,n)!=1;a++); //todo ull r=powmod(a,d,n); if(r==1||r==n-1)continue; int t=s; for(;t;t--){ r=powmod(r,2,n); if(r==n-1)break; } if(!t)return false; } return true; } void solve(ull n){ if(n%2==0){puts(n>2?"Yes":"No");return;} for(ull m=2;m<n;m*=2){ ull z=n-m; for(int k=1;k<61;k++){ ull r=k==1?z:expl(logl(n+1)/k)+1e-6; if(m+powmod(r,k,1LL<<62)!=n)continue; if(miller_rabin(r)){puts("Yes");return;} } } puts("No"); } int main(){ ull t,n; for(scanf("%llu",&t);t--;solve(n))scanf("%llu",&n); }