LightOJ - 1284（確率）

タイトル情報元：クリックして
タイトルの意味：ポイントの3次元空間には、スイッチの状態を持っています。オープン状態の期待を尋ねるのポイントの数。
光スイッチの数を計算した後、YをX、Zは、範囲1~100である。確かにポイントが算出された端ありません。2つだけの状態ので、関連する組み合わせの数を検討します。ポイントカット数回含まれているから計算することができます。
pは、そうすることによって選択された点の確率は、であるように選択間隔である1-iはII 2 IXが、繰り返すことができる1を減算することによって計算することができます。
二項定理（+ b）は^ n及び（AB）^ 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;
}