結果
問題 | No.2829 GCD Divination |
ユーザー | 👑 emthrm |
提出日時 | 2024-08-03 16:47:50 |
言語 | C++23 (gcc 12.3.0 + boost 1.83.0) |
結果 |
AC
|
実行時間 | 131 ms / 2,000 ms |
コード長 | 2,489 bytes |
コンパイル時間 | 3,187 ms |
コンパイル使用メモリ | 252,392 KB |
実行使用メモリ | 120,716 KB |
最終ジャッジ日時 | 2024-08-03 16:47:57 |
合計ジャッジ時間 | 6,293 ms |
ジャッジサーバーID (参考情報) |
judge5 / judge3 |
(要ログイン)
テストケース
テストケース表示入力 | 結果 | 実行時間 実行使用メモリ |
---|---|---|
testcase_00 | AC | 2 ms
6,812 KB |
testcase_01 | AC | 2 ms
6,944 KB |
testcase_02 | AC | 131 ms
120,716 KB |
testcase_03 | AC | 2 ms
6,940 KB |
testcase_04 | AC | 131 ms
120,648 KB |
testcase_05 | AC | 118 ms
104,936 KB |
testcase_06 | AC | 96 ms
88,156 KB |
testcase_07 | AC | 83 ms
79,752 KB |
testcase_08 | AC | 58 ms
54,252 KB |
testcase_09 | AC | 47 ms
45,776 KB |
testcase_10 | AC | 71 ms
69,068 KB |
testcase_11 | AC | 6 ms
6,944 KB |
testcase_12 | AC | 23 ms
24,676 KB |
testcase_13 | AC | 116 ms
108,232 KB |
testcase_14 | AC | 74 ms
73,972 KB |
testcase_15 | AC | 52 ms
54,272 KB |
testcase_16 | AC | 25 ms
27,256 KB |
testcase_17 | AC | 15 ms
16,256 KB |
testcase_18 | AC | 12 ms
12,736 KB |
testcase_19 | AC | 18 ms
20,136 KB |
testcase_20 | AC | 2 ms
6,944 KB |
testcase_21 | AC | 45 ms
46,264 KB |
testcase_22 | AC | 54 ms
54,696 KB |
testcase_23 | AC | 67 ms
64,464 KB |
testcase_24 | AC | 103 ms
95,780 KB |
testcase_25 | AC | 125 ms
117,140 KB |
testcase_26 | AC | 48 ms
51,228 KB |
testcase_27 | AC | 89 ms
87,632 KB |
testcase_28 | AC | 18 ms
20,964 KB |
testcase_29 | AC | 56 ms
55,052 KB |
testcase_30 | AC | 122 ms
112,744 KB |
testcase_31 | AC | 16 ms
16,456 KB |
testcase_32 | AC | 69 ms
67,508 KB |
testcase_33 | AC | 50 ms
49,284 KB |
testcase_34 | AC | 91 ms
87,456 KB |
ソースコード
#include <bits/stdc++.h> using namespace std; #define FOR(i,m,n) for(int i=(m);i<(n);++i) #define REP(i,n) FOR(i,0,n) using ll = long long; constexpr int INF = 0x3f3f3f3f; constexpr long long LINF = 0x3f3f3f3f3f3f3f3fLL; constexpr double EPS = 1e-8; constexpr int MOD = 998244353; // constexpr int MOD = 1000000007; constexpr int DY4[]{1, 0, -1, 0}, DX4[]{0, -1, 0, 1}; constexpr int DY8[]{1, 1, 0, -1, -1, -1, 0, 1}; constexpr int DX8[]{0, -1, -1, -1, 0, 1, 1, 1}; template <typename T, typename U> inline bool chmax(T& a, U b) { return a < b ? (a = b, true) : false; } template <typename T, typename U> inline bool chmin(T& a, U b) { return a > b ? (a = b, true) : false; } struct IOSetup { IOSetup() { std::cin.tie(nullptr); std::ios_base::sync_with_stdio(false); std::cout << fixed << setprecision(20); } } iosetup; template <bool GETS_ONLY_PRIME> std::vector<int> prime_sieve(const int n) { std::vector<int> smallest_prime_factor(n + 1), prime; std::iota(smallest_prime_factor.begin(), smallest_prime_factor.end(), 0); for (int i = 2; i <= n; ++i) { if (smallest_prime_factor[i] == i) [[unlikely]] prime.emplace_back(i); for (const int p : prime) { if (i * p > n || p > smallest_prime_factor[i]) break; smallest_prime_factor[i * p] = p; } } return GETS_ONLY_PRIME ? prime : smallest_prime_factor; } struct Divisor { const std::vector<int> smallest_prime_factor; explicit Divisor(const int n) : smallest_prime_factor(prime_sieve<false>(n)) {} std::vector<int> query(int n) const { std::vector<int> res{1}; while (n > 1) { const int prime_factor = smallest_prime_factor[n], d = res.size(); int tmp = 1; for (; n % prime_factor == 0; n /= prime_factor) { tmp *= prime_factor; for (int i = 0; i < d; ++i) { res.emplace_back(res[i] * tmp); } } } // std::sort(res.begin(), res.end()); return res; } }; int main() { int n; cin >> n; const Divisor divisor(n); vector<double> dp(n + 1, 0); for (const int d : divisor.query(n)) { if (d == 1) continue; const vector<int> ds = divisor.query(d); const int l = ds.size(); vector<int> num(l, 0); for (int i = l - 1; i >= 0; --i) { num[i] = d / ds[i]; FOR(j, i + 1, l) { if (ds[j] % ds[i] == 0) num[i] -= num[j]; } } dp[d] = d; FOR(i, 1, l - 1) dp[d] += num[i] * dp[ds[i]]; dp[d] /= d - 1; } cout << dp[n] << '\n'; return 0; }