概率期望--牛客noip赛前集训（提高组第三场）T1

传送门

理解题意就花了一个多小时···

然后其实也很好推的，首先肯定是一圈一圈清理的，也就是每个块都有一个清理时间的上限

然后我们就可枚举清理的时间，考虑一个块如果一定在k秒被清理，那么以它为中心的以k-1为半对角线那个菱形一定是全部堵塞的

但是这个概率算的时候不太好算

考场上就是因为算概率的方法比较复杂所以最后没有写完（我不会承认我打字慢的）

算的时候可以算>k步的概率，用1减去就是至少用k步的概率

再根据那个公式把<=k的概率全加起来就是恰好用k步的期望

而且算的时候可以一圈一圈的算

具体看代码注释吧

#include<iostream>
#include<cstdio>
#include<algorithm>
#include<cstring>
#include<cmath>
#define maxn 205
#define LL long long
using namespace std;
int n,m;
const int mod=1e9+7;
LL a[maxn][maxn];
 
inline int rd(){
    int x=0,f=1;char c=' ';
    while(c<'0' || c>'9') {if(c=='-')f=-1;c=getchar();}
    while(c<='9' && c>='0') x=x*10+c-'0',c=getchar();
    return x*f;
}
 
inline LL qpow(LL x,int k){
    LL ret=1;
    while(k){
        if(k&1) (ret*=x)%=mod;
        (x*=x)%=mod; k>>=1;
    } return ret;
}
 
int main(){
    n=rd(); m=rd();
    for(int i=1;i<=n;i++)
        for(int j=1;j<=m;j++){
            a[i][j]=rd(); int x=rd();
            a[i][j]=a[i][j]*qpow(x,mod-2)%mod;
        }
    for(int i=1;i<=n;i++)
        for(int j=1;j<=m;j++){
            int x1=min(i,n-i+1),x2=min(j,m-j+1);//最多用几步清除
            LL ans=0,x=a[i][j];
            for(int k=1;k<=min(x1,x2);k++){//枚举这个菱形的大小
                LL pi=1;//pi表示必须用>k步的概率
                for(int t=i-k;t<=i+k;t++){
                    int l=abs(i-t);
                    if(l<k)
                        (pi*=a[t][j+k-l]*a[t][j-k+l]%mod)%=mod;
                    else (pi*=a[t][j])%=mod;
                }
                (ans+=(1+mod-pi)*x%mod*k%mod)%=mod;//1-pi是最多用k步的概率
                (x*=pi)%=mod;//菱形扩大时里面的也要乘上
            }
            printf("%lld",ans);
            if(j==m) printf("\n"); else printf(" ");
        }
    return 0;
}