problem
题目描述
有两个可重集合 nnn,初始时 mmm 只包含一个 000,是给定的。
执行以下操作:
- 在 BBB 中随机选一个数 yyy,把 yyy 从 BBB 移动到 AAA。
- 给答案加上 AAA 的平均值。
- 若 BBB 非空,回到步骤 111。
求最后答案的期望。
输入格式
第一行两个整数 n,mn,mn,m,表示集合 BBB 可以分成 mmm 个大小为 nnn 的部分。
solution
考虑第一个被选择的数 xxx 的贡献,即 1n∗m⋅(12+13+...+1n∗m+1)⋅x\frac{1}{n*m}·(\frac{1}{2}+\frac{1}{3}+...+\frac{1}{n*m+1})·xn∗m1⋅(21+31+...+n∗m+11)⋅x 【从加入后每一轮都会产生一个贡献】
发现系数具有连续性,第 iii 次加入的数的贡献为 1n∗m⋅(1i+1+1i+2+...+1n∗m+1)⋅x\frac{1}{n*m}·(\frac{1}{i+1}+\frac{1}{i+2}+...+\frac{1}{n*m+1})·xn∗m1⋅(i+11+i+21+...+n∗m+11)⋅x
考虑枚举这个数是第几次被加入的,统计答案即可。
相同权值的数的所有情况贡献是完全相同的。
线性递推逆元,并处理逆元的后缀和即可。
code
#include <cstdio>
#define maxn 20000005
#define mod 998244353
int n, m;
int inv[maxn];int main() {freopen( "mos.in", "r", stdin );freopen( "mos.out", "w", stdout );scanf( "%d %d", &n, &m );inv[1] = 1;for( int i = 2;i <= n * m + 1;i ++ ) inv[i] = ( - mod / i * 1ll * inv[mod % i] % mod + mod ) % mod;int p = 0, f = 0;for( int i = n * m + 1;i > 1;i -- ) {f = ( 1ll * f + inv[i] ) % mod;p = ( 1ll * p + f ) % mod;} int ans = 0;for( int i = 1, x;i <= n;i ++ ) {scanf( "%d", &x );ans = ( 1ll * ans + 1ll * p * x % mod * m ) % mod;}printf( "%d\n", 1ll * ans * inv[n * m] % mod );
}