#include "bits/stdc++.h" using namespace std; #define FOR(i,j,k) for(int (i)=(j);(i)<(int)(k);++(i)) #define rep(i,j) FOR(i,0,j) #define each(x,y) for(auto &(x):(y)) #define mp make_pair #define mt make_tuple #define all(x) (x).begin(),(x).end() #define debug(x) cout<<#x<<": "<<(x)< pii; typedef vector vi; typedef vector vll; /* trigamma(z) = Σ(k_0_inf, B[k]/z^(k+1)) trigamma(z) = trigamma(z+1)+1/(z^2) なのでこれを繰り返すと trigamma(z+1e8)は初項以外は十分小さい。 なのでこれを繰り返して 十分大きなz'をつくる。 */ double trigamma(double z) { static double B[7] = {1.0 / 1,-1.0 / 2,1.0 / 6,0.0,-1.0 / 30,0.0,1.0 / 42}; double res = 0; for(int i = 0; i < 100000; ++i) { res += 1 / (z*z); z += 1; } double x = 1; rep(i, 7) { x /= z; res += B[i] * x; } return res; } int main(){ ios::sync_with_stdio(false); cin.tie(0); double x; while(cin >> x) { double ans = trigamma(1+x); cout << setprecision(20) << ans << endl; } }