結果
問題 | No.1973 Divisor Sequence |
ユーザー |
![]() |
提出日時 | 2021-08-21 12:24:06 |
言語 | C++17 (gcc 13.3.0 + boost 1.87.0) |
結果 |
AC
|
実行時間 | 296 ms / 2,000 ms |
コード長 | 957 bytes |
コンパイル時間 | 2,634 ms |
コンパイル使用メモリ | 199,672 KB |
最終ジャッジ日時 | 2025-01-24 01:20:44 |
ジャッジサーバーID (参考情報) |
judge5 / judge2 |
(要ログイン)
ファイルパターン | 結果 |
---|---|
sample | AC * 3 |
other | AC * 22 |
ソースコード
#include <bits/stdc++.h> using namespace std; #include <atcoder/modint> using mint = atcoder::modint1000000007; #define rep(i, x) for(int i = 0; i < (x); i++) using ll = long long; mint res = 1; int main() { int n; ll m; cin >> n >> m; int lim = sqrt(m); for(int i = 2; i <= lim; i++) { int c = 0; while(m % i == 0) m /= i, c++; if(c) { vector dp(n + 1, vector(c + 1, mint(0))); dp[0][0] = 1; for(int i = 0; i < n; i++) { mint cum = 0; for(int j = c; j >= 0; j--) { cum += dp[i][c - j]; dp[i + 1][j] = cum; } } res *= accumulate(dp.back().begin(), dp.back().end(), mint(0)); } } if(m > 1) { int c = 1; vector dp(n + 1, vector(c + 1, mint(0))); dp[0][0] = 1; for(int i = 0; i < n; i++) { mint cum = 0; for(int j = c; j >= 0; j--) { cum += dp[i][c - j]; dp[i + 1][j] = cum; } } res *= accumulate(dp.back().begin(), dp.back().end(), mint(0)); } cout << res.val() << endl; }