結果
問題 | No.1973 Divisor Sequence |
ユーザー | eve__fuyuki |
提出日時 | 2024-06-17 16:48:35 |
言語 | C++17(gcc12) (gcc 12.3.0 + boost 1.87.0) |
結果 |
AC
|
実行時間 | 541 ms / 2,000 ms |
コード長 | 1,067 bytes |
コンパイル時間 | 2,732 ms |
コンパイル使用メモリ | 207,624 KB |
実行使用メモリ | 198,496 KB |
最終ジャッジ日時 | 2024-06-17 16:48:43 |
合計ジャッジ時間 | 7,832 ms |
ジャッジサーバーID (参考情報) |
judge1 / judge2 |
(要ログイン)
ファイルパターン | 結果 |
---|---|
sample | AC * 3 |
other | AC * 22 |
ソースコード
#include <bits/stdc++.h> #include <atcoder/modint> using namespace std; using mint = atcoder::modint1000000007; void fast_io() { ios_base::sync_with_stdio(false); cin.tie(nullptr); } mint calc(int n, int cnt) { vector<vector<mint>> dp(n, vector<mint>(cnt + 1)), cum(n, vector<mint>(cnt + 2, 0)); for (int i = 0; i <= cnt; i++) { dp[0][i] = 1; cum[0][i + 1] = cum[0][i] + dp[0][i]; } for (int i = 1; i < n; i++) { for (int j = 0; j <= cnt; j++) { dp[i][j] = cum[i - 1][cnt - j + 1]; cum[i][j + 1] = cum[i][j] + dp[i][j]; } } return cum[n - 1][cnt + 1]; } int main() { fast_io(); int n; long long m; cin >> n >> m; mint ans = 1; for (long long i = 2; i * i <= m; i++) { if (m % i == 0) { int cnt = 0; while (m % i == 0) { m /= i; cnt++; } ans *= calc(n, cnt); } } if (m > 1) { ans *= calc(n, 1); } cout << ans.val() << endl; }