結果

問題 No.1659 Product of Divisors
ユーザー qib
提出日時 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
権限があれば一括ダウンロードができます

ソースコード

diff #

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)
0