题目来源:click
题意:一个三维空间的的点有开关的状态。问开的状态的点的个数期望。
x,y,z的范围都是1-100.肯定不是计算两端点之后计算开关灯的个数。由于只有两种状态考虑与组合数有关。可以从点进行切入计算被包含了几次。
p就是该点被选中的概率,区间选择为1-i i-x所以乘2 但是可能i-i被重复计算减去1即可。
之后利用二项式定理(a+b)^n与 (a-b)^n进行计算。只有被奇数个的点选择过的才是开的状态相加。
#include<cmath>
#include<iostream>
#include<cstdio>
#include<cstring>
#include<string>
#include<cstdlib>
#include<istream>
#include<vector>
#include<stack>
#include<map>
#include<algorithm>
#include<queue>
#define MAX_len 50100*4
using namespace std;
typedef long long ll;
double dp[110][105][105];
int C(int n)
{
return n*(n-1)/2;
}
double quickpow(double a,int n)
{
double res=1.0;
while(n)
{
if(n&1)
{
res*=a;
}
n>>=1;
a*=a;
}
return res;
}
int main()
{
int T;
scanf("%d",&T);
int yy=1;
while(T--)
{
memset(dp,0,sizeof(dp));
int x,y,z,k,i,j;
scanf("%d %d %d %d",&x,&y,&z,&k);
double ans=0;
for(i=1;i<=x;i++)
{
for(j=1;j<=y;j++)
{
for(int hh=1;hh<=z;hh++)
{
double p=0.0;
p=1.0*((x-i+1)*i*2.0-1)/(x*x*1.0);
p*=1.0*((y-j+1)*(j)*2.0-1)/(y*y*1.0);
p*=1.0*((z-hh+1)*(hh)*2.0-1)/(z*z*1.0);
ans+=(1-quickpow(1-2*p,k))*1.0/(2.0);
}
}
}
printf("Case %d: %.10lf\n",yy++,ans);
}
return 0;
}