結果
問題 |
No.1659 Product of Divisors
|
ユーザー |
![]() |
提出日時 | 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)