// TODO: 解答方針不明. #include using namespace std; const int MAX = 1e8; double dp[MAX + 1]; 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); dp[i] = dp[i - 1] + d; } // 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 * dp[MAX]; // 4. 出力 ~ 後処理. printf("%.17lf\n", ans); return 0; }