Problem:
Problem F: Ben Toh

As usual, those who called wolves get together on 8 p.m. at the supermarket. The thing they want is only one, a box lunch that is labeled half price. Scrambling for a few discounted box lunch, they fiercely fight every day. And those who are blessed by hunger and appetite the best can acquire the box lunch, while others have to have cup ramen or something with tear in their eyes.

A senior high school student, Sato, is one of wolves. A dormitry he lives doesn't serve a dinner, and his parents don't send so much money. Therefore he absolutely acquire the half-priced box lunch and save his money. Otherwise he have to give up comic books and video games, or begin part-time job.

Since Sato is an excellent wolf, he can acquire the discounted box lunch in 100% probability on the first day. But on the next day, many other wolves cooperate to block him and the probability to get a box lunch will be 50%. Even though he can get, the probability to get will be 25% on the next day of the day. Likewise, if he gets a box lunch on a certain day, the probability to get on the next day will be half. Once he failed to get a box lunch, probability to get would be back to 100%.

He continue to go to supermaket and try to get the discounted box lunch for n days. Please write a program to computes the expected value of the number of the discounted box lunches he can acquire.

Input

Input consists of several datasets.

Input for a single dataset is given as a single integer n.

Input terminates with a dataset where n = 0.

Output

For each dataset, write a line that contains an expected value. You may print any number of digits after the decimal point. Answers that have an error less than 1.0e-2 will be accepted.

Constraints

1 ≤ n ≤ 100,000

Sample Input

1
2
3
0

Output for the Sample Input

1.00000000
1.50000000
2.12500000

Solution in C++11 (judge verdict: Accepted)

#include <stdio.h>
#include <cctype>
#include <limits.h>
#include <cmath>
#include <complex>
#include <bitset>
#include <vector>
#include <map>
#include <set>
#include <stack>
#include <queue>
#include <cstring>
#include <string>
#include <sstream>
#include <algorithm>
#include <iomanip>
#include <iostream>

#define VARIABLE(x) cerr << #x << "=" << x << endl
#define BINARY(x) static_cast<bitset<16> >(x)
#define rep(i,n) for(int i=0;i<(int)(n);i++)
#define REP(i,m,n) for (int i=m;i<(int)(n);i++)
#define if_range(x, y, w, h) if (0<=(int)(x) && (int)(x)<(int)(w) && 0<=(int)(y) && (int)(y)<(int)(h))
#define is_in(x, y, w, h) (0<=(int)(x) && (int)(x)<(int)(w) && 0<=(int)(y) && (int)(y)<(int)(h))
#define ALL(a) (a).begin(),(a).end()
int dx[4]={0, -1, 1, 0}, dy[4]={-1, 0, 0, 1};
using namespace std;
typedef long long ll;
typedef vector<vector<int> > VVI;
const int INF = 1e9;
typedef pair<int, int> PII;

const int MAX_N = 100010;
const double THRESH_HOLD = 1e-15;

int main()
{
	int N;
	while (cin >> N, N) {
		double dp[MAX_N] = {0.0};
		dp[0] = 1.0;

		double prob[MAX_N], r_prob[MAX_N], tmp=1.0;
		prob[0] = r_prob[0] = 1.0;
		for (int i=1; i<N; i++) {
			tmp *= 0.5;
			prob[i] = prob[i-1] * tmp;
			r_prob[i] = prob[i-1] * (1.0 - tmp);
			// cerr << "r_prob[" << i << "] = " << r_prob[i] << endl;
		}

		for (int i=1; i<N; i++) {
			dp[i] += (i+1) * prob[i];
			dp[i] += i * r_prob[i];

			for (int j=i-2; j>=0; j--) {
				// cerr << "dp[i] += " << r_prob[i-j-1] * (dp[j] + i - j - 1) << endl;
				if (r_prob[i-j-1] < THRESH_HOLD) break;
				dp[i] += r_prob[i-j-1] * (dp[j] + i - j - 1);
			}
		}

		cout << fixed << setprecision(10);
		cout << dp[N-1] << endl;
	}
}