結果
| 問題 |
No.454 逆2乗和
|
| コンテスト | |
| ユーザー |
@abcde
|
| 提出日時 | 2019-06-01 21:10:06 |
| 言語 | C++14 (gcc 13.3.0 + boost 1.87.0) |
| 結果 |
AC
|
| 実行時間 | 277 ms / 2,000 ms |
| コード長 | 1,133 bytes |
| コンパイル時間 | 1,536 ms |
| コンパイル使用メモリ | 165,944 KB |
| 実行使用メモリ | 6,944 KB |
| 最終ジャッジ日時 | 2024-09-17 20:00:30 |
| 合計ジャッジ時間 | 11,980 ms |
|
ジャッジサーバーID (参考情報) |
judge1 / judge3 |
(要ログイン)
| ファイルパターン | 結果 |
|---|---|
| sample | AC * 3 |
| other | AC * 31 |
ソースコード
// TODO: 解答方針不明.
// TODO: MLE改善.
#include <bits/stdc++.h>
using namespace std;
const int MAX = 1e8;
double cur, bef;
int main() {
// 1. 入力情報取得.
double x;
scanf("%lf", &x);
// 2. 100000000 まで, 計算してみる.
for(int i = 1; i < MAX + 1; i++){
// 今回分を計算.
double d = 1.0 / (i + x + 0.0);
d /= (i + x + 0.0);
cur = bef + d;
// 前回の計算結果に保存.
bef = cur;
}
// 3. 倍率してみる(※とりあえず実装).
// x = 0: 1.6449340668482264 ÷ 1.64493405783457503 = 約 1.000000005479642984513倍
// x = 1.5: 0.490357756100235 ÷ 0.49035774607552685 = 約 1.000000020443662265419倍
// x = 1000: 0.0009995001666667166 ÷ 0.00099949016676597 = 約 1.000010005001628936956倍
// -> 1.000000005479642984513 × (1.0 + x * 1e-8) と仮定して, 倍率をかけてみる.
double eps = 1.000000005479642984513 * (1.0 + (x + 0.0) / MAX);
double ans = eps * cur;
// 4. 出力 ~ 後処理.
printf("%.17lf\n", ans);
return 0;
}
@abcde