Problem:
次のようなゲームを考える． 1 から n までの数が 1 つずつ書かれた n 枚のカードが k 組ある．これら kn 枚のカードをよくシャッフル（よく切ること）して， k 枚ずつの山を作り横一列に並べる．このようにしてできる n 個の山の左から i 番目の（k 枚のカードの）山を「山 i 」と呼ぶことにする．

ゲームは山 1 から始める．山の一番上のカード 1 枚を引き（引いたカードは元の山に戻さない），そのカードに書かれていた数が i だった場合には山 i の一番上のカード 1 枚を引く．このようにして，引いたカードに書かれていた数を番号とする山の一番上のカード 1 枚を引くことを繰り返し，すべての山にカードが無くなれば成功である．まだカードが残っている山があるのに，次にカードを引くべき山が無くなっていた場合は失敗である．

　途中で失敗した場合には，そのまま失敗で終了するか，または残ったカードの山をそのまま（山の番号もそのまま）にしてゲームを再開する．ゲームを再開する場合は，最初に引くカードはカードが残っている山のうちの一番左の山からとする（その山の一番上のカードが最初に引かれるカードとなる）．再開後も再開前と同様の方法でゲームを進め，すべての山にカードが無くなれば成功であり，まだカードが残っている山があるのに，次にカードを引くべき山が無くなった場合は失敗である．

このようなゲームの再開を最大 m 回まで行うものとする．ただし，m は 0 か 1 である．つまり， 1 回も再開しないか， 1 回だけ再開するかのいずれかである．ゲーム開始前のシャッフルの仕方によりカードの初期配置は異なる．当然，カードの初期配置により，再開せずに成功することもあれば，再開して成功することも，再開して失敗することもある．十分シャッフルしているので，どの初期配置も全て同じ確率で現れるものと考えることにして，再開が m 回以内で成功する確率 p を求めたい．この確率 p を小数で表し，小数第 r 位まで求めて出力するプログラムを作りなさい．ただし，次の条件を満たすように出力すること.

十分大きい正整数 K を取ると p×10K が 整数となる場合， 小数部は途中から 0 が続くが，その 0 も出力すること． 例えば， p = 3/8 = 0.375 の場合， r = 5 なら 0.37500 と出力し， r = 2 なら 0.37 と出力する． p = 1.0 の場合も同様に， 例えば r = 3 なら 1.000 と出力すること．

例えば 0.150000… は循環小数 0.1499999… として表すこともできるが， このような場合， 前者の表し方を用いる．

入力ファイルの 1 行目には整数 n，k，m，r がこの順に空白を区切り文字として書いてある． 1 ≦ n ≦ 10000， 1 ≦ k ≦ 100， m = 0 または m = 1， 1 ≦ r ≦ 10000 である．

入力例１

入力例２

入力例３

2 1 0 5

3 1 1 3

2 2 1 3

　

出力例１

出力例２

出力例３

0.50000

0.833

1.000

入力

入力は複数のデータセットからなる．n, k, m, r がすべて 0 のとき入力が終了する．データセットの数は 5 を超えない．

出力

データセットごとに、指定通りに p を１行に出力する．

入力例

2 1 0 5
3 1 1 3
2 2 1 3
0 0 0 0

出力例

0.50000
0.833
1.000

各データセットについて，指定通りに出力した p の後に改行を入れること．

上記問題文と自動審判に使われるデータは、情報オリンピック日本委員会が作成し公開している問題文と採点用テストデータです。

Solution in C++ (judge verdict: Time Limit Exceeded)

#include <cstdio>
#include <cstring>
#include <cmath>
#include <vector>
#include <list>
#include <string>
#include <set>
#include <map>
#include <queue>
#include <stack>
#include <algorithm>
#include <iostream>
#include <sstream>
using namespace std;

typedef long long int64;

int64 gcd(int64 a, int64 b)
{
	if(a > b)
		swap(a, b);
		
	if(a == 0LL)
		return b;
		
	return gcd(b % a, a);
}

int64 lcm(int64 a, int64 b)
{
	return a / gcd(a, b) * b;
}

struct number {

	int64 x, y;	// x / y
	
	number() : x(0LL), y(1LL) {}
	number(int64 x, int64 y) : x(x), y(y) {}
	
	void reg()
	{
		int64 d = gcd(x, y);
		x /= d;
		y /= d;
	}
	
	void print(int64 r)
	{
		int64 cx = x;
		
		printf("%lld.", x / y);
		cx = (x - x / y) * 10LL;
	
		for(int64 i = 0; i < r; ++i) {
		
			int64 d = cx / y;
			cx = (cx - d * y) * 10LL;
			
			printf("%lld", d);
		}
	}
	
	number operator+(const number &rhs)
	{
		int64 ny = lcm(this->y, rhs.y);
		int64 nx = this->x * (ny / this->y) + rhs.x * (ny / rhs.y);
		
		number ret(nx, ny);
		
		ret.reg();
		
		return ret;
	}
	
	number operator*(const number &rhs)
	{
		int64 ny = this->y * rhs.y;
		int64 nx = this->x * rhs.x;
		
		number ret(nx, ny);
		
		ret.reg();
		
		return ret;
	}
};


number calc1(int64 n, int64 k)
{
	return number(1LL, n);
}

number calc2(int64 n, int64 k)
{
	if(n == 1)
		return number(0LL, 1LL);

	number a(n - 1LL, n), b(1LL, n);

	number ret = calc1(n - 1LL, k) * a + calc2(n - 1LL, k) * b;

	/*printf("[debug] calc2(%lld, %lld) = ", n, k);
	ret.print(5);
	printf("\n");*/

	return ret;
}

number calc(int64 n, int64 m, int64 k, int64 r)
{
	number ans = calc1(n, k);
	
	if(m == 1)
		ans = ans + calc2(n, k);
	
	ans.print(r);
	printf("\n");
}

int main()
{
	while(true) {
	
		int n, k, m, r;
		
		scanf("%d%d%d%d", &n, &k, &m, &r);
		if(n == 0)
			break;
	
		calc(n, m, k, r);
	}

	return 0;
}

/* はやく人間になりたい */
/* 人権を得たい */
/* nullmineralが書きましたが */
/* 責任はまったくとりません */