题目传送门

https://loj.ac/problem/3119

现在 BZOJ 的管理员已经不干活了吗,CTS(C)2019 和 NOI2019 的题目到现在还没与传上去。

果然还是 LOJ 好。

题目

恰好有 \(k\) 个极大数不太好求,我们还是转化成二项式反演。

然后就变成了给定一个点的集合 \(S\),求钦定 \(S\) 中的点是极大点的方案数。可以发现 \(S\) 中的点因为必须要保证没有一维的坐标相同,所以到底是哪些点是不重要的,有用的只有 \(|S|\)。所以问题转化为钦定了 \(k\) 个点必须是极大点的方案数。

我们可以发现一个性质,如果我们把所有的极大点的数字从小到大考虑的话,那么每次选完和当前个点有关的的面上的数字以后,剩下的点上的数字如何安排实际上是一个 \((n-1)\times(m-1)\times(l-1),k-1\) 的子问题了。

我们需要先选出 \(k\) 个被我们钦定的点的排列,它们应该满足没有一维的坐标是相同的,于是就很显然是 \(\binom nk\binom mk\binom lk(k!)^3\)。至于为什么是排列,上面说了,只有从小到大考虑,才满足上面的性质。

然后我们需要给被我们选出来的 \(k\) 个点所在的所有面上的所有位置、剩下来的所有位置安排数字了。剩下来的 \(nml-(n-k)(m-k)(l-k)\) 很好安排,先从所有的 \(nml\) 个数字中选出来那么多,然后自由排列一下就好了。这个是 \(\binom{nml}{(n-k)(m-k)(l-k)}((n-k)(m-k)(l-k))!\)。

主要的难点在于被选出的 \(k\) 个点所在的面上的数字。我们有 \(nmk - (n-k)(m-k)(l-k)\) 个备选数字。因为我们是从小到大排列的,所以被最后选出来的位置上的数字应该是最大的。同时,它还会带走 \((n-k+1)(m-k+1)(l-k+1) - (n-k)(m-k)(l-k)\) 个最后选它的时候的子立方体中和它同一面的数字。此时,因为它是最大的,所以无论选哪些数字都不会影响前一位置的选择。因此它有 \(\frac{(nml-(n-k)(m-k)(l-k)-1)!}{(nml-(n-k+1)(m-k+1)(l-k+1))!}\) 中选择方法。

然后选完以后,倒数第二个被选出来的位置上的数字,应该是剩下的数字中最大的。而剩下的数字一共有 \(nml-(n-k+1)(m-k+1)(l-k+1)\) 种,于是和之前最后一个一样一样考虑就可以了。

因此给被选出的 \(k\) 个点所在的面上安排数字的总方案数为

\[\prod_{i=1}^k \frac{(nml-(n-i)(m-i)(l-i)-1)!}{(nml-(n-i+1)(m-i+1)(l-i+1))!}
\]

可以发现在阶乘上的我们可以消掉一部分,于是就是

\[(nml-(n-k)(m-k)(l-k)-1)!\prod_{i=1}^{k-1} \frac 1{(nml-(n-i+1)(m-i+1)(l-i+1))}
\]

于是把上面的东西都整合一下,钦定了 \(k\) 个点必须是极大点的方案数为

\[\binom nk\binom mk\binom lk(k!)^3\binom{nml}{(n-k)(m-k)(l-k)}((n-k)(m-k)(l-k))!(nml-(n-k)(m-k)(l-k)-1)!\prod_{i=1}^{k-1} \frac 1{(nml-(n-i+1)(m-i+1)(l-i+1))}
\]

\[=\binom nk\binom mk\binom lk(k!)^3\frac{(nml)!}{((n-k)(m-k)(l-k))!(nml-(n-k)(m-k)(l-k))!}((n-k)(m-k)(l-k))!(nml-(n-k)(m-k)(l-k)-1)!\prod_{i=1}^{k-1} \frac 1{(nml-(n-i+1)(m-i+1)(l-i+1))}
\]

\[= \binom nk\binom mk\binom lk(k!)^3{(nml)!}\prod_{i=1}^{k} \frac 1{(nml-(n-i+1)(m-i+1)(l-i+1))}
\]

令 \(f(k)\) 等于上面那个钦定了 \(k\) 个点必须是极大点的方案数,\(g(k)\) 表示恰好 \(k\) 个点的方案数,那么显然有

\[f(k)=\sum_{i=k}^n\binom ikg(i)
\]

所以直接二项式反演就可以了。

因为是求期望,最后记得除掉 \((nml)!\)。


