BZOJ 2969 浅谈矩阵期望动态规划面积分布式转移

世界真的很大
这道题还是想了好一会儿
。。。思路都是对的。。忽略了K为0的情况。。
以后凡是乘方还是写快速幂吧，又快还能回避K为0的情况

看题先：

description：

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

input：

第一行是整数K，W，H

output：

一行，为答案，四舍五入保留到整数。

首先做期望一定要把 “什么什么的期望” 就当成“什么什么”来想
然后考虑把大的期望拆分成小的期望，利用期望的线性性
这道题首先看到期望被染黑的格子数，等价于有多少个格子期望被染黑，等价于每个格子最后被染黑的概率 * 1累加
每个格子被染黑只需要在K次里面一次被染黑就行了
等价于1减去K次一次都没有被染黑的概率

在一次里面没有被染黑，就是说选择的两个点都在这个点的同一侧，都在左边，都在右边，都在上边，都在下边

考虑选一次在左边的概率，就是左边的面积/总的面积，而又由于两次可以选同一个点，所以直接就是选一次的平方就好

但是这样会算重，比如同在左边的情况和同在上边的情况就算重了同在左上角的情况

减去这四个角的概率就好

然后连续K次选择都没选到，就是一次没有被选到的K次方
1减去这个就是至少一次被选到的概率
累加就好。
我是算出一次的再乘以K-1遍，这样要是K为0就会多算一次
改成快速幂就A了

完整代码：

#include<stdio.h>
#include<algorithm>
using namespace std;

int K,W,H;
double ans=0,n,m;

double quickmub(double a,int b)
{
    double rt=1.0;
    while(b)
    {
        if(b&1)
            rt=rt*a;
        a=a*a;
        b>>=1;
    }
    return rt;
}

double getans(double x,double y)
{
    double rt=0;
    rt+=(x-1)/n * (x-1)/n;
    rt+=(n-x)/n * (n-x)/n;
    rt+=(y-1)/m * (y-1)/m;
    rt+=(m-y)/m * (m-y)/m;
    rt-=(x-1)*(y-1)/(n*m) * (x-1)*(y-1)/(n*m);
    rt-=(n-x)*(m-y)/(n*m) * (n-x)*(m-y)/(n*m);
    rt-=(x-1)*(m-y)/(n*m) * (x-1)*(m-y)/(n*m);
    rt-=(n-x)*(y-1)/(n*m) * (n-x)*(y-1)/(n*m);
    return 1-quickmub(rt,K);
}

int main()
{
    scanf("%d%d%d",&K,&W,&H);
    n=W*1.0,m=H*1.0;
    for(int i=1;i<=W;i++)
        for(int j=1;j<=H;j++)
            ans+=getans(i,j);
    printf("%0.0lf\n",ans);
    return 0;
}
/*
EL PSY CONGROO
*/

嗯，就是这样