結果
問題 |
No.2365 Present of good number
|
ユーザー |
|
提出日時 | 2023-06-30 23:13:38 |
言語 | PyPy3 (7.3.15) |
結果 |
WA
|
実行時間 | - |
コード長 | 1,332 bytes |
コンパイル時間 | 201 ms |
コンパイル使用メモリ | 82,560 KB |
実行使用メモリ | 57,856 KB |
最終ジャッジ日時 | 2024-07-07 11:08:31 |
合計ジャッジ時間 | 2,440 ms |
ジャッジサーバーID (参考情報) |
judge1 / judge5 |
(要ログイン)
ファイルパターン | 結果 |
---|---|
sample | AC * 2 |
other | AC * 17 WA * 22 |
ソースコード
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) * k2, MOD) three = pow(3, pow(2, amari // 2, MOD) * k3, MOD) else: two = pow(2, pow(2, (amari + 1) // 2, MOD) * k3, MOD) three = pow(3, pow(2, amari // 2, MOD) * k2, MOD) # print(amari, two, three, two * three, pow(2, amari // 2, MOD)) N = (two * three) % MOD break else: N = c print(N % MOD)