蒟蒻ABC F求条

#include<bits/stdc++.h>
using namespace std;
const int N=39;
int n,t,m,k;
long double f[N][N],g[N];
int main(){
	cin>>n>>t>>m>>k;
	f[t][k]=1;
	for(int i=t-1;i>=0;i--){
		f[i][k]=1;
		for(int j=k-1;j>=0;j--){
			g[0]=n*f[i+1][j];
			for(int c=1;c<=m;c++){
				g[c]=0;
				for(int d=1;d<=min(c,k-j);d++) g[c]=max(g[c],g[c-d]+f[i+1][j+d]-f[i+1][j]);
			}
			f[i][j]=g[m]/n;
		}
	}
	printf("%.9Lf",f[0][0]);
	return 0;
}