結果

問題 No.2954 Calculation of Exponentiation
ユーザー leafirbyleafirby
提出日時 2024-11-08 22:48:23
言語 C++14
(gcc 12.3.0 + boost 1.83.0)
結果
AC  
実行時間 2 ms / 2,000 ms
コード長 3,224 bytes
コンパイル時間 2,151 ms
コンパイル使用メモリ 177,960 KB
実行使用メモリ 5,248 KB
最終ジャッジ日時 2024-11-08 22:48:39
合計ジャッジ時間 3,134 ms
ジャッジサーバーID
(参考情報)
judge2 / judge1
このコードへのチャレンジ
(要ログイン)

テストケース

テストケース表示
入力 結果 実行時間
実行使用メモリ
testcase_00 AC 2 ms
5,248 KB
testcase_01 AC 2 ms
5,248 KB
testcase_02 AC 2 ms
5,248 KB
testcase_03 AC 2 ms
5,248 KB
testcase_04 AC 2 ms
5,248 KB
testcase_05 AC 2 ms
5,248 KB
testcase_06 AC 2 ms
5,248 KB
testcase_07 AC 2 ms
5,248 KB
testcase_08 AC 2 ms
5,248 KB
testcase_09 AC 2 ms
5,248 KB
testcase_10 AC 2 ms
5,248 KB
testcase_11 AC 2 ms
5,248 KB
testcase_12 AC 2 ms
5,248 KB
testcase_13 AC 2 ms
5,248 KB
testcase_14 AC 2 ms
5,248 KB
testcase_15 AC 2 ms
5,248 KB
testcase_16 AC 2 ms
5,248 KB
testcase_17 AC 2 ms
5,248 KB
testcase_18 AC 2 ms
5,248 KB
testcase_19 AC 2 ms
5,248 KB
testcase_20 AC 2 ms
5,248 KB
testcase_21 AC 2 ms
5,248 KB
testcase_22 AC 2 ms
5,248 KB
testcase_23 AC 2 ms
5,248 KB
testcase_24 AC 2 ms
5,248 KB
testcase_25 AC 2 ms
5,248 KB
testcase_26 AC 2 ms
5,248 KB
testcase_27 AC 2 ms
5,248 KB
testcase_28 AC 2 ms
5,248 KB
testcase_29 AC 2 ms
5,248 KB
testcase_30 AC 2 ms
5,248 KB
権限があれば一括ダウンロードができます

ソースコード

diff #

#include <bits/stdc++.h>
#pragma GCC target("avx2")
#pragma GCC optimize("O3")
#pragma GCC optimize("unroll-loops")
using namespace std;
#define ll long
#define rep(i, m, n) for (ll i = m; i < n; i++)
#define rrep(i, m, n) for (ll i = m - 1; i >= n; i--)
#define pb(n) push_back(n)
#define UE(N) N.erase(unique(N.begin(), N.end()), N.end());
#define Sort(n) sort(n.begin(), n.end())
#define Rev(n) reverse(n.begin(), n.end())
#define HOUT(S) cout << setprecision(50) << S << endl
#define pint pair<ll, ll>
#define tint tuple<ll, ll, ll>
#define endl "\n"
#define mod 1000000007LL
#define MAX 100010LL
#define chmax(a, b) a = (((a) < (b)) ? (b) : (a))
#define chmin(a, b) a = (((a) > (b)) ? (b) : (a))

template <typename T> vector<T> make_v(size_t a) { return vector<T>(a); }

template <typename T, typename... Ts> auto make_v(size_t a, Ts... ts) {
    return vector<decltype(make_v<T>(ts...))>(a, make_v<T>(ts...));
}

template <typename T, typename U, typename... V>
typename enable_if<is_same<T, U>::value != 0>::type fill_v(U& u, const V... v) {
    u = U(v...);
}

template <typename T, typename U, typename... V>
typename enable_if<is_same<T, U>::value == 0>::type fill_v(U& u, const V... v) {
    for (auto& e : u) fill_v<T>(e, v...);
}

ll gcd(ll a, ll b) {
    if(b > a) swap(a, b);
    if(b == 0) return a;
    else return gcd(b, a % b);
}

map<ll, int> prime_factor(ll n) {
    map<ll, int> res;
    for (ll i = 2; i * i <= n; i++) {
        while (n % i == 0) {
            res[i]++;
            n /= i;
        }
    }
    if (n != 1) res[n] = 1;
    return res;
}

int main() {
    cin.tie(nullptr);
    ios::sync_with_stdio(false);
    string A, B;
    cin >> A >> B;
    Rev(A);
    if(B == "0.0000") cout << "Yes" << endl; 
    else if(A.substr(0, 5) != "0000." && B[0] != '-') cout << "No" << endl;
    else if(B[0] != '-'){
        Rev(B);
        ll a = 0;
        rep(i, 5, A.size()) a += (A[i] - '0') * pow(10, i - 5);
        bool f = true;
        ll b = 0;
        rep(i, 0, 4) b += (B[i] - '0') * pow(10, i);
        rep(i, 5, B.size()) b += (B[i] - '0') * pow(10, i);
        ll g = gcd(b, 10000);
        ll m = 10000 / g;
        auto mp = prime_factor(a);
        
        for(auto v : mp) {
            if(v.second % m != 0) f = false;
        }
        if(f) cout << "Yes" << endl;
        else cout << "No" << endl;
        //cout << a << " " << b << " " << g << " " << m << endl;
    }
    else { 
        if(A.size() != 6 || A[5] != '0') cout << "No" << endl;
        else {
            ll a = 0;
            rep(i, 0, 4) a += (A[i] - '0') * pow(10, i);
            bool f = true;
            if(10000 % a != 0) f = false;
            else {
                a = 10000 / a;
                Rev(B);
                ll b = 0;
                rep(i, 0, 4) b += (B[i] - '0') * pow(10, i);
                rep(i, 5, B.size() - 1) b += (B[i] - '0') * pow(10, i);
                ll g = gcd(b, 10000);
                ll m = 10000 / g;
                auto mp = prime_factor(a);

                for(auto v : mp) {
                    if(v.second % m != 0) f = false;
                }
            }
            if(f) cout << "Yes" << endl;
            else cout << "No" << endl;
        }
    }
}
0