結果

問題 No.1659 Product of Divisors
ユーザー titiatitia
提出日時 2021-08-28 01:44:36
言語 Python3
(3.12.2 + numpy 1.26.4 + scipy 1.12.0)
結果
TLE  
実行時間 -
コード長 616 bytes
コンパイル時間 111 ms
コンパイル使用メモリ 12,416 KB
実行使用メモリ 21,504 KB
最終ジャッジ日時 2024-11-21 09:08:20
合計ジャッジ時間 51,644 ms
ジャッジサーバーID
（参考情報） 		judge2 / judge3
このコードへのチャレンジ
（要ログイン）

テストケース

テストケース表示
入力 結果 実行時間
実行使用メモリ
testcase_00 AC 30 ms
16,000 KB
testcase_01 AC 84 ms
17,568 KB
testcase_02 AC 35 ms
16,128 KB
testcase_03 AC 33 ms
21,504 KB
testcase_04 TLE -
testcase_05 AC 1,351 ms
21,248 KB
testcase_06 AC 619 ms
16,000 KB
testcase_07 AC 150 ms
21,248 KB
testcase_08 AC 163 ms
17,696 KB
testcase_09 TLE -
testcase_10 TLE -
testcase_11 TLE -
testcase_12 TLE -
testcase_13 AC 31 ms
16,000 KB
testcase_14 TLE -
testcase_15 TLE -
testcase_16 TLE -
testcase_17 TLE -
testcase_18 TLE -
testcase_19 TLE -
testcase_20 TLE -
testcase_21 TLE -
testcase_22 TLE -
testcase_23 TLE -
testcase_24 TLE -
権限があれば一括ダウンロードができます

ソースコード

diff # 

import sys
input = sys.stdin.readline
import math 

N,K=map(int,input().split())
mod=10**9+7


x=N
L=int(math.sqrt(x))

FACT=dict()

for i in range(2,L+2):
    while x%i==0:
        FACT[i]=FACT.get(i,0)+1
        x=x//i

if x!=1:
    FACT[x]=FACT.get(x,0)+1

def calc(x):
    DP=[0]*(x+1)
    DP[-1]=1

    for i in range(K):
        NDP=[0]*(x+1)

        for j in range(x+1):
            if DP[j]!=0:
                for k in range(j+1):
                    NDP[k]+=DP[j]
                    NDP[k]%=mod
        DP=NDP

    return sum(DP)

ANS=1
for f in FACT:
    ANS*=calc(FACT[f])
    ANS%=mod

print(ANS%mod)
0