結果
| 問題 |
No.1659 Product of Divisors
|
| コンテスト | |
| ユーザー |
prussian_coder
|
| 提出日時 | 2022-03-27 21:05:10 |
| 言語 | PyPy3 (7.3.15) |
| 結果 |
RE
|
| 実行時間 | - |
| コード長 | 821 bytes |
| コンパイル時間 | 492 ms |
| コンパイル使用メモリ | 82,408 KB |
| 実行使用メモリ | 847,832 KB |
| 最終ジャッジ日時 | 2024-11-06 13:56:07 |
| 合計ジャッジ時間 | 3,652 ms |
|
ジャッジサーバーID (参考情報) |
judge4 / judge5 |
(要ログイン)
| ファイルパターン | 結果 |
|---|---|
| sample | AC * 2 |
| other | AC * 8 RE * 5 MLE * 1 -- * 9 |
ソースコード
N,K=map(int,input().split())
ans=1
mod=10**9+7
if N==1:
print(1)
exit()
def factorization(n):
arr = []
temp = n
for i in range(2, int(-(-n**0.5//1))+1):
if temp%i==0:
cnt=0
while temp%i==0:
cnt+=1
temp //= i
arr.append([i, cnt])
if temp!=1:
arr.append([temp, 1])
if arr==[]:
arr.append([n, 1])
return arr
M=K+100
fact=[1 for _ in range(M+1)]
factinv=[1 for _ in range(M+1)]
a=1
b=1
p=10**9+7
for i in range(M):
a*=(i+1)
a%=p
b*=pow(i+1,p-2,p)
b%=p
fact[i+1],factinv[i+1]=a,b
def comb(n,k,p):
return fact[n]*(factinv[k]*factinv[n-k])%p
L=factorization(N)
for l in L:
a=l[1]
s=0
for i in range(a+1):
s+=comb(i+K-1,K-1,mod)
ans*=s
ans%=mod
print(ans)
prussian_coder