結果
問題 | No.1973 Divisor Sequence |
ユーザー |
![]() |
提出日時 | 2022-06-10 22:09:12 |
言語 | PyPy3 (7.3.15) |
結果 |
TLE
|
実行時間 | - |
コード長 | 652 bytes |
コンパイル時間 | 266 ms |
コンパイル使用メモリ | 82,292 KB |
実行使用メモリ | 76,800 KB |
最終ジャッジ日時 | 2024-09-21 06:23:09 |
合計ジャッジ時間 | 9,562 ms |
ジャッジサーバーID (参考情報) |
judge2 / judge5 |
(要ログイン)
ファイルパターン | 結果 |
---|---|
sample | AC * 3 |
other | AC * 20 TLE * 1 -- * 1 |
ソースコード
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) for i in range(k+1): for j in range(k+1): if i+j<=k: dp_n[j]+=dp[i] 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)