結果
問題 | No.1973 Divisor Sequence |
ユーザー |
![]() |
提出日時 | 2022-06-06 22:07:18 |
言語 | C++14 (gcc 13.3.0 + boost 1.87.0) |
結果 |
WA
|
実行時間 | - |
コード長 | 1,183 bytes |
コンパイル時間 | 1,876 ms |
コンパイル使用メモリ | 176,256 KB |
実行使用メモリ | 178,940 KB |
最終ジャッジ日時 | 2024-09-21 05:50:16 |
合計ジャッジ時間 | 6,385 ms |
ジャッジサーバーID (参考情報) |
judge4 / judge2 |
(要ログイン)
ファイルパターン | 結果 |
---|---|
sample | AC * 2 WA * 1 |
other | WA * 22 |
ソースコード
#include <bits/stdc++.h>using namespace std;#define MOD 1000000007vector<pair<long long, long long>> fact(long long N) {vector<pair<long long, long long> > res;for (long long a = 2; a * a <= N; ++a) {if (N % a != 0) continue;long long ex = 0;while (N % a == 0) {++ex;N /= a;}res.push_back({ a, ex });}if (N != 1) res.push_back({ N, 1 });return res;}int main() {long long N, M; cin >> N >> M;auto factor = fact(M);long long N_ = factor.size();long long ans = 1;for (int q = 0; q < N_; q++) {vector<vector<long long>> dp(N + 1, vector<long long>(factor[q].second + 1, 0));dp[0][0] = 1;for (int i = 1; i <= N; i++) {for (int j = factor[q].second; j >= 0; j--) {for (int k = 0; k <= j; k++) {dp[i][j] = (dp[i][j] + dp[i - 1][k]) % MOD;}}}long long sum = 0;for (int i = 0; i <= factor[q].second; i++) {sum = (sum + dp[N][i]) % MOD;}ans = ans * sum % MOD;}cout << ans << endl;}