結果

問題 No.3072 Sum of sqrt(x)
ユーザー 👑 hos.lyrichos.lyric
提出日時 2020-09-12 16:33:52
言語 C++14
(gcc 12.3.0 + boost 1.83.0)
結果
WA  
実行時間 -
コード長 1,865 bytes
コンパイル時間 1,123 ms
コンパイル使用メモリ 97,764 KB
実行使用メモリ 6,944 KB
最終ジャッジ日時 2024-06-10 15:31:35
合計ジャッジ時間 57,922 ms
ジャッジサーバーID
（参考情報） 		judge2 / judge5
このコードへのチャレンジ
（要ログイン）

テストケース

テストケース表示
入力 結果 実行時間
実行使用メモリ
testcase_00 AC 2 ms
5,248 KB
testcase_01 WA -
testcase_02 AC 3 ms
5,376 KB
testcase_03 AC 3 ms
5,376 KB
testcase_04 AC 2 ms
5,376 KB
testcase_05 AC 3 ms
5,376 KB
testcase_06 WA -
testcase_07 WA -
testcase_08 AC 495 ms
5,376 KB
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 WA -
testcase_23 WA -
testcase_24 WA -
testcase_25 WA -
testcase_26 TLE -
testcase_27 WA -
testcase_28 WA -
testcase_29 WA -
権限があれば一括ダウンロードができます

ソースコード

diff # 

#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];
  int N;
  for (; ~scanf("%d", &N); ) {
    Int upp = 0;
    long double low = 0.0;
    for (int i = 0; i < N; ++i) {
      Int x;
      scanf("%lld", &x);
      if (x > 0) {
        /*
          t = sqrt(x) - k
          f(t) = (t + k)^2 - x
          f'(t) = 2 (t + k)
          t <- t - ((t + k)^2 - x) / (2 (t + k)) = (t^2 + x - k^2) / (2 (t + k))
        */
        const long double y = sqrtl((long double)x);
        const Int k = (Int)y;
        long double t = k;
        t = (t * t + (x - k * k)) / (2.0L * (t + k));
        t = (t * t + (x - k * k)) / (2.0L * (t + k));
        t = (t * t + (x - k * k)) / (2.0L * (t + k));
        t = (t * t + (x - k * k)) / (2.0L * (t + k));
        t = (t * t + (x - k * k)) / (2.0L * (t + k));
// cerr<<k<<" "<<t<<endl;
        upp += k;
        low += t;
      }
      sprintf(buf, "%.18Lf", upp + low);
      buf[18] = '\0';
      puts(buf);
    }
  }
  return 0;
}
0