結果
問題 | No.1973 Divisor Sequence |
ユーザー | magsta |
提出日時 | 2021-08-20 12:21:47 |
言語 | C++14 (gcc 13.2.0 + boost 1.83.0) |
結果 |
AC
|
実行時間 | 309 ms / 2,000 ms |
コード長 | 1,186 bytes |
コンパイル時間 | 1,722 ms |
コンパイル使用メモリ | 175,612 KB |
実行使用メモリ | 179,100 KB |
最終ジャッジ日時 | 2023-10-21 06:16:48 |
合計ジャッジ時間 | 4,947 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 | 101 ms
18,988 KB |
testcase_05 | AC | 32 ms
16,680 KB |
testcase_06 | AC | 106 ms
25,520 KB |
testcase_07 | AC | 22 ms
7,048 KB |
testcase_08 | AC | 57 ms
14,900 KB |
testcase_09 | AC | 88 ms
16,944 KB |
testcase_10 | AC | 112 ms
20,376 KB |
testcase_11 | AC | 15 ms
5,492 KB |
testcase_12 | AC | 86 ms
20,940 KB |
testcase_13 | AC | 27 ms
8,436 KB |
testcase_14 | AC | 60 ms
17,872 KB |
testcase_15 | AC | 151 ms
22,012 KB |
testcase_16 | AC | 147 ms
26,048 KB |
testcase_17 | AC | 93 ms
17,472 KB |
testcase_18 | AC | 6 ms
4,472 KB |
testcase_19 | AC | 122 ms
22,072 KB |
testcase_20 | AC | 22 ms
6,256 KB |
testcase_21 | AC | 132 ms
30,536 KB |
testcase_22 | AC | 122 ms
22,156 KB |
testcase_23 | AC | 256 ms
179,100 KB |
testcase_24 | AC | 309 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++) { long long rui = 0; for (int j = factor[q].second; j >= 0; j--) { rui = (rui + dp[i - 1][factor[q].second - j]) % MOD; dp[i][j] = rui; } } 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; }