yongpython 回答篇幅：简单明了

```pascal program solve; var y, left, right, mid, x0: double; iter: integer; ans: int64; function S(x: double): double; var k: integer; sum: double; begin sum := 0.0; for k := 0 to 15 do sum := sum + 1.0 / (x + k); S := sum; end; begin readln(y); if y <= 0.0 then begin writeln(0); exit; end; left := 1e-15; right := 1e9; for iter := 1 to 80 do begin mid := (left + right) * 0.5; if S(mid) > y then left := mid else right := mid; end; x0 := (left + right) * 0.5; ans := trunc(1.0 / x0); writeln(ans); end. ``` [2025-12-17 17:56:02 | AI问答 | 238点数解答]