結果

問題 No.2954 Calculation of Exponentiation
ユーザー ArcAkiArcAki
提出日時 2024-11-09 02:35:33
言語 PyPy3
(7.3.15)
結果
AC  
実行時間 84 ms / 2,000 ms
コード長 1,295 bytes
コンパイル時間 354 ms
コンパイル使用メモリ 82,048 KB
実行使用メモリ 76,436 KB
最終ジャッジ日時 2024-11-09 02:35:37
合計ジャッジ時間 3,551 ms
ジャッジサーバーID
（参考情報） 		judge2 / judge3
このコードへのチャレンジ
（要ログイン）

テストケース

テストケース表示
入力 結果 実行時間
実行使用メモリ
testcase_00 AC 57 ms
62,336 KB
testcase_01 AC 59 ms
61,824 KB
testcase_02 AC 59 ms
61,952 KB
testcase_03 AC 59 ms
61,952 KB
testcase_04 AC 57 ms
61,696 KB
testcase_05 AC 59 ms
61,952 KB
testcase_06 AC 58 ms
61,952 KB
testcase_07 AC 57 ms
62,336 KB
testcase_08 AC 58 ms
62,080 KB
testcase_09 AC 75 ms
62,080 KB
testcase_10 AC 57 ms
61,952 KB
testcase_11 AC 59 ms
61,696 KB
testcase_12 AC 58 ms
61,952 KB
testcase_13 AC 59 ms
61,824 KB
testcase_14 AC 59 ms
62,208 KB
testcase_15 AC 84 ms
76,436 KB
testcase_16 AC 77 ms
76,288 KB
testcase_17 AC 58 ms
61,952 KB
testcase_18 AC 59 ms
61,568 KB
testcase_19 AC 80 ms
76,200 KB
testcase_20 AC 57 ms
61,952 KB
testcase_21 AC 58 ms
61,952 KB
testcase_22 AC 63 ms
61,888 KB
testcase_23 AC 59 ms
61,952 KB
testcase_24 AC 58 ms
62,208 KB
testcase_25 AC 83 ms
76,176 KB
testcase_26 AC 59 ms
62,208 KB
testcase_27 AC 58 ms
61,824 KB
testcase_28 AC 60 ms
61,952 KB
testcase_29 AC 59 ms
61,568 KB
testcase_30 AC 60 ms
61,952 KB
権限があれば一括ダウンロードができます

ソースコード

diff # 

from sys import stdin
from random import randint
from collections import deque, defaultdict as dd
from copy import deepcopy
input = stdin.readline
MOD = 998244353
INF = 1 << 60
xor = randint(100, INF)


def extended_gcd(a, b):
    if b == 0:
        return a, 1, 0
    else:
        g, x, y = extended_gcd(b, a % b)
        return g, y, x - (a // b) * y


def mod_inverse(a, m):
    _, x, _ = extended_gcd(a, m)
    return (x % m + m) % m


def fast_mod_pow(x, p):
    res = 1
    t = x
    z = p
    while z > 0:
        if z % 2 == 1:
            res = (res * t)
        t = (t * t)
        z //= 2
    return res


def main():
    a, b = map(float, input().split())
    div = 10000
    a = int(a*div)
    b = int(b*div)
    if not b:
        print("Yes")
        return
    g1 = extended_gcd(a, div)[0]
    if not (div//g1==1 and b > 0 or a//g1==1 and b < 0):
        print("No")
        return
    if b < 0:
        a = div//g1
    else:
        a //= g1
    l = 0
    r = 1<<30
    g2 = div//extended_gcd(b, div)[0]
    while l+1 < r:
        m = (l+r)//2
        if fast_mod_pow(m, g2) <= a:
            l = m
        else:
            r = m
    if fast_mod_pow(l, g2)==a or fast_mod_pow(r, g2)==a:
        print("Yes")
    else:
        print("No")


if __name__ == "__main__":
    main()
0