結果
問題 | No.1791 Repeat Multiplication |
ユーザー | commy |
提出日時 | 2022-09-11 04:10:54 |
言語 | C++17 (gcc 12.3.0 + boost 1.83.0) |
結果 |
AC
|
実行時間 | 936 ms / 3,000 ms |
コード長 | 3,386 bytes |
コンパイル時間 | 940 ms |
コンパイル使用メモリ | 88,764 KB |
実行使用メモリ | 32,464 KB |
最終ジャッジ日時 | 2023-08-18 08:02:38 |
合計ジャッジ時間 | 21,508 ms |
ジャッジサーバーID (参考情報) |
judge13 / judge14 |
(要ログイン)
テストケース
テストケース表示入力 | 結果 | 実行時間 実行使用メモリ |
---|---|---|
testcase_00 | AC | 2 ms
4,380 KB |
testcase_01 | AC | 2 ms
4,380 KB |
testcase_02 | AC | 2 ms
4,380 KB |
testcase_03 | AC | 2 ms
4,380 KB |
testcase_04 | AC | 1 ms
4,376 KB |
testcase_05 | AC | 2 ms
4,380 KB |
testcase_06 | AC | 2 ms
4,376 KB |
testcase_07 | AC | 1 ms
4,376 KB |
testcase_08 | AC | 62 ms
4,836 KB |
testcase_09 | AC | 56 ms
4,608 KB |
testcase_10 | AC | 53 ms
4,644 KB |
testcase_11 | AC | 50 ms
4,640 KB |
testcase_12 | AC | 29 ms
4,376 KB |
testcase_13 | AC | 832 ms
30,008 KB |
testcase_14 | AC | 840 ms
30,732 KB |
testcase_15 | AC | 613 ms
23,324 KB |
testcase_16 | AC | 592 ms
22,620 KB |
testcase_17 | AC | 896 ms
31,524 KB |
testcase_18 | AC | 701 ms
25,968 KB |
testcase_19 | AC | 816 ms
29,864 KB |
testcase_20 | AC | 641 ms
24,188 KB |
testcase_21 | AC | 810 ms
29,284 KB |
testcase_22 | AC | 870 ms
31,808 KB |
testcase_23 | AC | 643 ms
24,608 KB |
testcase_24 | AC | 844 ms
30,584 KB |
testcase_25 | AC | 764 ms
28,248 KB |
testcase_26 | AC | 857 ms
31,208 KB |
testcase_27 | AC | 691 ms
25,976 KB |
testcase_28 | AC | 889 ms
32,336 KB |
testcase_29 | AC | 905 ms
32,464 KB |
testcase_30 | AC | 896 ms
32,444 KB |
testcase_31 | AC | 2 ms
4,380 KB |
testcase_32 | AC | 904 ms
32,268 KB |
testcase_33 | AC | 903 ms
32,448 KB |
testcase_34 | AC | 936 ms
32,380 KB |
ソースコード
#include <algorithm> #include <iostream> #include <numeric> #include <string> #include <tuple> #include <utility> #include <vector> #include <limits> #define rep(i, a, b) for (int i = int(a); i < int(b); i++) using namespace std; using ll = long long int; // NOLINT using P = pair<ll, ll>; // clang-format off #ifdef _DEBUG_ #define dump(...) do{ cerr << __LINE__ << ":\t" << #__VA_ARGS__ << " = "; debug_print(__VA_ARGS__); } while(false) template<typename T, typename... Ts> void debug_print(const T &t, const Ts &...ts) { cerr << t; ((cerr << ", " << ts), ...); cerr << endl; } #else #define dump(...) do{ } while(false) #endif template<typename T> vector<T> make_v(size_t a, T b) { return vector<T>(a, b); } template<typename... Ts> auto make_v(size_t a, Ts... ts) { return vector<decltype(make_v(ts...))>(a, make_v(ts...)); } template<typename T> bool chmin(T &a, const T& b) { if (a > b) {a = b; return true; } return false; } template<typename T> bool chmax(T &a, const T& b) { if (a < b) {a = b; return true; } return false; } template<typename T, typename... Ts> void print(const T& t, const Ts&... ts) { cout << t; ((cout << ' ' << ts), ...); cout << '\n'; } template<typename... Ts> void input(Ts&... ts) { (cin >> ... >> ts); } template<typename T> istream &operator,(istream &in, T &t) { return in >> t; } struct Inf { template<typename T> constexpr operator T() { return numeric_limits<T>::max() / 2; } }; // clang-format on class Eratosthenes { vector<int> plist; public: Eratosthenes(int sz) : plist(sz, -1) { for (int i = 2; i < sz; i++) { if (plist[i] == -1) { for (int j = i; j < sz; j += i) { if (plist[j] == -1) { plist[j] = i; } } } } } bool is_prime(int n) { return plist[n] == n; } vector<int> factor(int n) { vector<int> res; while (n > 1) { res.push_back(plist[n]); n /= plist[n]; } return res; } vector<int> divisor(int n) { vector<pair<int, int>> runlength; for (auto e : this->factor(n)) { if (runlength.size() == 0 || runlength.back().first != e) { runlength.emplace_back(e, 1); } else { runlength.back().second++; } } vector<int> res; auto dfs = [&](auto &f, unsigned int idx, int d) { if (idx == runlength.size()) { res.push_back(d); return; } for (int i = 0, m = 1; i <= runlength[idx].second; i++, m *= runlength[idx].first) { f(f, idx + 1, d * m); } }; dfs(dfs, 0, 1); // sort(res.begin(), res.end()); return res; } }; int main() { cin.tie(nullptr); ios::sync_with_stdio(false); int n, q; cin, n, q; vector<ll> ans1(n + 1, 0), ans2(n + 1, 1); ans1[1] = 1; Eratosthenes erat(n + 1); for (int i = 1; i <= n; i++) { for (int j = 2 * i; j <= n; j += i) { ans1[j] += ans1[i]; } } for (int i = n; i >= 1; i--) { for (auto e : erat.divisor(i)) { ans2[e] += ans2[i]; } } rep(i, 0, q) { int x; cin, x; print(ans1[x] * ans2[x] / 2); } return 0; }