下面是代码。\(\prod\limits_{i=1}^{k} \frac 1{(nml-(n-i+1)(m-i+1)(l-i+1))}\) 可以线性预处理,因此总的时间复杂度为 \(O(n)\)。

#include<bits/stdc++.h>

#define fec(i, x, y) (int i = head[x], y = g[i].to; i; i = g[i].ne, y = g[i].to)
#define dbg(...) fprintf(stderr, __VA_ARGS__)
#define File(x) freopen(#x".in", "r", stdin), freopen(#x".out", "w", stdout)
#define fi first
#define se second
#define pb push_back template<typename A, typename B> inline char smax(A &a, const B &b) {return a < b ? a = b , 1 : 0;}
template<typename A, typename B> inline char smin(A &a, const B &b) {return b < a ? a = b , 1 : 0;} typedef long long ll; typedef unsigned long long ull; typedef std::pair<int, int> pii; template<typename I>
inline void read(I &x) {
int f = 0, c;
while (!isdigit(c = getchar())) c == '-' ? f = 1 : 0;
x = c & 15;
while (isdigit(c = getchar())) x = (x << 1) + (x << 3) + (c & 15);
f ? x = -x : 0;
} const int N = 5e6 + 7;
const int P = 998244353; int n, m, l, k;
int a[N], s[N], f[N]; inline int smod(int x) { return x >= P ? x - P : x; }
inline void sadd(int &x, const int &y) { x += y; x >= P ? x -= P : x; }
inline int fpow(int x, int y) {
int ans = 1;
for (; y; y >>= 1, x = (ll)x * x % P) if (y & 1) ans = (ll)ans * x % P;
return ans;
} int fac[N], inv[N], ifac[N];
inline void ycl(const int &n = ::n) {
fac[0] = 1; for (int i = 1; i <= n; ++i) fac[i] = (ll)fac[i - 1] * i % P;
inv[1] = 1; for (int i = 2; i <= n; ++i) inv[i] = (ll)(P - P / i) * inv[P % i] % P;
ifac[0] = 1; for (int i = 1; i <= n; ++i) ifac[i] = (ll)ifac[i - 1] * inv[i] % P;
} inline int C(int x, int y) {
if (x < y) return 0;
return (ll)fac[x] * ifac[y] % P * ifac[x - y] % P;
} inline void work() {
ycl(std::max(m, l));
int nml = (ll)n * m % P * l % P, ii = 1;
a[0] = 1;
for (int i = 1; i <= n; ++i) a[i] = (nml + P - (ll)(n - i) * (m - i) % P * (l - i) % P) % P, ii = (ll)ii * a[i] % P;
ii = fpow(ii, P - 2);
for (int i = n; i; --i) s[i] = ii, ii = (ll)ii * a[i] % P;
for (int i = 1; i <= n; ++i) f[i] = (ll)C(n, i) * C(m, i) % P * C(l, i) % P * fpow(fac[i], 3) % P * s[i] % P;
int ans = 0;
for (int i = k; i <= n; ++i)
if ((i - k) & 1) sadd(ans, P - (ll)C(i, k) * f[i] % P);
else sadd(ans, (ll)C(i, k) * f[i] % P);
printf("%d\n", ans);
} inline void init() {
read(n), read(m), read(l), read(k);
if (n > m) std::swap(n, m);
if (n > l) std::swap(n, l);
} int main() {
#ifdef hzhkk
freopen("hkk.in", "r", stdin);
#endif
int T;
read(T);
while (T--) {
init();
work();
}
fclose(stdin), fclose(stdout);
return 0;
}

