結果

問題 No.2829 GCD Divination
ユーザー 👑 emthrmemthrm
提出日時 2024-08-03 16:36:19
言語 C++23
(gcc 12.3.0 + boost 1.83.0)
結果
TLE  
実行時間 -
コード長 1,408 bytes
コンパイル時間 3,359 ms
コンパイル使用メモリ 252,308 KB
実行使用メモリ 81,664 KB
最終ジャッジ日時 2024-08-03 16:36:52
合計ジャッジ時間 31,320 ms
ジャッジサーバーID
（参考情報） 		judge4 / judge5
このコードへのチャレンジ
（要ログイン）

テストケース

テストケース表示
入力 結果 実行時間
実行使用メモリ
testcase_00 AC 2 ms
5,248 KB
testcase_01 AC 2 ms
5,376 KB
testcase_02 AC 989 ms
81,664 KB
testcase_03 AC 2 ms
5,376 KB
testcase_04 AC 1,895 ms
81,432 KB
testcase_05 TLE -
testcase_06 TLE -
testcase_07 TLE -
testcase_08 AC 1,294 ms
37,376 KB
testcase_09 AC 1,067 ms
31,480 KB
testcase_10 AC 772 ms
47,064 KB
testcase_11 AC 20 ms
5,376 KB
testcase_12 AC 233 ms
17,360 KB
testcase_13 AC 955 ms
73,216 KB
testcase_14 AC 682 ms
50,304 KB
testcase_15 AC 787 ms
37,248 KB
testcase_16 AC 257 ms
19,328 KB
testcase_17 AC 137 ms
12,160 KB
testcase_18 AC 130 ms
9,608 KB
testcase_19 AC 180 ms
14,592 KB
testcase_20 AC 2 ms
5,376 KB
testcase_21 AC 751 ms
31,788 KB
testcase_22 AC 397 ms
37,416 KB
testcase_23 AC 749 ms
44,032 KB
testcase_24 AC 899 ms
65,152 KB
testcase_25 AC 1,006 ms
79,232 KB
testcase_26 AC 440 ms
35,240 KB
testcase_27 AC 1,242 ms
59,520 KB
testcase_28 AC 155 ms
15,104 KB
testcase_29 AC 758 ms
37,760 KB
testcase_30 AC 1,283 ms
76,132 KB
testcase_31 AC 105 ms
12,160 KB
testcase_32 AC 562 ms
45,824 KB
testcase_33 AC 447 ms
33,920 KB
testcase_34 AC 1,173 ms
59,264 KB
権限があれば一括ダウンロードができます

ソースコード

diff # 

#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 <typename T>
std::vector<T> divisor(const T n) {
  std::vector<T> res;
  T i = 1;
  for (; i * i < n; ++i) {
    if (n % i == 0) [[unlikely]] {
      res.emplace_back(i);
      res.emplace_back(n / i);
    }
  }
  if (i * i == n && n % i == 0) res.emplace_back(i);
  std::sort(res.begin(), res.end());
  return res;
}

int main() {
  int n; cin >> n;
  vector<double> dp(n + 1, 0);
  for (const int d : divisor(n)) {
    if (d == 1) continue;
    dp[d] = d;
    FOR(m, 1, d) dp[d] += dp[gcd(d, m)];
    dp[d] /= d - 1;
  }
  cout << dp[n] << '\n';
  return 0;
}
0