結果
問題 | No.2954 Calculation of Exponentiation |
ユーザー |
👑 ![]() |
提出日時 | 2024-11-08 21:51:48 |
言語 | C++23 (gcc 13.3.0 + boost 1.87.0) |
結果 |
AC
|
実行時間 | 2 ms / 2,000 ms |
コード長 | 2,452 bytes |
コンパイル時間 | 3,833 ms |
コンパイル使用メモリ | 258,136 KB |
実行使用メモリ | 5,248 KB |
最終ジャッジ日時 | 2024-11-08 21:51:53 |
合計ジャッジ時間 | 3,675 ms |
ジャッジサーバーID (参考情報) |
judge5 / judge4 |
(要ログイン)
ファイルパターン | 結果 |
---|---|
sample | AC * 3 |
other | AC * 28 |
ソースコード
#include <bits/stdc++.h>using namespace std;#define FOR(i,m,n) for(int i=(m);i<(n);++i)#define REP(i,n) FOR(i,0,n)using ll = long long;constexpr int INF = 0x3f3f3f3f;constexpr long long LINF = 0x3f3f3f3f3f3f3f3fLL;constexpr double EPS = 1e-8;constexpr int MOD = 998244353;// constexpr int MOD = 1000000007;constexpr int DY4[]{1, 0, -1, 0}, DX4[]{0, -1, 0, 1};constexpr int DY8[]{1, 1, 0, -1, -1, -1, 0, 1};constexpr int DX8[]{0, -1, -1, -1, 0, 1, 1, 1};template <typename T, typename U>inline bool chmax(T& a, U b) { return a < b ? (a = b, true) : false; }template <typename T, typename U>inline bool chmin(T& a, U b) { return a > b ? (a = b, true) : false; }struct IOSetup {IOSetup() {std::cin.tie(nullptr);std::ios_base::sync_with_stdio(false);std::cout << fixed << setprecision(20);}} iosetup;template <bool GETS_ONLY_PRIME>std::vector<int> prime_sieve(const int n) {std::vector<int> smallest_prime_factor(n + 1), prime;std::iota(smallest_prime_factor.begin(), smallest_prime_factor.end(), 0);for (int i = 2; i <= n; ++i) {if (smallest_prime_factor[i] == i) [[unlikely]] prime.emplace_back(i);for (const int p : prime) {if (i * p > n || p > smallest_prime_factor[i]) break;smallest_prime_factor[i * p] = p;}}return GETS_ONLY_PRIME ? prime : smallest_prime_factor;}template <typename T>std::vector<std::pair<T, int>> prime_factorization(T n) {std::vector<std::pair<T, int>> res;for (T i = 2; i * i <= n; ++i) {if (n % i == 0) [[unlikely]] {int exponent = 0;for (; n % i == 0; n /= i) {++exponent;}res.emplace_back(i, exponent);}}if (n > 1) res.emplace_back(n, 1);return res;}pair<ll, ll> normalize(string s) {s.erase(prev(s.end(), 5));const ll num = stoll(s), den = 10000, g = gcd(num, den);return {num / g, den / g};}int main() {string a, b; cin >> a >> b;if (b == "0.0000") {cout << "Yes\n";return 0;}auto [num_a, den_a] = normalize(a);auto [num_b, den_b] = normalize(b);if (num_b < 0) {swap(num_a, den_a);num_b = -num_b;}if (den_a != 1) {cout << "No\n";return 0;}if (num_a == 1) {cout << (den_b == 1 ? "Yes\n" : "No\n");return 0;}int ex = 0;for (const int v : prime_factorization(num_a) | ranges::views::values) {ex = gcd(ex, v);}cout << (ex % den_b == 0 ? "Yes\n" : "No\n");return 0;}