#include using namespace std; //(1 - p[0] - p[1]) + p[1] * 2 + p[0] * (1 - p[0]) * 2 + 3 * p[0] * p[0] = x //p[1] = x + p[0] - 1 - p[0] * (1 - p[0]) * 2 - 3 * p[0] * p[0] //=x-p[0]*p[0]-1 //p[0] * dp[1]+ 1= dp[2] //p[0] * (dp[2] + 1) + p[1] * (dp[1] + 1) + (1 - p[0] - p[1]) * (dp[0] + 1) //1+p[1]+p[0]*(p[0]+1)=dp[3] //p[0]*dp[3]+p[1]*dp[2]+p[2]*dp[1]+1=dp[4] // double dp[1000010]; int main() { cin.tie(0); ios::sync_with_stdio(false); dp[1] = 1.0000000000000000; dp[2] = 1.0833333333333333; dp[3] = 1.2569444444444444; dp[4] = 1.5353009259259260; dp[5] = 1.6915991512345676; dp[6] = 2.0513639724794235; double p[6]; double sum = 0.0; for (int i = 2; i <= 6; i++) { p[i - 2] = dp[i] - 1.0; for (int j = 0; j < i - 2; j++) { p[i - 2] -= p[j] * dp[i - j - 1]; } sum += p[i - 2]; } p[5] = 1 - sum; for (int i = 7; i <= 1000000; i++) { dp[i] = 1; for (int j = 0; j < 6; j++) { dp[i] += dp[i - j - 1] * p[j]; } } int t; cin >> t; for (int i = 0; i < t; i++) { int n; cin >> n; cout << fixed << setprecision(10) << dp[n] << endl; } return 0; }