結果
問題 | No.2954 Calculation of Exponentiation |
ユーザー |
![]() |
提出日時 | 2024-11-08 22:23:44 |
言語 | C++14 (gcc 13.3.0 + boost 1.87.0) |
結果 |
AC
|
実行時間 | 3 ms / 2,000 ms |
コード長 | 2,649 bytes |
コンパイル時間 | 1,774 ms |
コンパイル使用メモリ | 167,348 KB |
実行使用メモリ | 5,248 KB |
最終ジャッジ日時 | 2024-11-08 22:23:48 |
合計ジャッジ時間 | 2,811 ms |
ジャッジサーバーID (参考情報) |
judge2 / judge1 |
(要ログイン)
ファイルパターン | 結果 |
---|---|
sample | AC * 3 |
other | AC * 28 |
ソースコード
#include <bits/stdc++.h>using namespace std;#define CPP_STR(x) CPP_STR_I(x)#define CPP_CAT(x, y) CPP_CAT_I(x, y)#define CPP_STR_I(args...) #args#define CPP_CAT_I(x, y) x##y#define ASSERT(expr...) assert((expr))using i8 = int8_t;using u8 = uint8_t;using i16 = int16_t;using u16 = uint16_t;using i32 = int32_t;using u32 = uint32_t;using i64 = int64_t;using u64 = uint64_t;using f32 = float;using f64 = double;// }}}constexpr i64 INF = 1'010'000'000'000'000'017LL;constexpr i64 MOD = 998244353LL;constexpr f64 EPS = 1e-12;constexpr f64 PI = 3.14159265358979323846;#define M5 100007#define M9 1000000000#define F first#define S second// util {{{#define FOR(i, start, end) for (i64 i = (start), CPP_CAT(i, xxxx_end) = (end); i < CPP_CAT(i, xxxx_end); ++i)#define REP(i, n) FOR(i, 0, n)#define all(x) (x).begin(), (x).end()#define ll long long int#define VI vector<ll>#define VVI vector<VI>#define ISD true#define debug(x) \if (ISD) \cout << #x << ": " << x << endltemplate <typename T, typename U, typename Comp = less<>>bool chmax(T &xmax, const U &x, Comp comp = {}){if (comp(xmax, x)){xmax = x;return true;}return false;}template <typename T, typename U, typename Comp = less<>>bool chmin(T &xmin, const U &x, Comp comp = {}){if (comp(x, xmin)){xmin = x;return true;}return false;}long long GCD(long long a, long long b){if (b == 0)return a;elsereturn GCD(b, a % b);}long long safe_pow(long long a, long long b){long long res = 1;for (long long i = 0; i < b; i++){double dres = res;dres *= a;if (dres > 2e18){return 2e18;}res *= a;}return res;}void solve(){double A, B;cin >> A >> B;if (A == 1 || B == 0){cout << "Yes" << endl;return;}ll T = 10000;ll LA = A * T;ll LB = B * T;if (B < 0){ll G = GCD(LA, T);if (LA / G != 1){cout << "No" << endl;}else{cout << "Yes" << endl;}return;}if (LA / T * T != LA){cout << "No" << endl;return;}LA /= T;ll G = GCD(LB, T);T /= G;REP(i, 1000000){ll NX = safe_pow(i, T);if (NX == LA){cout << "Yes" << endl;return;}if (NX > LA){cout << "No" << endl;return;}}}int main(){solve();}