結果

問題 No.2645 Sum of Divisors?
ユーザー KKT89KKT89
提出日時 2024-11-05 01:48:16
言語 C++17
(gcc 12.3.0 + boost 1.83.0)
結果
AC  
実行時間 268 ms / 2,000 ms
コード長 1,635 bytes
コンパイル時間 3,064 ms
コンパイル使用メモリ 219,908 KB
実行使用メモリ 54,132 KB
最終ジャッジ日時 2024-11-05 01:48:23
合計ジャッジ時間 6,800 ms
ジャッジサーバーID
（参考情報） 		judge4 / judge1
このコードへのチャレンジ
（要ログイン）

テストケース

テストケース表示
入力 結果 実行時間
実行使用メモリ
testcase_00 AC 2 ms
6,816 KB
testcase_01 AC 2 ms
6,816 KB
testcase_02 AC 161 ms
52,956 KB
testcase_03 AC 2 ms
6,820 KB
testcase_04 AC 2 ms
6,816 KB
testcase_05 AC 2 ms
6,820 KB
testcase_06 AC 266 ms
52,460 KB
testcase_07 AC 268 ms
53,844 KB
testcase_08 AC 2 ms
6,816 KB
testcase_09 AC 2 ms
6,820 KB
testcase_10 AC 2 ms
6,820 KB
testcase_11 AC 2 ms
6,820 KB
testcase_12 AC 2 ms
6,816 KB
testcase_13 AC 2 ms
6,824 KB
testcase_14 AC 2 ms
6,816 KB
testcase_15 AC 2 ms
6,816 KB
testcase_16 AC 2 ms
6,816 KB
testcase_17 AC 2 ms
6,816 KB
testcase_18 AC 3 ms
6,820 KB
testcase_19 AC 2 ms
6,816 KB
testcase_20 AC 5 ms
6,820 KB
testcase_21 AC 4 ms
6,816 KB
testcase_22 AC 6 ms
6,820 KB
testcase_23 AC 11 ms
6,816 KB
testcase_24 AC 30 ms
10,672 KB
testcase_25 AC 23 ms
11,176 KB
testcase_26 AC 37 ms
10,088 KB
testcase_27 AC 49 ms
17,728 KB
testcase_28 AC 99 ms
28,800 KB
testcase_29 AC 261 ms
53,552 KB
testcase_30 AC 109 ms
29,004 KB
testcase_31 AC 165 ms
54,132 KB
testcase_32 AC 213 ms
53,384 KB
testcase_33 AC 113 ms
28,836 KB
testcase_34 AC 200 ms
52,816 KB
権限があれば一括ダウンロードができます

ソースコード

diff # 

#pragma GCC optimize("Ofast")
#include <bits/stdc++.h>
using namespace std;
typedef long long int ll;
typedef unsigned long long int ull;

mt19937_64 rng(chrono::steady_clock::now().time_since_epoch().count());
ll myRand(ll B) { return (ull)rng() % B; }
inline double time() {
    return static_cast<long double>(chrono::duration_cast<chrono::nanoseconds>(chrono::steady_clock::now().time_since_epoch()).count()) * 1e-9;
}

typedef long double ld;
ld Harmonic(ll n) {
    static vector<ld> dp;
    if (dp.empty()) {
        dp.resize(1000);
        for (int i = 1; i < 1000; ++i) {
            dp[i] = dp[i - 1] + ld(1) / i;
        }
    }
    if (n < 1000) return dp[n];
    return logl(n) + 0.57721566490153286060651209 + ld(1) / (2LL * n) - ld(1) / 12 / n / n + ld(1) / 120 / n / n / n / n;
}

// {x,{l,r}} : n/i = x (l <= i <= r) を{l,r}の昇順に返す
template <typename T> vector<pair<T, pair<T, T>>> quotient_range(T n) {
    vector<pair<T, pair<T, T>>> res;
    T m = 1;
    while (m * m <= n) {
        res.push_back({n / m, {m, m}});
        m += 1;
    }
    for (int i = m; i >= 1; --i) {
        T l = n / (i + 1) + 1;
        T r = n / i;
        if (l <= r and res.back().second.second < l) res.push_back({n / l, {l, r}});
    }
    return res;
}

int main() {
    cin.tie(nullptr);
    ios::sync_with_stdio(false);
    ll n;
    cin >> n;
    auto q = quotient_range(n);
    ld res = 0;
    for (auto [x, lr] : q) {
        ld h1 = Harmonic(lr.second);
        ld h2 = Harmonic(lr.first - 1);
        ld hx = Harmonic(x);
        res += (h1 - h2) * hx;
    }
    cout << fixed << setprecision(10) << res << '\n';
}
0