#552516 (C++14) No.8072 Sum of sqrt(x)

• 問題: No.8072 Sum of sqrt(x)
• ユーザー: 👑 hos.lyric
• 提出日時: 2020-09-12 18:02:22
• 言語: C++14 (gcc 12.3.0 + boost 1.83.0)
• 結果: AC
• 実行時間: 1,370 ms / 2,000 ms
• コード長: 1,538 bytes
• コンパイル時間: 3,107 ms
• コンパイル使用メモリ: 97,640 KB
• 実行使用メモリ: 6,944 KB
• 最終ジャッジ日時: 2024-06-10 17:07:43
• 合計ジャッジ時間: 124,463 ms
• ジャッジサーバーID （参考情報）: judge5 / judge1

入力	結果	実行時間; 実行使用メモリ
testcase_00	AC	2 ms; 5,248 KB
testcase_01	AC	3 ms; 5,376 KB
testcase_02	AC	2 ms; 5,376 KB
testcase_03	AC	3 ms; 5,376 KB
testcase_04	AC	4 ms; 5,376 KB
testcase_05	AC	3 ms; 5,376 KB
testcase_06	AC	3 ms; 5,376 KB
testcase_07	AC	1,114 ms; 5,376 KB
testcase_08	AC	386 ms; 5,376 KB
testcase_09	AC	528 ms; 5,376 KB
testcase_10	AC	1,159 ms; 5,376 KB
testcase_11	AC	1,187 ms; 5,376 KB
testcase_12	AC	1,165 ms; 5,376 KB
testcase_13	AC	1,317 ms; 5,376 KB
testcase_14	AC	1,277 ms; 5,376 KB
testcase_15	AC	1,181 ms; 5,376 KB
testcase_16	AC	1,253 ms; 6,944 KB
testcase_17	AC	1,134 ms; 6,940 KB
testcase_18	AC	1,212 ms; 6,940 KB
testcase_19	AC	1,267 ms; 5,376 KB
testcase_20	AC	1,370 ms; 6,940 KB
testcase_21	AC	1,214 ms; 5,376 KB
testcase_22	AC	1,276 ms; 5,376 KB
testcase_23	AC	1,225 ms; 5,376 KB
testcase_24	AC	1,221 ms; 5,376 KB
testcase_25	AC	1,186 ms; 5,376 KB
testcase_26	AC	1,224 ms; 5,376 KB
testcase_27	AC	1,113 ms; 6,944 KB
testcase_28	AC	1,110 ms; 5,376 KB
testcase_29	AC	1,164 ms; 5,376 KB

ソースコード

#include <cassert>
#include <cmath>
#include <cstdint>
#include <cstdio>
#include <cstdlib>
#include <cstring>
#include <algorithm>
#include <bitset>
#include <complex>
#include <deque>
#include <functional>
#include <iostream>
#include <map>
#include <numeric>
#include <queue>
#include <set>
#include <sstream>
#include <string>
#include <unordered_map>
#include <unordered_set>
#include <utility>
#include <vector>

using namespace std;

using Int = long long;

template <class T1, class T2> ostream &operator<<(ostream &os, const pair<T1, T2> &a) { return os << "(" << a.first << ", " << a.second << ")"; };
template <class T> void pv(T a, T b) { for (T i = a; i != b; ++i) cerr << *i << " "; cerr << endl; }
template <class T> bool chmin(T &t, const T &f) { if (t > f) { t = f; return true; } return false; }
template <class T> bool chmax(T &t, const T &f) { if (t < f) { t = f; return true; } return false; }



int main() {
  char buf[110];
  fgets(buf, sizeof(buf), stdin);
  const int N = atoi(buf);
  Int upp = 0;
  long double low = 0.0, err = 0.0;
  for (int i = 0; i < N; ++i) {
    fgets(buf, sizeof(buf), stdin);
    const Int x = atoll(buf);
    if (x > 0) {
      /*
        sqrt(x) - k = (x - k^2) / (sqrt(x) + k)
      */
      const long double y = sqrtl((long double)x);
      const Int k = (Int)y;
      const long double z = (x - k * k) / (y + k);
      upp += k;
      const long double tmp = low + z;
      err += z - (tmp - low);
      low = tmp;
    }
    printf("%.17Lg\n", upp + low + err);
  }
  return 0;
}