LightOJ-1284 Lights inside 3D Grid

单独计算每个点被选中奇数次的期望，加起来就是总亮灯数的期望

对于一次选取，x轴位置 i 的点，被选中的概率就是1.0减去不被选中的概率（两个点都在 i 一侧）

P(x,i)=1.0-1.0*((i-1)*(i-1)+(x-i)*(x-i))/(x*x)

三维坐标下被选中的概率就是三个维度概率乘起来p=P(x,i)*P(y,j)*P(z,k);

假设选取n次，f(n)为选取n次被选中奇数次的概率，g(n)为选区n次被选中偶数次的概率

f(n)=p*g(n-1)+(1-p)*f(n-1)

g(n)=p*f(n-1)+(1-p)*g(n-1)

f(n)+g(n)=1

得f(n)=(1-2p)*f(n-1)+p，f(1)=p

展开后为f(n)=(1-2p)^(n-1)*p+(1-2p)^(n-2)*p+......+(1-2p)*p+p

就是一个等比数列，第一项为p，用等比求和公式得f(n)=0.5-0.5*(1-2p)^n

f(n)就是这个点灯亮的期望

#include<iostream>
#include<cstdio>
#include<cstring>
#include<algorithm>
#include<cmath>
using namespace std;
typedef long long ll;
double mypow(double a,int b)
{
	double ret=1.0;
	while(b)
	{
		if(b&1) ret*=a;
		a*=a;
		b>>=1;
	}
	return ret;
}
double P(int X,int x)
{
	return 1.0-1.0*((x-1)*(x-1)+(X-x)*(X-x))/(X*X);
}
int main()
{
	int T,kase=0;
	int x,y,z,K;
	scanf("%d",&T);
	while(T--)
	{
		scanf("%d%d%d%d",&x,&y,&z,&K);
		double ans=0.0;
		for(int i=1;i<=x;i++)
			for(int j=1;j<=y;j++)
				for(int k=1;k<=z;k++)
				{
					double p=P(x,i)*P(y,j)*P(z,k);
					ans+=0.5-0.5*mypow(1.0-2*p,K);
				}
		printf("Case %d: %.6f\n",++kase,ans);
	}
	return 0;
}