LOJ3119. 「CTS2019 | CTSC2019」随机立方体 二项式反演的更多相关文章

  1. 「CTS2019 | CTSC2019」随机立方体 解题报告

    「CTS2019 | CTSC2019」随机立方体 据说这是签到题,但是我计数学的实在有点差,这里认真说一说. 我们先考虑一些事实 如果我们在位置\((x_0,y_0,z_0)\)钦定了一个极大数\( ...

  2. LOJ 3119: 洛谷 P5400: 「CTS2019 | CTSC2019」随机立方体

    题目传送门:LOJ #3119. 题意简述: 题目说的很清楚了. 题解: 记恰好有 \(i\) 个极大的数的方案数为 \(\mathrm{cnt}[i]\),则答案为 \(\displaystyle\ ...

  3. LOJ #3119. 「CTS2019 | CTSC2019」随机立方体 组合计数+二项式反演

    好神的一道计数题呀. code: #include <cstdio> #include <algorithm> #include <cstring> #define ...

  4. LOJ #3119「CTS2019 | CTSC2019」随机立方体 (容斥)

    博客链接 里面有个下降幂应该是上升幂 还有个bk的式子省略了k^3 CODE 蛮短的 #include <bits/stdc++.h> using namespace std; const ...

  5. 【LOJ】#3119. 「CTS2019 | CTSC2019」随机立方体

    题解 用容斥,算至少K个极大值的方案数 我们先钦定每一维的K个数出来,然后再算上排列顺序是 \(w_{k} = \binom{n}{k}\binom{m}{k}\binom{l}{k}(k!)^3\) ...

  6. Loj #3124. 「CTS2019 | CTSC2019」氪金手游

    Loj #3124. 「CTS2019 | CTSC2019」氪金手游 题目描述 小刘同学是一个喜欢氪金手游的男孩子. 他最近迷上了一个新游戏,游戏的内容就是不断地抽卡.现在已知: - 卡池里总共有 ...

  7. 「CTS2019 | CTSC2019」氪金手游 解题报告

    「CTS2019 | CTSC2019」氪金手游 降 智 好 题 ... 考场上签到失败了,没想容斥就只打了20分暴力... 考虑一个事情,你抽中一个度为0的点,相当于把这个点删掉了(当然你也只能抽中 ...

  8. LOJ 3120: 洛谷 P5401: 「CTS2019 | CTSC2019」珍珠

    题目传送门:LOJ #3120. 题意简述: 称一个长度为 \(n\),元素取值为 \([1,D]\) 的整数序列是合法的,当且仅当其中能够选出至少 \(m\) 对相同元素(不能重复选出元素). 问合 ...

  9. @loj - 3120@ 「CTS2019 | CTSC2019」珍珠

    目录 @description@ @solution@ @accepted code@ @details@ @description@ 有 \(n\) 个在范围 \([1, D]\) 内的整数均匀随机 ...

随机推荐

  1. SQL Server 中用While循环替代游标Cursor的解决方案

    在编写SQL批处理或存储过程代码的过程中,经常会碰到有些业务逻辑的处理,需要对满足条件的数据记录逐行进行处理,这个时候,大家首先想到的方案大部分是用“游标”进行处理. 举个例子,在订单管理系统中,客服 ...

  2. 程序猿看了要xiao了

    一. 程序猿问科比:“你为什么这么成功? ”科比:“你知道洛杉矶凌晨四点是什么样子吗? ”程序猿:“知道,一般那个时候我还在写代码,怎么了?”科比:“额…….” 二. 女神:你能让这个论坛的人都吵起来 ...

  3. FastDFS介绍(一)

    1.简介 FastDFS对文件进行管理,功能包括:文件存储.文件同步.文件访问(文件上传.文件下载.文件删除)等,解决了大容量文件存储的问题,特别适合以文件为载体的在线服务,如相册网站.文档网站.图片 ...

  4. 检测代理IP匿名程度的方法,很实用

    做网络的基本都知道代理,这个是肯定的,不管是用花刺还是猎手的网页代理,还是直接VPN的通道代理,代理有着不用说大家也知道的重要性.不管是做CPA还是做点击亦或者投票,代理都能帮我们一下,虽然帮的忙不大 ...

  5. .NET Core:目录

    ylbtech-.NET Core:目录 1.返回顶部 1. https://dotnet.microsoft.com/ 2. 2.返回顶部   3.返回顶部   4.返回顶部   5.返回顶部   ...

  6. 使用 Select2 下拉框实现复选

    使用 Select2 下拉框实现复选 <!DOCTYPE html> <html> <head> <meta http-equiv="Content ...

  7. ubuntu搭建jdk+jenkins

    第一步,安装jdk(如果已安装,直接进行第二步) 1.下载 jdk-8u172-linux-x64.tar.gz       点此下载   2.解压 tar -zxvf jdk-8u172-linux ...

  8. curl 中关于 CURLINFO_HEADER_SIZE 的 BUG 定位及修复

    curl 官方下载页面 CentOS7 默认安装的 curl 版本太低了,需要升级为最新版. 1. 问题描述 对接了一个接口,用来下载 PDF 文件.使用 curl 下载后,文件老是报错无法打开.接口 ...

  9. Mac获取Jenkins管理员初始密码

    前言 最近在配置jenkins环境,但是启动jenkins后,进入jenkins解锁页时,需要自己获取初始密码. 尝试在访达中输入地址搜索,结果无该文件.后来百度上查看多篇文章后,终于获取到了初始密码 ...

  10. linux点滴记录

    以下均为在Ubuntu下实践操作 更改DNS //编辑文件 - “/etc/resolv.conf”,打开“终端应用程序”-“附件” - “终端”,在终端里输入下面的命令: sudo nano /et ...