#include <iostream>
#include <cstdio>
#include <cstring>
#include <algorithm>
#define LL long long
using namespace std;
int main(){
LL n,m; int k;
int T,icase=0;
scanf("%d",&T);
while(T--){
scanf("%lld%lld%d",&n,&m,&k);
double output=0;
for(LL i=1;i<=n;i++){
for(LL j=1;j<=m;j++){
LL ans=(i-1)*(i-1)*m*m;
ans+=(n-i)*(n-i)*m*m;
ans+=n*n*(j-1)*(j-1);
ans+=n*n*(m-j)*(m-j);
ans-=(i-1)*(i-1)*(j-1)*(j-1);
ans-=(n-i)*(n-i)*(j-1)*(j-1);
ans-=(i-1)*(i-1)*(m-j)*(m-j);
ans-=(n-i)*(n-i)*(m-j)*(m-j);
double p=ans*1.0/(n*n*m*m);
double tmp=1.0;
for(int c=1;c<=k;c++)
tmp*=p;
output+=(1.0-tmp);
}
}
printf("Case #%d: %.0f\n",++icase,(output));
}
return 0;
}
原文:http://www.cnblogs.com/jie-dcai/p/4539423.html