結果
| 問題 | No.1659 Product of Divisors |
| コンテスト | |
| ユーザー |
|
| 提出日時 | 2022-12-23 02:03:29 |
| 言語 | PyPy3 (7.3.15) |
| 結果 |
AC
|
| 実行時間 | 58 ms / 2,000 ms |
| コード長 | 621 bytes |
| 記録 | |
| コンパイル時間 | 266 ms |
| コンパイル使用メモリ | 82,124 KB |
| 実行使用メモリ | 59,440 KB |
| 最終ジャッジ日時 | 2024-11-18 03:51:43 |
| 合計ジャッジ時間 | 2,644 ms |
|
ジャッジサーバーID (参考情報) |
judge3 / judge2 |
(要ログイン)
| ファイルパターン | 結果 |
|---|---|
| sample | AC * 2 |
| other | AC * 23 |
ソースコード
MOD = 10 ** 9 + 7
n, k = map(int, input().split())
m = 100
inv = [1 for _ in range(m + 1)]
for x in range(2, m + 1):
inv[x] = (- (MOD // x) * inv[MOD % x]) % MOD
def rep_comb(n, r):
ans = 1
for x in range(r):
ans = (ans * (n + r - 1 - x)) % MOD
ans = (ans * inv[x + 1]) % MOD
return ans
cur = n
divs = {}
d = 2
while d * d <= n:
while cur % d == 0:
divs[d] = divs.get(d, 0) + 1
cur //= d
d += 1
if cur > 1:
divs[cur] = divs.get(cur, 0) + 1
ans = 1
for _, v in divs.items():
cnt = 0
for r in range(v + 1):
cnt = (cnt + rep_comb(k, r)) % MOD
ans = (ans * cnt) % MOD
print(ans)