結果

問題 No.1659 Product of Divisors
ユーザー prussian_coderprussian_coder
提出日時 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
このコードへのチャレンジ
（要ログイン）

テストケース

テストケース表示
入力 結果 実行時間
実行使用メモリ
testcase_00 AC 37 ms
53,208 KB
testcase_01 AC 56 ms
61,492 KB
testcase_02 AC 37 ms
53,668 KB
testcase_03 AC 39 ms
53,136 KB
testcase_04 AC 109 ms
62,404 KB
testcase_05 AC 127 ms
63,212 KB
testcase_06 AC 121 ms
63,356 KB
testcase_07 AC 121 ms
63,652 KB
testcase_08 AC 54 ms
58,616 KB
testcase_09 RE -
testcase_10 RE -
testcase_11 RE -
testcase_12 RE -
testcase_13 AC 41 ms
52,456 KB
testcase_14 RE -
testcase_15 MLE -
testcase_16 -- -
testcase_17 -- -
testcase_18 -- -
testcase_19 -- -
testcase_20 -- -
testcase_21 -- -
testcase_22 -- -
testcase_23 -- -
testcase_24 -- -
権限があれば一括ダウンロードができます

ソースコード

diff # 

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)
0