結果
問題 | No.2365 Present of good number |
ユーザー |
|
提出日時 | 2023-06-30 23:16:37 |
言語 | PyPy3 (7.3.15) |
結果 |
AC
|
実行時間 | 32 ms / 2,000 ms |
コード長 | 895 bytes |
コンパイル時間 | 280 ms |
コンパイル使用メモリ | 82,188 KB |
実行使用メモリ | 53,740 KB |
最終ジャッジ日時 | 2024-07-07 11:11:33 |
合計ジャッジ時間 | 2,447 ms |
ジャッジサーバーID (参考情報) |
judge2 / judge5 |
(要ログイン)
ファイルパターン | 結果 |
---|---|
sample | AC * 2 |
other | AC * 39 |
ソースコード
#int(input()) #map(int, input().split()) #list(map(int, input().split())) N, K = map(int, input().split()) mod = 10 ** 9 + 7 def prime_factorize(n): a = [] while n % 2 == 0: a.append(2) n //= 2 f = 3 while f * f <= n: if n % f == 0: a.append(f) n //= f else: f += 2 if n != 1: a.append(n) return a def f(A): p = prime_factorize(A) for i in range(len(p)): p[i] += 1 return p while K: u = f(N) K -= 1 N = 1 for y in u: N *= y if u[-1] <= 4: break if K != 0: p3 = u.count(3) p2 = 2 * u.count(4) a, b = divmod(K, 2) p2 = (p2 * pow(2, a, mod-1)) % (mod-1) p3 = (p3 * pow(2, a, mod-1)) % (mod-1) if b: p2, p3 = 2*p3, p2 ans = (pow(2, p2, mod) * pow(3, p3, mod)) % mod else: ans = N print(ans % mod)