bzoj 2969: 矩形粉刷 期望

题意

为了庆祝新的一年到来，小M决定要粉刷一个大木板。大木板实际上是一个W*H的方阵。小M得到了一个神奇的工具，这个工具只需要指定方阵中两个格子，就可以把这两格子为对角的，平行于木板边界的一个子矩形全部刷好。小M乐坏了，于是开始胡乱地使用这个工具。
假设小M每次选的两个格子都是完全随机的（方阵中每个格子被选中的概率是相等的），而且小M使用了K次工具，求木板上被小M粉刷过的格子个数的期望值是多少。
1 ≤ W, H ≤ 1000， 0 ≤ K ≤ 100

分析

分别算出每个格子的贡献然后加起来即可。

代码

#include<iostream>
#include<cstdio>
#include<cstdlib>
#include<cstring>
#include<algorithm>
#include<cmath>
#include <iomanip>
using namespace std;

int n,m,k;

int main()
{
    scanf("%d%d%d",&k,&n,&m);
    double ans=0;
    for (int i=1;i<=n;i++)
        for (int j=1;j<=m;j++)
        {
            double w=(double)1.0*2*i*j/n/m*(n-i+1)*(m-j+1)/n/m
                    +(double)1.0*2*i*(m-j+1)/n/m*j*(n-i+1)/n/m
                    -(double)1.0*2*i*(n-i+1)/n/m/n/m
                    -(double)1.0*2*j*(m-j+1)/n/m/n/m
                    +(double)1.0*1*1/n/m/n/m;
            ans+=1-pow(1-w,k);
        }
    cout<<fixed<<setprecision(0)<<ans<<endl;
    return 0;
}