結果
問題 | 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 = 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 def calc(n,k): dp=[0]*(k+1) dp[0]=1 for _ in range(n): dp_n=[0]*(k+1) cum=0 for i in range(k+1): cum+=dp[i] dp_n[k-i]=cum dp=[i%mod for i in dp_n] return sum(dp)%mod mod=10**9+7 n,m=map(int,input().split()) ans=1 if m==1: print(1) exit() f=fact(m) for _,k in f: ans*=calc(n,k) ans%=mod print(ans)