LightOJ - 1284 Lights inside 3D Grid（数学期望+组合数学）

题目链接：Lights inside 3D Grid - LightOJ 1284 - Virtual Judge (ppsucxtt.cn)

题意：给你一个三维矩阵，矩阵上的每一个点有一个灯，你可以操作k次，每次操作选取两个点，在这两个点之间的灯状态都会反转（原来亮着的现在灭了，原来灭着的现在亮了），问k次操作后亮着的灯的数目期望值是多少。

分析：这个数目期望值显然是每个灯最后亮着的概率的和，现在问题就转化为了每个灯最后亮着的概率是多少，最后亮代表了这个灯在k次操作内被操作了奇数次，想解决这个问题又面临着一个直接问题，那就是每个位置的灯每次被操作的概率是多少，我们现在来看一下位置在（ i , j , k ）的灯每次被操作的概率是多少，我们每次可以任意选取两个点的坐标，因为x，y，z的选取是相互独立的，所以我们可以分别求取这三个方向的概率，最后乘起来即可，那我们先以x为例，x方向的坐标选取是在（1，X）区间的，那怎么样才能包含i呢？直接考虑情况太多了，我们就来考虑他的对立面，如何才能不包含i呢？是不是两个点都在区间（1，i-1）或者两个点都在区间（i+1，X），所以横坐标包含i的概率就是 1- ( ( i - 1 ) * ( i - 1 ) + ( x - i ) * ( x - i ) ) / ( x * x ),y方向和z方向同理，最后将三个方向的概率乘起来就得到了包含这个点的概率了，好了每个点被操作的概率求出来了，接下来就要看每个点被操作奇数次的概率了，分析，我写在了一张纸上，大家可以参考一下：

 这样，我们就求出来了操作k次后，每个灯亮着的概率，最后直接把每个灯亮着的概率加起来就得到答案了

下面是代码：

#include<cstdio>
#include<cstring>
#include<algorithm>
#include<iostream>
#include<queue>
#include<cmath>
using namespace std;
double f(int i,int x)
{
	return (1-1.0*((i-1)*(i-1)+(x-i)*(x-i))/(x*x));
}
int main()
{
	int T;
	cin>>T;
	for(int o=1;o<=T;o++)
	{
		int x,y,z,K;
		double ans=0;
		scanf("%d%d%d%d",&x,&y,&z,&K);
		for(int i=1;i<=x;i++)
		for(int j=1;j<=y;j++)
		for(int k=1;k<=z;k++)
		{
			double p=f(i,x)*f(j,y)*f(k,z);//(i,j,k)每次被操作的概率 
			ans+=(1-pow(1-2*p,K))/2; 
		}
		printf("Case %d: %lf\n",o,ans);
	}
	return 0;
}