結果
問題 |
No.1659 Product of Divisors
|
ユーザー |
|
提出日時 | 2021-08-27 22:03:11 |
言語 | Python3 (3.13.1 + numpy 2.2.1 + scipy 1.14.1) |
結果 |
AC
|
実行時間 | 133 ms / 2,000 ms |
コード長 | 752 bytes |
コンパイル時間 | 79 ms |
コンパイル使用メモリ | 12,672 KB |
実行使用メモリ | 10,752 KB |
最終ジャッジ日時 | 2024-11-21 02:27:02 |
合計ジャッジ時間 | 2,470 ms |
ジャッジサーバーID (参考情報) |
judge4 / judge2 |
(要ログイン)
ファイルパターン | 結果 |
---|---|
sample | AC * 2 |
other | AC * 23 |
ソースコード
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 if N == 1: print(1) exit() p = factorization(N) L = len(p) ans = 1 def H(n,r): tmp = 1 for i in range(r): tmp *= n+r-1-i tmp %= mod tmp *= pow(r-i,-1,mod) tmp %= mod return tmp for i in range(L): p_num = p[i][1] S_p = 1 for j in range(1,p_num+1): S_p += H(K,j) S_p %= mod ans *= S_p ans %= mod print(ans%mod)