結果
問題 | No.1973 Divisor Sequence |
ユーザー |
![]() |
提出日時 | 2022-06-10 22:10:54 |
言語 | PyPy3 (7.3.15) |
結果 |
AC
|
実行時間 | 671 ms / 2,000 ms |
コード長 | 628 bytes |
コンパイル時間 | 214 ms |
コンパイル使用メモリ | 82,352 KB |
実行使用メモリ | 76,468 KB |
最終ジャッジ日時 | 2024-09-21 07:30:17 |
合計ジャッジ時間 | 6,344 ms |
ジャッジサーバーID (参考情報) |
judge1 / judge2 |
(要ログイン)
ファイルパターン | 結果 |
---|---|
sample | AC * 3 |
other | AC * 22 |
ソースコード
def fact(n):arr = []temp = nfor i in range(2, int(-(-n**0.5//1))+1):if temp%i==0:cnt=0while temp%i==0:cnt+=1temp //= iarr.append([i, cnt])if temp!=1:arr.append([temp, 1])if arr==[]:arr.append([n, 1])return arrdef calc(n,k):dp=[0]*(k+1)dp[0]=1for _ in range(n):dp_n=[0]*(k+1)cum=0for i in range(k+1):cum+=dp[i]dp_n[k-i]=cumdp=[i%mod for i in dp_n]return sum(dp)%modmod=10**9+7n,m=map(int,input().split())ans=1if m==1:print(1)exit()f=fact(m)for _,k in f:ans*=calc(n,k)ans%=modprint(ans)