//Δ_H
#include<bits/stdc++.h>
using namespace std;
typedef long long LL;
typedef double DB;
const int N = 1111;
const int K = 11;
const int W = 1<<9;
int n,k,w;
DB p[W][K][K],f[N][W][K],s[N];
void solve(int m,int t,DB*p){
	int a[K],i,j,l,o,x;
	DB f[K],g[K],nf[K],ng[K];
	memset(a,0,sizeof(a));
	a[0]=k;
	for(i=1;i<k;i++)
		if(m&(1<<(i-1)))
			a[i]=k-i;
	a[k]=k+t;
	for(i=0;i<=k;i++){
		if(a[i]){
			for(l=a[i];l<a[i]+k;l++){
				for(o=0;o<=k;o++)
					f[o]=0,g[o]=0,nf[o]=0,ng[o]=0;
				f[0]=(DB)1/k;
				for(j=0;j<=k;j++){
					if(a[j]&&j!=i){
						x=max(a[j]+k-max(l+1,a[j]),0);
						for(o=0;o<=k;o++)
							ng[o]+=f[o]*x/k;
						x=(a[j]<=l&&l+1<=a[j]+k);
						for(o=1;o<=k;o++)
							nf[o]+=f[o-1]*x/k,ng[o]+=g[o-1]*x/k;
						x=max(min(l,a[j]+k)-a[j],0);
						for(o=0;o<=k;o++)
							nf[o]+=f[o]*x/k,ng[o]+=g[o]*x/k;
						for(o=0;o<=k;o++)
							f[o]=nf[o],g[o]=ng[o],nf[o]=0,ng[o]=0;
					}
				}
				for(o=1;o<=k;o++)
					p[i]+=f[o]/(o+1);
				for(o=0;o<=k;o++)
					p[i]+=g[o]/(o+1);
			}
		}
	}//cout<<m<<' '<<t<<endl;
	//for(i=0;i<=k;i++)cout<<i<<' '<<p[i]<<endl;cout<<endl;
}
int main(){
	int i,j,l,o,x,y,z;
	scanf("%d%d",&n,&k);
	w=1<<(k-1);
	for(i=0;i<w;i++)
		for(j=1;j<=k;j++)
			solve(i,j,p[i][j]);
	x=0;
	for(i=1;i<=k-1;i++)
		if(n-1-i>=1)
			x|=1<<(i-1);
	f[n-1][x][1]=1;
	for(l=n;l>=1;l--){
		for(i=w-1;i>=0;i--){
			for(j=1;j<=k;j++){
				z=max(l-k,0)+__builtin_popcount(i)+1;//cout<<l<<' '<<i<<' '<<j<<' '<<f[l][i][j]<<endl;
				for(o=1;o<=k-1;o++){
					if(i&(1<<(o-1))){
						f[l][i^(1<<(o-1))][j]+=f[l][i][j]*p[i][j][o];
						s[l-o]+=f[l][i][j]*p[i][j][o]*z;
					}
				}
				x=i,y=l;
				while(!(x&1)&&y>=1){
					y--;
					x>>=1;
					if(y-(k-1)>=1)
						x|=1<<(k-2);
				}
				if(y){
					y--;
					x>>=1;
					if(y-(k-1)>=1)
						x|=1<<(k-2);
					f[y][x][min(l+j-y,k)]+=f[l][i][j]*p[i][j][0];
					f[y][x][min(l-y,k)]+=f[l][i][j]*p[i][j][k];
				}
				s[l]+=f[l][i][j]*p[i][j][0]*z;
				s[l+j]+=f[l][i][j]*p[i][j][k]*z;
			}
		}
	}
	for(i=1;i<=n;i++)
		printf("%.12lf\n",s[i]+1);
	return 0;
}