結果
| 問題 |
No.1659 Product of Divisors
|
| コンテスト | |
| ユーザー |
tktk_snsn
|
| 提出日時 | 2021-08-27 23:34:31 |
| 言語 | Python3 (3.13.1 + numpy 2.2.1 + scipy 1.14.1) |
| 結果 |
AC
|
| 実行時間 | 134 ms / 2,000 ms |
| コード長 | 705 bytes |
| コンパイル時間 | 202 ms |
| コンパイル使用メモリ | 12,416 KB |
| 実行使用メモリ | 10,880 KB |
| 最終ジャッジ日時 | 2024-11-21 05:27:08 |
| 合計ジャッジ時間 | 2,703 ms |
|
ジャッジサーバーID (参考情報) |
judge2 / judge1 |
(要ログイン)
| ファイルパターン | 結果 |
|---|---|
| sample | AC * 2 |
| other | AC * 23 |
ソースコード
def binomial_coefficient(n, r):
if r < 0 or r > n:
return 0
res = 1
div = 1
r = min(r, n - r)
for i in range(r):
res = res * (n - i) % mod
div = div * (i + 1) % mod
return res * pow(div, mod-2, mod) % mod
def prime_factorization(n): # 素因数分解
res = []
for i in range(2, int(n ** 0.5) + 1):
if n % i == 0:
cnt = 0
while n % i == 0:
cnt += 1
n //= i
yield cnt
if n > 1:
yield 1
comb = binomial_coefficient
mod = 10 ** 9 + 7
N, K = map(int, input().split())
ans = 1
for f in prime_factorization(N):
ans *= comb(K + f, K)
ans %= mod
print(ans)
tktk_snsn