結果

問題 No.2365 Present of good number
ユーザー yabit
提出日時 2023-06-30 23:40:16
言語 PyPy3
(7.3.15)
結果
AC  
実行時間 45 ms / 2,000 ms
コード長 1,340 bytes
コンパイル時間 240 ms
コンパイル使用メモリ 82,216 KB
実行使用メモリ 59,384 KB
最終ジャッジ日時 2024-07-07 11:31:25
合計ジャッジ時間 2,796 ms
ジャッジサーバーID
（参考情報） 		judge4 / judge5
このコードへのチャレンジ
（要ログイン）
ファイルパターン 結果
sample AC * 2
other AC * 39
権限があれば一括ダウンロードができます

ソースコード

diff # 

def factorization(n):  # 素因数分解 [値, かけられた数]
    arr = []
    temp = n
    for i in range(2, int(-(-(n**0.5) // 1)) + 1):
        if temp % i == 0:
            cnt = 0
            while temp % i == 0:
                cnt += 1
                temp //= i
            arr.append([i, cnt])
    if temp != 1:
        arr.append([temp, 1])
    if arr == []:
        arr.append([n, 1])
    return arr


N, K = map(int, input().split())
MOD = 10**9 + 7
for i in range(K):
    # print(i, N, factorization(N))
    ft = factorization(N)
    c = 1
    flag = 1
    for n, k in ft:
        if n != 2 and n != 3:
            flag = 0
        c *= (n + 1) ** k
    if flag:
        amari = K - i
        if len(ft) == 1:
            if ft[0][0] == 2:
                ft.append([3, 0])
            else:
                ft.insert(0, [2, 0])
        k2, k3 = ft[0][1], ft[1][1]
        if amari % 2 == 0:
            two = pow(2, pow(2, amari // 2, MOD-1) * k2, MOD)
            three = pow(3, pow(2, amari // 2, MOD-1) * k3, MOD)
        else:
            two = pow(2, pow(2, (amari + 1) // 2, MOD-1) * k3, MOD)
            three = pow(3, pow(2, amari // 2, MOD-1) * k2, MOD)
        # print(amari, two, three, two * three, pow(2, amari // 2, MOD))
        N = (two * three) % MOD
        break
    else:
        N = c
print(N % MOD)
0