結果

問題 No.1659 Product of Divisors
ユーザー qibqib
提出日時 2022-12-23 02:03:29
言語 PyPy3
(7.3.15)
結果
AC  
実行時間 61 ms / 2,000 ms
コード長 621 bytes
コンパイル時間 233 ms
コンパイル使用メモリ 82,432 KB
実行使用メモリ 58,112 KB
最終ジャッジ日時 2024-04-29 04:08:18
合計ジャッジ時間 2,562 ms
ジャッジサーバーID
（参考情報） 		judge1 / judge2
このコードへのチャレンジ
（要ログイン）

テストケース

テストケース表示
入力 結果 実行時間
実行使用メモリ
testcase_00 AC 44 ms
52,096 KB
testcase_01 AC 49 ms
57,472 KB
testcase_02 AC 42 ms
52,096 KB
testcase_03 AC 42 ms
51,968 KB
testcase_04 AC 58 ms
57,600 KB
testcase_05 AC 56 ms
57,856 KB
testcase_06 AC 45 ms
57,856 KB
testcase_07 AC 49 ms
57,984 KB
testcase_08 AC 61 ms
57,856 KB
testcase_09 AC 61 ms
57,984 KB
testcase_10 AC 60 ms
57,600 KB
testcase_11 AC 60 ms
57,728 KB
testcase_12 AC 60 ms
57,984 KB
testcase_13 AC 40 ms
51,584 KB
testcase_14 AC 61 ms
57,856 KB
testcase_15 AC 46 ms
57,344 KB
testcase_16 AC 45 ms
58,112 KB
testcase_17 AC 46 ms
57,472 KB
testcase_18 AC 46 ms
57,856 KB
testcase_19 AC 46 ms
57,472 KB
testcase_20 AC 54 ms
57,984 KB
testcase_21 AC 57 ms
57,344 KB
testcase_22 AC 55 ms
57,728 KB
testcase_23 AC 59 ms
57,344 KB
testcase_24 AC 53 ms
57,344 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