結果

問題 No.2829 GCD Divination
ユーザー HaruiHarui
提出日時 2024-08-02 22:04:48
言語 C++23
(gcc 12.3.0 + boost 1.83.0)
結果
WA  
実行時間 -
コード長 1,068 bytes
コンパイル時間 1,268 ms
コンパイル使用メモリ 99,216 KB
実行使用メモリ 6,948 KB
最終ジャッジ日時 2024-08-02 22:04:51
合計ジャッジ時間 2,437 ms
ジャッジサーバーID
(参考情報)
judge1 / judge5
このコードへのチャレンジ
(要ログイン)

テストケース

テストケース表示
入力 結果 実行時間
実行使用メモリ
testcase_00 AC 1 ms
6,812 KB
testcase_01 AC 1 ms
6,940 KB
testcase_02 AC 1 ms
6,940 KB
testcase_03 AC 1 ms
6,940 KB
testcase_04 WA -
testcase_05 WA -
testcase_06 WA -
testcase_07 WA -
testcase_08 WA -
testcase_09 WA -
testcase_10 WA -
testcase_11 WA -
testcase_12 WA -
testcase_13 WA -
testcase_14 WA -
testcase_15 WA -
testcase_16 WA -
testcase_17 WA -
testcase_18 WA -
testcase_19 WA -
testcase_20 WA -
testcase_21 WA -
testcase_22 AC 2 ms
6,944 KB
testcase_23 WA -
testcase_24 WA -
testcase_25 WA -
testcase_26 WA -
testcase_27 WA -
testcase_28 WA -
testcase_29 WA -
testcase_30 WA -
testcase_31 WA -
testcase_32 WA -
testcase_33 WA -
testcase_34 WA -
権限があれば一括ダウンロードができます

ソースコード

diff #

#include <iostream>
#include <iomanip>
#include <cassert>
#include <map>
using namespace std;
using ll = long long;
template<class T> inline bool chmax(T& a, const T& b) {if (a<b) {a=b; return true;} return false;}
template<class T> inline bool chmin(T& a, const T& b) {if (b<a) {a=b; return true;} return false;}
const int INTINF = 1000001000;
const int INTMAX = 2147483647;
const ll LLMAX = 9223372036854775807;
const ll LLINF = 1000000000000000000;

int main() {
  ios::sync_with_stdio(0); cin.tie(0); cout.tie(0);
  ll N; cin >> N;
  map<ll,double> memo;

  memo[1] = 0.0;

  auto solve = [&memo] (auto self, ll n) -> double {
    if (memo.contains(n)) return memo.at(n);

    double divsum = 0;
    for (int d=2; d*d <= n; d++) {
      if (n % d == 0) {
        divsum += self(self, d);

        if (d != n/d) {
          divsum += self(self, n/d);
        }
      }
    }

    double ret = double(n) / double(n-1) + divsum / double(n-1);
    memo[n] = ret;
    return ret;
  };
  cout << fixed << setprecision(10);
  cout << solve(solve,N) << endl;
  return 0;
}
0