結果

問題 No.2795 Perfect Number
ユーザー tama04tama04
提出日時 2024-06-28 22:16:57
言語 C++14
(gcc 13.3.0 + boost 1.87.0)
結果
AC  
実行時間 2 ms / 2,000 ms
コード長 1,323 bytes
コンパイル時間 1,835 ms
コンパイル使用メモリ 171,392 KB
実行使用メモリ 5,376 KB
最終ジャッジ日時 2024-06-28 22:17:00
合計ジャッジ時間 2,435 ms
ジャッジサーバーID
（参考情報） 		judge4 / judge5
このコードへのチャレンジ
（要ログイン）
ファイルパターン 結果
sample AC * 3
other AC * 35
権限があれば一括ダウンロードができます

ソースコード

diff # 

#include <bits/stdc++.h>
using namespace std;
typedef long long ll;
typedef long double ld;
template<class T> inline bool chmin(T& a, T b) { if (a > b) { a = b; return true; } else { return false; } }
template<class T> inline bool chmax(T& a, T b) { if (a < b) { a = b; return true; } else { return false; } }
#define rep(i, s, n) for(int i = (s); i < n; i++) 
#define reps(i, s, n) for(int i = (s); i <= n; i++) 
#define rrep(i, n, s) for(int i = n; i >= (s); i--)
#define all(x) (x).begin(), (x).end()
#define print(V) for (int z = 0; z < (V).size(); z++) { cout << V[z] << " "; } cout << endl;
#define MOD 998244353
#define MOD2 1000000007
#define INF 1LL << 60
#define PI acos(-1)

ll culc(ll k) {
	ll res = 1;
	rep(i, 0, k - 1) {
		res *= 2;
	}

	return res;
}

bool is_prime(ll n) {
	if (n == 1) return false;
	set<ll> S;
	for (ll i = 2; i * i <= n; i++) {
		if (n % i == 0) {
			S.insert(i);
		}
	}
	if (S.empty()) {
		return true;
	}
	else {
		return false;
	}
}


int main(void) {
	ll n;
	cin >> n;

	if (n % 2 != 0) {
		cout << "No" << endl;
		return 0;
	}

	ll k = 2, p = -1;
	bool flag = false;
	while (p < n) {
		p = culc(k);

		if (n == p * (p * 2 - 1) && is_prime(p * 2 - 1)) {
			flag = true;
			break;
		}
		k++;
	}

	if (flag) {
		cout << "Yes" << endl;
	}
	else {
		cout << "No" << endl;
	}

	return 0;
}
0