#552474 (C++14) No.3072 Sum of sqrt(x)

• 問題: No.3072 Sum of sqrt(x)
• ユーザー: 👑 hos.lyric
• 提出日時: 2020-09-12 16:38:37
• 言語: C++14 (gcc 12.3.0 + boost 1.83.0)
• 結果: WA
• 実行時間: -
• コード長: 2,016 bytes
• コンパイル時間: 944 ms
• コンパイル使用メモリ: 97,756 KB
• 実行使用メモリ: 5,376 KB
• 最終ジャッジ日時: 2024-06-10 15:39:25
• 合計ジャッジ時間: 90,584 ms
• ジャッジサーバーID （参考情報）: judge1 / judge5

入力	結果	実行時間; 実行使用メモリ
testcase_00	AC	2 ms; 5,248 KB
testcase_01	AC	2 ms; 5,376 KB
testcase_02	AC	2 ms; 5,376 KB
testcase_03	AC	2 ms; 5,376 KB
testcase_04	AC	3 ms; 5,376 KB
testcase_05	AC	2 ms; 5,376 KB
testcase_06	AC	2 ms; 5,376 KB
testcase_07	WA	-
testcase_08	AC	456 ms; 5,376 KB
testcase_09	AC	823 ms; 5,376 KB
testcase_10	AC	1,801 ms; 5,376 KB
testcase_11	AC	1,818 ms; 5,376 KB
testcase_12	AC	1,821 ms; 5,376 KB
testcase_13	AC	1,903 ms; 5,376 KB
testcase_14	AC	1,860 ms; 5,376 KB
testcase_15	AC	1,787 ms; 5,376 KB
testcase_16	AC	1,821 ms; 5,376 KB
testcase_17	AC	1,755 ms; 5,376 KB
testcase_18	AC	1,821 ms; 5,376 KB
testcase_19	AC	1,865 ms; 5,376 KB
testcase_20	AC	1,933 ms; 5,376 KB
testcase_21	AC	1,746 ms; 5,376 KB
testcase_22	AC	1,803 ms; 5,376 KB
testcase_23	AC	1,732 ms; 5,376 KB
testcase_24	AC	1,806 ms; 5,376 KB
testcase_25	AC	1,868 ms; 5,376 KB
testcase_26	AC	1,787 ms; 5,376 KB
testcase_27	AC	1,732 ms; 5,376 KB
testcase_28	WA	-
testcase_29	AC	1,729 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;
  for (int i = 0; i < N; ++i) {
    fgets(buf, sizeof(buf), stdin);
    const Int x = atoll(buf);
    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));
      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;
}