結果
問題 | No.1659 Product of Divisors |
ユーザー |
![]() |
提出日時 | 2021-09-14 09:03:01 |
言語 | PyPy3 (7.3.15) |
結果 |
WA
|
実行時間 | - |
コード長 | 800 bytes |
コンパイル時間 | 214 ms |
コンパイル使用メモリ | 81,920 KB |
実行使用メモリ | 252,240 KB |
最終ジャッジ日時 | 2024-06-26 21:32:07 |
合計ジャッジ時間 | 7,222 ms |
ジャッジサーバーID (参考情報) |
judge3 / judge2 |
(要ログイン)
ファイルパターン | 結果 |
---|---|
sample | AC * 2 |
other | AC * 6 WA * 2 RE * 15 |
ソースコード
N,K=map(int, input().split()) mod=10**9+7 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=0 for i in range(b+1): cnt+=cmb(K-1+i,K-1) ans*=cnt ans%=mod print(ans)