HDU 5245

题目大意：

每次随机选择两个点，便把这两个点之间形成的子矩阵上的每一个方块涂色，问随机选择k次，整个m*n的矩阵中有多少个小方块被涂上了颜色

这道题不难，但自己智商实在捉急，一直想不出来。。。

因为这里n,m<=500,所以总共250000个方块，我们可以考虑的是每一个方块在随机选择1次后被染色的概率 p[i][j]

那么k次后就变成了 1-(1-p[i][j])^k的概率了，我们将所有概率相加就得到了总共的染色块数

至于怎么计算被染色概率就是可以将整个矩形块分解，然后去计算所有不能包括当前点的所有点对组合，求出不被染色的次数

总的次数n*m*n*m 会超int

 #include <iostream>
 #include <cstdio>
 #include <cstring>
 #include <algorithm>
 #include <queue>
 using namespace std;
 #define N 505
 #define ll long long
 double p[N][N];
 int n , m , k;
 //计算i,j位置的正方形小块不包含的取点总数
 ll cal(int i , int j)
 {
     ll c1 = i- , c2=(i-)*(m-j) , c3=m-j , c4=(n-i)*(m-j);
     ll c5 = n-i , c6=(n-i)*(j-) , c7=j- , c8=(i-)*(j-);
     ll ret = ;
     ret+=c1*(c1+c2+c8)+c2*(c2+c1+c8+c3+c4)+c3*(c2+c3+c4)+c4*(c2+c3+c4+c5+c6)+c5*(c4+c5+c6)+c6*(c4+c5+c6+c7+c8)+c7*(c6+c7+c8)+c8*(c6+c7+c8+c1+c2);
   //  if(ret<0 || ret>n*m*n*m) cout<<i<<" "<<j<<" "<<ret<<endl;
     return ret;
 }

 double q_pow(double a , int b)
 {
     double ret=;
     while(b){
         if(b&) ret*=a;
         a*=a;
         b>>=;
     }
     return ret;
 }

 int main()
 {
     #ifndef ONLINE_JUDGE
         freopen("a.in" , "r" , stdin);
     #endif
     int T , cas=;
     scanf("%d" , &T);
     while(T--)
     {
         scanf("%d%d%d" , &n , &m , &k);
         double ret = ;
         ll all = (ll)n*m*n*m;
        // cout<<all<<endl;
         for(int i= ; i<=n ; i++)
             for(int j= ; j<=m ; j++){
                 p[i][j] = cal(i,j)*1.0/all;
                // if(p[i][j]>=1||p[i][j]<0) cout<<i<<" "<<j<<" "<<p[i][j]<<endl;
                 p[i][j] = q_pow(p[i][j] , k);
                 ret += -p[i][j];
             }
         printf("Case #%d: %.0f\n" , ++cas , ret);

     }
     return ;
 }