結果

問題 No.1659 Product of Divisors
ユーザー taiga0629kyoprotaiga0629kyopro
提出日時 2021-08-27 22:01:13
言語 PyPy3
(7.3.15)
結果
AC  
実行時間 54 ms / 2,000 ms
コード長 614 bytes
コンパイル時間 188 ms
コンパイル使用メモリ 82,356 KB
実行使用メモリ 59,908 KB
最終ジャッジ日時 2024-11-21 02:23:44
合計ジャッジ時間 2,030 ms
ジャッジサーバーID
（参考情報） 		judge1 / judge3
このコードへのチャレンジ
（要ログイン）

テストケース

テストケース表示
入力 結果 実行時間
実行使用メモリ
testcase_00 AC 36 ms
52,628 KB
testcase_01 AC 40 ms
59,908 KB
testcase_02 AC 37 ms
52,976 KB
testcase_03 AC 34 ms
52,468 KB
testcase_04 AC 49 ms
59,424 KB
testcase_05 AC 45 ms
57,784 KB
testcase_06 AC 37 ms
59,172 KB
testcase_07 AC 38 ms
57,972 KB
testcase_08 AC 53 ms
58,612 KB
testcase_09 AC 53 ms
57,936 KB
testcase_10 AC 53 ms
57,592 KB
testcase_11 AC 54 ms
59,908 KB
testcase_12 AC 54 ms
58,812 KB
testcase_13 AC 36 ms
52,800 KB
testcase_14 AC 54 ms
58,224 KB
testcase_15 AC 38 ms
58,516 KB
testcase_16 AC 53 ms
59,308 KB
testcase_17 AC 38 ms
58,568 KB
testcase_18 AC 38 ms
58,736 KB
testcase_19 AC 37 ms
58,008 KB
testcase_20 AC 47 ms
58,908 KB
testcase_21 AC 50 ms
58,000 KB
testcase_22 AC 46 ms
57,968 KB
testcase_23 AC 51 ms
57,764 KB
testcase_24 AC 44 ms
58,640 KB
権限があれば一括ダウンロードができます

ソースコード

diff # 

from _collections import defaultdict
n,k=map(int,input().split())

dic=defaultdict(int)
p=2

nod=n
while p**2<=n:
    if nod%p==0:
        dic[p]+=1
        nod//=p
    else:
        p+=1
if nod>1:dic[nod]+=1


ans=1
mod=10**9+7

def cnb(n,r):
    if n<0 or r<0 or r>n:
        return 0
    res=1
    res2=1
    for i in range(1,r+1):
        res*=n-i+1
        res%=mod
        res2*=i
        res2%=mod
    ans=res*pow(res2,mod-2,mod)
    ans%=mod
    return ans

for p in dic:
    subans=0
    for a in range(dic[p]+1):
        subans+=cnb(k-1+a,a)
        subans%=mod
    ans*=subans
    ans%=mod

print(ans)
0