結果
問題 | No.1973 Divisor Sequence |
ユーザー | magsta |
提出日時 | 2022-06-06 22:09:37 |
言語 | C++14 (gcc 12.3.0 + boost 1.83.0) |
結果 |
AC
|
実行時間 | 1,770 ms / 2,000 ms |
コード長 | 1,202 bytes |
コンパイル時間 | 1,760 ms |
コンパイル使用メモリ | 175,704 KB |
実行使用メモリ | 179,100 KB |
最終ジャッジ日時 | 2023-10-21 06:17:00 |
合計ジャッジ時間 | 6,478 ms |
ジャッジサーバーID (参考情報) |
judge10 / judge14 |
(要ログイン)
テストケース
テストケース表示入力 | 結果 | 実行時間 実行使用メモリ |
---|---|---|
testcase_00 | AC | 2 ms
4,348 KB |
testcase_01 | AC | 2 ms
4,348 KB |
testcase_02 | AC | 173 ms
30,064 KB |
testcase_03 | AC | 11 ms
5,540 KB |
testcase_04 | AC | 102 ms
18,988 KB |
testcase_05 | AC | 32 ms
16,680 KB |
testcase_06 | AC | 110 ms
25,520 KB |
testcase_07 | AC | 22 ms
7,048 KB |
testcase_08 | AC | 59 ms
14,900 KB |
testcase_09 | AC | 90 ms
16,944 KB |
testcase_10 | AC | 115 ms
20,376 KB |
testcase_11 | AC | 16 ms
5,492 KB |
testcase_12 | AC | 89 ms
20,940 KB |
testcase_13 | AC | 28 ms
8,436 KB |
testcase_14 | AC | 74 ms
17,872 KB |
testcase_15 | AC | 155 ms
22,012 KB |
testcase_16 | AC | 149 ms
26,048 KB |
testcase_17 | AC | 96 ms
17,472 KB |
testcase_18 | AC | 7 ms
4,472 KB |
testcase_19 | AC | 128 ms
22,072 KB |
testcase_20 | AC | 22 ms
6,256 KB |
testcase_21 | AC | 134 ms
30,536 KB |
testcase_22 | AC | 124 ms
22,156 KB |
testcase_23 | AC | 1,770 ms
179,100 KB |
testcase_24 | AC | 344 ms
46,308 KB |
ソースコード
#include <bits/stdc++.h> using namespace std; #define MOD 1000000007 vector<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 <= factor[q].second - 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; }