結果
| 問題 | No.1659 Product of Divisors |
| コンテスト | |
| ユーザー |
brthyyjp
|
| 提出日時 | 2021-08-29 12:36:09 |
| 言語 | PyPy3 (7.3.15) |
| 結果 |
AC
|
| 実行時間 | 57 ms / 2,000 ms |
| コード長 | 778 bytes |
| 記録 | |
| コンパイル時間 | 315 ms |
| コンパイル使用メモリ | 82,176 KB |
| 実行使用メモリ | 57,728 KB |
| 最終ジャッジ日時 | 2024-11-22 04:58:51 |
| 合計ジャッジ時間 | 2,490 ms |
|
ジャッジサーバーID (参考情報) |
judge3 / judge4 |
(要ログイン)
| ファイルパターン | 結果 |
|---|---|
| sample | AC * 2 |
| other | AC * 23 |
ソースコード
import math
def factorize(n):
d = {}
temp = int(math.sqrt(n))+1
for i in range(2, temp):
while n%i== 0:
n //= i
if i in d:
d[i] += 1
else:
d[i] = 1
if d == {}:
d[n] = 1
else:
if n in d:
d[n] += 1
elif n != 1:
d[n] =1
return d
def cmb2(n, r, mod):
res = 1
temp = 1
for k in range(1, r+1):
res *= (n-k+1)
temp *= k
res %= mod
temp %= mod
res *= pow(temp,(mod-2),mod)
res %= mod
return res
n, k = map(int, input().split())
if n == 1:
print(1)
exit()
mod = 10**9+7
D = factorize(n)
ans = 1
for p, q in D.items():
ans *= cmb2(q+k, q, mod)
ans %= mod
print(ans)
brthyyjp