結果

問題 No.1659 Product of Divisors
ユーザー qibqib
提出日時 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
このコードへのチャレンジ
（要ログイン）

テストケース

テストケース表示
入力 結果 実行時間
実行使用メモリ
testcase_00 AC 37 ms
52,160 KB
testcase_01 AC 41 ms
57,916 KB
testcase_02 AC 38 ms
52,436 KB
testcase_03 AC 38 ms
52,632 KB
testcase_04 AC 52 ms
58,644 KB
testcase_05 AC 48 ms
57,968 KB
testcase_06 AC 41 ms
58,724 KB
testcase_07 AC 45 ms
58,028 KB
testcase_08 AC 57 ms
58,916 KB
testcase_09 AC 57 ms
58,192 KB
testcase_10 AC 56 ms
57,840 KB
testcase_11 AC 56 ms
58,416 KB
testcase_12 AC 56 ms
58,312 KB
testcase_13 AC 36 ms
52,980 KB
testcase_14 AC 56 ms
58,756 KB
testcase_15 AC 43 ms
58,048 KB
testcase_16 AC 40 ms
58,256 KB
testcase_17 AC 41 ms
58,436 KB
testcase_18 AC 39 ms
58,968 KB
testcase_19 AC 41 ms
59,176 KB
testcase_20 AC 48 ms
58,440 KB
testcase_21 AC 52 ms
59,440 KB
testcase_22 AC 49 ms
58,876 KB
testcase_23 AC 58 ms
58,696 KB
testcase_24 AC 50 ms
59,252 KB
権限があれば一括ダウンロードができます

ソースコード

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