結果

問題 No.577 Prime Powerful Numbers
コンテスト
ユーザー はまやんはまやん
提出日時 2017-10-14 14:33:57
言語 C++14
(gcc 13.3.0 + boost 1.87.0)
結果
TLE  
実行時間 -
コード長 5,016 bytes
コンパイル時間 1,739 ms
コンパイル使用メモリ 169,284 KB
実行使用メモリ 16,840 KB
最終ジャッジ日時 2024-11-17 18:40:35
合計ジャッジ時間 29,408 ms
ジャッジサーバーID
（参考情報） 		judge3 / judge4
このコードへのチャレンジ
（要ログイン）
ファイルパターン 結果
sample AC * 1
other AC * 2 TLE * 8
権限があれば一括ダウンロードができます

ソースコード

diff # 

#include<bits/stdc++.h>
#define rep(i,a,b) for(int i=a;i<b;i++)
#define rrep(i,a,b) for(int i=a;i>=b;i--)
#define fore(i,a) for(auto &i:a)
#pragma GCC optimize ("-O3")
using namespace std; void _main(); int main() { cin.tie(0); ios::sync_with_stdio(false); _main(); }
//---------------------------------------------------------------------------------------------------
#ifdef _MSC_VER
#pragma push_macro("long")
#undef long
#ifdef _WIN32
inline unsigned int __builtin_ctz(unsigned int x) { unsigned long r; _BitScanForward(&r, x); return r; }
inline unsigned int __builtin_clz(unsigned int x) { unsigned long r; _BitScanReverse(&r, x); return 31 - r; }
inline unsigned int __builtin_ffs(unsigned int x) { unsigned long r; return _BitScanForward(&r, x) ? r + 1 : 0; }
inline unsigned int __builtin_popcount(unsigned int x) { return __popcnt(x); }
inline unsigned int hidword(unsigned long long x) { return static_cast<unsigned int>(x >> 32); }
inline unsigned int lodword(unsigned long long x) { return static_cast<unsigned int>(x & 0xFFFFFFFF); }
inline unsigned long long __builtin_ctzll(unsigned long long x) { return lodword(x) ? __builtin_ctz(lodword(x)) : __builtin_ctz(hidword(x)) + 32; }
inline unsigned long long __builtin_clzll(unsigned long long x) { return hidword(x) ? __builtin_clz(hidword(x)) : __builtin_clz(lodword(x)) + 32; }
inline unsigned long long __builtin_ffsll(unsigned long long x) { return lodword(x) ? __builtin_ffs(lodword(x)) : hidword(x) ? __builtin_ffs(hidword(x)) + 32 : 0; }
inline unsigned long long __builtin_popcountll(unsigned long long x) { return __builtin_popcount(lodword(x)) + __builtin_popcount(hidword(x)); }
#endif // _WIN32
#pragma pop_macro("long")
#endif // _MSC_VER
typedef long long ll;
typedef __int128 int128_t;
//typedef ll int128_t;
int128_t powmod(int128_t x, int128_t y, int128_t p) { // O(log y)
    assert(0 <= x and x < p);
    assert(0 <= y);
    int128_t z = 1;
    for (int128_t i = 1; i <= y; i <<= 1) {
        if (y & i) z = z * x % p;
        x = x * x % p;
    }
    return z;
}
int K = 10;
int getrand(int l, int r) { // [l, r]
    static uint32_t y = time(NULL);
    y ^= (y << 13); y ^= (y >> 17);
    y ^= (y << 5);
    return y % (r - l + 1) + l;
}
ll getrandll(ll l, ll r) { // [l,r]
    ll a = getrand(0, 999999999);
    ll b = getrand(0, 999999999);

    ll x = a * 1000000000LL + b;
    assert(0 <= x);

    return x % (r - l + 1) + l;
}
ll pr[15] = { 2,3,5,7,11,13,17,19,23,29,31,37,41,43,47 };
bool isprime(ll n) {
    if (n == 2) return true;
    if (n == 1 || n % 2 == 0) return false;
    const ll d = (n - 1) >> __builtin_ctzll(n - 1);
    rep(_, 0, 15) {
        //ll a = getrandll(1, n - 2);
        ll a = pr[_];
        if (a == n) { return true; }
        if (a >= n) break;
        int128_t t = d;
        int128_t y = powmod(a, t, n);
        while (t != n - 1 and y != 1 and y != n - 1) {
            y = y * y % n;
            t *= 2;
        }
        if (y != n - 1 and t % 2 == 0) return false;
    }
    return true;
}
/*---------------------------------------------------------------------------------------------------
　　　　　　　　　　　 ∧＿∧  
　　　　　 ∧＿∧ 　（´<_｀ ）　 Welcome to My Coding Space!
　　　　 （ ´_ゝ`）　/　 ⌒i     
　　　　／　　　＼　 　  |　|     
　　　 /　　 /￣￣￣￣/　　|  
　 ＿_(__ﾆつ/　    ＿/ .| .|＿＿＿＿  
　 　　　＼/＿＿＿＿/　（u　⊃  
---------------------------------------------------------------------------------------------------*/




typedef long long ll; typedef long double ld;
#define INF 1LL<<60
ld lginf = -1;
inline ll mul(ll a, ll b) {
    if (a == 0 or b == 0) return 0;
    if (a*b < a or a*b < b) return INF;
    return a*b;
}
inline ll fastpow(ll x, ll n) {
    ll ret = 1;
    while (0 < n) {
        if ((n % 2) == 0) x = mul(x,x), n >>= 1;
        else ret = mul(ret,x), --n;
    }
    return ret;
}









ll N;
#define INF 1LL<<60
#define yes "Yes"
#define no "No"
//---------------------------------------------------------------------------------------------------
int isprimePow(ll n) {
    int k = 1;
    int up = n;
    for (ll q = 2; q < n; q *= 2) {
        ll ng = 0, ok = up;
        while (ng + 1 != ok) {
            ll x = (ng + ok) / 2;
            if (n <= fastpow(x, k)) ok = x;
            else ng = x;
        }

        if (fastpow(ok, k) == n and isprime(ok)) return 1;
        k++;
        up = ok;
    }
    return 0;
}
//---------------------------------------------------------------------------------------------------
string solve() {
    cin >> N;
    if (N <= 2) return no;

    if (N % 2 == 0) {
        if (4 <= N) return yes;
        else return no;
    }
    
    for (ll q = 2; q < N; q *= 2) if (isprimePow(N - q)) return yes;
    return no;
}
//---------------------------------------------------------------------------------------------------
void _main() {
    int Q; cin >> Q;
    rep(q, 0, Q) cout << solve() << endl;
}
0