Educational Codeforces Round 109 (Rated for Div. 2) E. Assimilation IV（组合数学，概率轮）

题目链接：https://codeforces.com/contest/1525/problem/E

题目大意：

有个城市，个点，表示城市i到点j的距离。

现在有轮，每一轮可以选择一个未建造过雕塑的城市，建造一个雕塑。雕塑的初始影响力为1，建造后每一轮影响力都会加1，若城市的雕塑的影响力大于等于城市到某个点的距离，那么这个点就被雕塑的影响力覆盖。

现在每一轮随机选择一个城市建筑雕塑，求最后被雕塑影响力覆盖的点的期望个数。

题解：

若我们直接去计算某个点被覆盖的方案数，由于一个点可能被多个城市覆盖，计算的难度比较大。

我们可以反过来求，一个点不被覆盖的方案数，总方案数减去不被覆盖的方案数，就是被覆盖的方案数。总方案数显然是

接下我们具体来看一个点，一定不被覆盖的方案数如何计算。

我们将城市按照距点的距离从小到大排序记为，

我们首先看第1个城市，那么它只能选择在最后个顺序里建造，可以选择的方案数为

然后看第2个城市，那么它只能选择在最后个顺序里建造，由于，因此第1个城市已经占了一个位置，所以还可以选择的方案数为

由此，第个城市，可以选择的方案数为

每个城市可选择的方案数乘起来，就是这个点不被覆盖的方案数。

计算出这个关键的数据，我们就可以根据期望的数学定义，计算最终的答案。

代码如下：

#include<bits/stdc++.h>

using namespace std;
const int nn =5100;
const int inff = 0x3fffffff;
const double eps = 1e-8;
typedef long long LL;
const double pi = acos(-1.0);
const LL mod = 998244353;
int n,m;
LL POW(LL x,LL y)
{
    LL ret=1;
    while(y)
    {
        if(y&1)
            ret=(ret*x)%mod;
        x=(x*x)%mod;
        y/=2;
    }
    return ret;
}
int d[25][51000];
int main()
{
    scanf("%d%d",&n,&m);
    for(int i=1;i<=n;i++)
    {
        for(int j=1;j<=m;j++)
            scanf("%d",&d[i][j]);
    }
    LL jc=1;
    for(int i=1;i<=n;i++)
        jc=(jc*i)%mod;
    LL tem=0;
    for(int i=1;i<=m;i++)
    {
        vector<int>ve;
        for(int j=1;j<=n;j++)
            ve.push_back(d[j][i]);
        sort(ve.begin(),ve.end());
        LL num=1;
        for(int j=0;j<n;j++)
        {
            LL chose = max(0,ve[j]-1-j);
            num=(num*chose)%mod;
        }
        tem=(tem+num)%mod;
    }
    LL ans=(m*jc-tem)%mod;
    ans=(ans+mod)%mod;
    ans=(ans*POW(jc,mod-2))%mod;
    cout<<ans<<endl;
    return 0;
}