結果

問題 No.577 Prime Powerful Numbers
ユーザー square1001square1001
提出日時 2017-08-16 16:18:30
言語 C++14
(gcc 13.2.0 + boost 1.83.0)
結果
WA  
(最新)
AC  
(最初)
実行時間 -
コード長 1,917 bytes
コンパイル時間 674 ms
コンパイル使用メモリ 70,444 KB
実行使用メモリ 4,380 KB
最終ジャッジ日時 2023-08-11 04:43:32
合計ジャッジ時間 2,319 ms
ジャッジサーバーID
（参考情報） 		judge12 / judge14
このコードへのチャレンジ
（要ログイン）

テストケース

テストケース表示
入力 結果 実行時間
実行使用メモリ
testcase_00 AC 4 ms
4,380 KB
testcase_01 WA -
testcase_02 AC 3 ms
4,380 KB
testcase_03 AC 79 ms
4,376 KB
testcase_04 AC 7 ms
4,380 KB
testcase_05 AC 308 ms
4,376 KB
testcase_06 AC 61 ms
4,380 KB
testcase_07 AC 340 ms
4,376 KB
testcase_08 AC 69 ms
4,376 KB
testcase_09 AC 65 ms
4,380 KB
testcase_10 AC 2 ms
4,376 KB
権限があれば一括ダウンロードができます

ソースコード

diff # 

#include <cmath>
#include <iostream>
using namespace std;
long long power(long long a, int b) {
	long long ret = 1;
	while (b) {
		if (b & 1) ret *= a;
		a *= a;
		b >>= 1;
	}
	return ret;
}
long long root(long long a, int b) {
	// return maximum x that x^b ≤ a (assume a ≤ 10^18)
	if (b == 1) return a;
	int l = 0, r = (int)(exp((log(10) * 18 + log(3)) / b));
	while (r - l > 1) {
		int m = (l + r) >> 1;
		if (power(m, b) > a) r = m;
		else l = m;
	}
	return l;
}
long long modmul(long long a, long long b, long long m) {
	long long ret = 0; a %= m; b %= m;
	while (a) {
		if (a & 1) {
			ret += b;
			if (ret >= m) ret -= m;
		}
		b += b;
		if (b >= m) b -= m;
		a >>= 1;
	}
	return ret;
}
long long modpowl(long long a, long long b, long long m) {
	long long ret = 1;
	while (b) {
		if (b & 1) ret = modmul(ret, a, m);
		a = modmul(a, a, m);
		b >>= 1;
	}
	return ret;
}
bool issprp(long long x, int a) {
	int s = 0; long long d = x - 1;
	while (!(d & 1)) d >>= 1, s++;
	long long z = modpowl(a, d, x);
	if (z == 1) return true;
	for (int i = 0; i < s; i++) {
		if (z == x - 1) return true;
		z = modmul(z, z, x);
	}
	return false;
}
bool isprime(long long x) {
	if (x == 2) return true;
	if (x <= 1 || x % 2 == 0) return false;
	int a[] = { 2, 3, 5, 7, 11, 13, 17, 19, 23 };
	for (int i = 0; i < 9; i++) {
		if (x == a[i]) return true;
		if (!issprp(x, a[i])) return false;
	}
	return true;
}
bool isprimepower(long long x) {
	long long mul = 3;
	for (int i = 1; mul <= x; i++) {
		long long r = root(x, i);
		if (power(r, i) == x && isprime(r)) return true;
		mul *= 3;
	}
	return false;
}
bool solve(long long x) {
	if (x == 2) return false;
	if (x % 2 == 0) return true;
	for (long long b = 2; b < x; b *= 2) {
		if (isprimepower(x - b)) return true;
	}
	return false;
}
int main() {
	int Q;
	cin >> Q;
	while (Q--) {
		long long N;
		cin >> N;
		cout << (solve(N) ? "Yes" : "No") << endl;
	}
	return 0;
}
0