結果
| 問題 | No.1659 Product of Divisors |
| コンテスト | |
| ユーザー |
timi
|
| 提出日時 | 2021-09-14 09:41:48 |
| 言語 | PyPy3 (7.3.15) |
| 結果 |
AC
|
| 実行時間 | 232 ms / 2,000 ms |
| コード長 | 899 bytes |
| 記録 | |
| コンパイル時間 | 314 ms |
| コンパイル使用メモリ | 82,132 KB |
| 実行使用メモリ | 252,536 KB |
| 最終ジャッジ日時 | 2024-06-26 22:12:28 |
| 合計ジャッジ時間 | 6,517 ms |
|
ジャッジサーバーID (参考情報) |
judge2 / judge1 |
(要ログイン)
| ファイルパターン | 結果 |
|---|---|
| sample | AC * 2 |
| other | AC * 23 |
ソースコード
N,K=map(int, input().split())
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
n=10**6
inv_t=[0]+[1]
for i in range(2,n):
inv_t+=[inv_t[mod%i]*(mod-int(mod/i))%mod]
kai=[1,1]
rev_kai=[1,inv_t[1]]
for i in range(2,n):
kai.append(kai[-1]*i%mod)
rev_kai.append(rev_kai[-1]*inv_t[i]%mod)
def cmb(n,r):
return kai[n]*rev_kai[r]*rev_kai[n-r]%mod
A=factorization(N)
#print(A)
ans=1
for a,b in A:
cnt=1
d=1
#素因数aをi回使い、分配するとき
for i in range(1,b+1):
d*=(K+i-1)
d//=(i)
cnt+=d
ans*=cnt
ans%=mod
print(ans)
timi