hdu5245(概率)

题意大致是：进行K次染色，每次染色会随机选取一个以(x1,y1),(x2,y2)为一组对角的子矩阵进行染色，求K次染色后染色面积的期望值(四舍五入)。

显而易见的是，当K=1时，期望被染色的面积会等于每个1*1的方块被染色的期望累加之和。

假设K=1时即只染色一次时，位于第x行第y列的方块被染色的概率为A[x,y]

在K次操作后被染色的期望假设为P[x,y]，可以用

P[x,y]=1-(1-A[x,y])^k

来计算。

所以就转为了求出A[x,y]。把所有的情况跑出来

由题目描述，一次染色中可能的操作有n^2*m^2种

计算A[x,y]时，我们可以把整个矩阵做如下拆分

当前计算的方块为[x,y]，即图中编号为5的部分

将其他部分拆分成图上8个区域，则可得到以下关系

对于一种染色方案能够覆盖方块[x,y]时
①[x1,y1]取在区域1内时，[x2,y2]可以在5、6、8、9四个区域内任取;
②[x1,y1]取在区域2内时，[x2,y2]可以在4、5、6、7、8、9六个区域内任取;
③[x1,y1]取在区域3内时，[x2,y2]可以在4、5、7、8四个区域内任取;
④[x1,y1]取在区域4内时，[x2,y2]可以在2、3、5、6、8、9六个区域内任取;
⑤[x1,y1]取在区域5内时，[x2,y2]可以在所有区域内任取;
⑥[x1,y1]取在区域6内时，[x2,y2]可以在1、2、4、5、7、8六个区域内任取;
⑦[x1,y1]取在区域7内时，[x2,y2]可以在2、3、5、6四个区域内任取;
⑧[x1,y1]取在区域8内时，[x2,y2]可以在1、2、3、4、5、6六个区域内任取;
⑨[x1,y1]取在区域1内时，[x2,y2]可以在1、2、4、5四个区域内任取;

按照这个关系，即可推出A[x,y]的表达式。

P.S.:本题因为计算过程中会出现n^2*m^2大小的计算，因此需要注意int溢出的问题

参考：http://www.itkeyword.com/doc/7955111441277446327/hdu-5245-joyful

#include<bits/stdc++.h>

using namespace std;
typedef long long ll;
typedef pair<ll,int> P;
#define fi first
#define se second
#define INF 0x3f3f3f3f
#define clr(x,y) memset(x,y,sizeof x)
#define PI acos(-1.0)
#define ITER set<int>::iterator
const int Mod = 1e9 + 7;

const int maxn = 500 + 10;
double fun(double x,int k)
{
    double ret = 1.0;
    for(int i = 1;i <= k;i ++)ret *= x;
    return ret;
}
int main()
{
    int n,m,k;int Tcase;scanf("%d",&Tcase);
    for(int ii = 1;ii <= Tcase; ii ++)
    {
        scanf("%d%d%d",&n,&m,&k);
        double ans = 0;
        for(int i = 1;i <= n;i ++)
        {
            for(int j = 1;j <= m;j ++)
            {
                double p = 0;
                p += 1ll * (i - 1) * (j - 1) * (n - i + 1) * (m - j + 1);
                p += 1ll * (i - 1) * (n - i + 1) * m;
                p += 1ll * (i - 1) * (m - j) * (n - i + 1) * j;

                p += 1ll * (j - 1) * (m - j + 1)* n;
                p += 1ll * n * m;
                p += 1ll * (m - j) * (j ) * n;

                p += 1ll *(n - i) * (j - 1) * i * (m - j + 1);
                p += 1ll * (n - i) * i * m;
                p += 1ll * (n - i) * (m - j) * i * j;
//                cout << i << " " << j << " " << p << endl;
                p = p / n/n/m/m;
                ans += 1 - fun(1 - p,k);
            }
        }
        printf("Case #%d: %d\n",ii,(int)round(ans));
    }
    return 0;
}