uva 11605 概率 化简公式

#include<cstdio>
#include<cstring>
#include<cmath>
using namespace std;
double cal(int a,int b){
	return 1.0*(2*a*(b-a+1)-1)/(b*b);
}
double solve(double p,int k){
	return (1-pow(1-2*p,k))/2;
}
int main(){
	int t=1,T,m,n,p,k;
	scanf("%d",&T);
	while(T--){
		double res=0;
		scanf("%d%d%d%d",&m,&n,&p,&k);
		for(int i=1;i<=m;i++){
			double t1=cal(i,m);
			for(int j=1;j<=n;j++){
				double t2=cal(j,n);
				for(int q=1;q<=p;q++){
					double t3=cal(q,p);
					res+=solve(t1*t2*t3,k);
				}
			}
		}
		printf("Case %d: %.10f\n",t++,res);
	}
}