hdu5245 Joyful（概率dp）

hdu5245

题目

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

思路

http://blog.csdn.net/Tc_To_Top/article/details/47681573
考虑某个小方格不被染色的概率，就是所选的矩形区域在他的上下左右的情况，不过四个角会被多算，所以减去，乘以k次表示染k次色，再用1减去就是被染色的概率，累加所有小方格就是结果。

代码

#include <cstdio>
#include <cstdlib>
#include <cstring>
#include <algorithm>
#include <cmath>
#include <queue>

using namespace std;
typedef long long ll;

ll m,n,k;

int main()
{
    int T,kase=1;
    scanf("%d",&T);
    while(T--)
    {
        scanf("%I64d %I64d %I64d",&m,&n,&k);
        ll sum=(ll)m*m*n*n;
        double num=0;
        for(int i=1; i<=m; i++)
            for(int j=1; j<=n; j++)
            {
                ll ans=0;
                ans+=(ll)(j-1)*(j-1)*m*m;
                ans+=(ll)(n-j)*(n-j)*m*m;
                ans+=(ll)(i-1)*(i-1)*n*n;
                ans+=(ll)(m-i)*(m-i)*n*n;
                ans-=(ll)(j-1)*(j-1)*(i-1)*(i-1);
                ans-=(ll)(n-j)*(n-j)*(m-i)*(m-i);
                ans-=(ll)(j-1)*(j-1)*(m-i)*(m-i);
                ans-=(ll)(i-1)*(i-1)*(n-j)*(n-j);
                double p=1.0*ans/sum;
                p=pow(p,k);
                num+=1.0-p;
            }
        printf("Case #%d: %d\n", kase++, (int)round(num));
    }
    return 0;
}