Lights inside 3D Grid（概率&期望&二项式）

题目：https://vjudge.z180.cn/problem/LightOJ-1284#author=0

题意：给一个三维空间xyz，每一个位置有一个灯，初始灯是灭的，给一个k，进行k回合选择，，每次随机选出俩个点，这俩个点形成的长方体空间内的每一个灯都要变换状态，包括选出的点及其边缘，问经过k回合后，还亮着的灯的数量的期望。

题解：计算出每一个灯还亮着的期望，他们的和就是总的亮着的期望。首先计算每个点在一个回合内被选到的概率：计算被选到的概率的时候xyz相互独立不影响，分别计算出概率相乘就是被包含的概率：1-((xi-1)(xi-1)+(X-xi)(X-xi))/(X*X),即计算不包含这个点的概率再用1减去即可，，计算出每个点被选中的概率后，再计算其期望：灯初始是灭的，所以当被包含次数是奇数次的时候它就是亮的，即在k回合中所有被选中奇数次的和就是这个点的期望。
这里用到一个公式：
二项分布：

当q=1-p时

#include <algorithm>
#include <bitset>
#include <cmath>
#include <cstdio>
#include <cstdlib>
#include <cstring>
#include <deque>
#include <functional>
#include <iostream>
#include <map>
#include <queue>
#include <set>
#include <stack>
#include <string>
#include <vector>
//#include <unordered_map>
//#include <unordered_set>
//#include <bits/stdc++.h>
//#define int long long
#define pb push_back
#define PII pair<int, int>
#define mpr make_pair
#define ms(a, b) memset((a), (b), sizeof(a))
#define x first
#define y second
typedef long long ll;
const int inf = 0x3f3f3f3f;
const ll INF = 0x3f3f3f3f3f3f3f3f;
const int N = 1e6 + 7;
using namespace std;

int main(int argc, char const *argv[]) {
    int T;
    cin >> T;
    int Case = 0;
    while (T--) {
        printf("Case %d: ", ++Case);
        int x, y, z, k;
        cin >> x >> y >> z >> k;
        double ans = 0;
        for (int i = 1; i <= x; i++) {
            for (int j = 1; j <= y; j++) {
                for (int f = 1; f <= z; f++) {
                    //p用来计算概率
                    double p = 1 - ((i - 1) * (i - 1) + (x - i) * (x - i)) *
                                       1.0 / (x * x);
                    p = p * (1 - ((j - 1) * (j - 1) + (y - j) * (y - j)) * 1.0 /
                                     (y * y));
                    p = p * (1 - ((f - 1) * (f - 1) + (z - f) * (z - f)) * 1.0 /
                                     (z * z));
                    //用公式计算期望
                    ans = ans + (1 - pow(1 - 2 * p, k)) / 2.0;
                }
            }
        }
        printf("%.10lf\n", ans);
    }
    return 0;
}