結果
問題 | No.2645 Sum of Divisors? |
ユーザー |
![]() |
提出日時 | 2024-02-19 22:42:09 |
言語 | C++23 (gcc 13.3.0 + boost 1.87.0) |
結果 |
AC
|
実行時間 | 126 ms / 2,000 ms |
コード長 | 1,156 bytes |
コンパイル時間 | 4,893 ms |
コンパイル使用メモリ | 309,280 KB |
実行使用メモリ | 19,428 KB |
最終ジャッジ日時 | 2024-09-29 02:33:22 |
合計ジャッジ時間 | 6,102 ms |
ジャッジサーバーID (参考情報) |
judge3 / judge4 |
(要ログイン)
ファイルパターン | 結果 |
---|---|
sample | AC * 3 |
other | AC * 32 |
ソースコード
#include <atcoder/all>#include <bits/stdc++.h>using ll = long long;using ull = unsigned long long;#define rep(i, n) for(int i = 0; i < (int)(n); i++)#define REP(i, m, n) for(int i = (int)(m); i < (int)(n); i++)using namespace std;using namespace atcoder;using mint = modint998244353;const int inf = 1000000007;const ll longinf = 1ll << 60;int main() {ios::sync_with_stdio(false);cin.tie(0);cout.tie(0);ll n;cin >> n;vector<long double> sum(1000001);REP(i, 1, 1000001) {sum[i] = sum[i - 1] + 1.0 / i;}auto calc = [&](ll x) {if(x <= 1000000)return sum[x];elsereturn log((long double)x) + 0.5772156649015;};long double ans = 0;REP(i, 1, 1000001) {if(i > n / i) {break;}ans += calc(i) * (calc(n / i) - calc(n / (i + 1)));// cout << i << " " << cnt << endl;}REP(i, 1, 1000001) {if(i >= n / i) {break;}ans += calc(n / i) / i;// cout << i << " " << endl;}cout << fixed << setprecision(15) << ans << endl;return 0;}