LightOJ - 1284(概率)

题目来源：click
题意：一个三维空间的的点有开关的状态。随机选取两个点之间的所有点都改一次，问k次开的状态的点的个数期望。
x,y,z的范围都是1-100.肯定不是计算两端点之后计算开关灯的个数。由于只有两种状态考虑与组合数有关。可以从点进行切入计算被包含了几次。
p就是该点被选中的概率，区间选择为1-i i-x所以乘2 但是可能i-i被重复计算减去1即可。
奇数次的才是开，C(x,k) * p^x * (1-p)^(k-x)，当x为奇数的时候成立， 之后利用二项式定理(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;
}