ABC216H - Random Robots(容斥,状压DP)
题面
有
K
K
K 个机器人初始分别位于数轴上
x
1
,
x
2
,
.
.
.
,
x
K
x_1,x_2,...,x_{K}
x1,x2,...,xK 的整点位置。
接下来会经历
N
N
N 秒,每一秒都会发生如下事件:
- 每个机器人分别有一半的概率停住不动,有一般的概率往坐标轴正方向移动一单位距离。每个机器人的移动是同时进行的。
问机器人互相不碰撞的概率是多少。对
998244353
998244353
998244353 取模。
2
≤
K
≤
10
,
1
≤
N
≤
1000
,
0
≤
x
1
<
x
2
<
.
.
.
<
x
K
≤
1000
.
2\leq K\leq10,1\leq N\leq1000,0\leq x_1<x_2<...<x_K\leq1000\,.
2≤K≤10,1≤N≤1000,0≤x1<x2<...<xK≤1000.
题解
听说要用
L
G
V
\tt LGV
LGV ?听他吹吧
还是离得挺远的,只是用了下
L
G
V
\tt LGV
LGV 引理的核心——容斥原理罢了。
把问题转化到二维平面,每个机器人在
x
x
x 轴方向照常,在
y
y
y 轴方向每秒都前进一格。看起来是不是很像
L
V
G
\tt LVG
LVG 的问题了?但这样完全没用。
我们把最终的机器人序列从左到右记录为
P
P
P,为了利用容斥原理,我们定义一个不合法条件为某个逆序对是否存在,因为存在某个逆序对,可以等价于两个机器人相撞,而最终答案不受影响。可以用
L
G
V
LGV
LGV 引理的路径交换思想加深理解。
形式化地,我们可以把答案表示为:
∑
P
,
Q
(
−
1
)
σ
(
P
)
⋅
(
1
2
)
N
K
⋅
∏
i
=
1
K
C
(
N
,
Q
i
−
x
i
)
\sum_{P,Q} (-1)^{σ(P)}\cdot(\frac{1}{2})^{NK}\cdot\prod_{i=1}^K {\rm C}(N,Q_i-x_i)
P,Q∑(−1)σ(P)⋅(21)NK⋅i=1∏KC(N,Qi−xi)
序列
Q
Q
Q 是最终的位置序列,是递增的,中间的
(
1
2
)
N
K
(\frac{1}{2})^{NK}
(21)NK 是一个方案的概率,右边用组合数算出了每种
(
P
,
Q
)
(P,Q)
(P,Q) 的情况数。
现在,我们面前有两个艰难的选择:壹,用状态压缩DP来模拟排列
P
P
P 的形成,同时用DP模拟序列
Q
Q
Q 的形成;贰,无比简单地直接枚举排列,然后计算别的,序列
Q
Q
Q 、组合数……
后者,
n
!
n!
n! 直接达到了
3
e
6
3e6
3e6 级别,后面的没法算了。所以我们选择前者。定义
d
p
[
S
]
[
j
]
dp[S][j]
dp[S][j] 表示当前已经往排列中(从左往右)填入的机器人集合为
S
S
S,同时最靠右的机器人最终位置
≤
j
\leq j
≤j 的带权方案数。加入一个机器人,可以通过
S
S
S 算出前面比它标号大的个数,维护逆序对,而且它的最终位置要比前面的都靠右,所以很好维护,我们令
d
(
S
,
i
)
d(S,i)
d(S,i) 表示
S
S
S 中比
i
i
i 标号大的机器人个数,有如下转移:
d
p
[
S
]
[
j
]
=
d
p
[
S
]
[
j
−
1
]
+
∑
i
∈
S
(
−
1
)
d
(
S
,
i
)
d
p
[
∁
S
i
]
[
j
−
1
]
⋅
C
(
N
,
j
−
x
i
)
⋅
(
1
2
)
N
dp[S][j]=dp[S][j-1]~~+~~\sum_{i\in S}(-1)^{d(S,i)}dp[\complement_S\,i][j-1]\cdot {\rm C}(N,j-x_i)\cdot (\frac{1}{2})^{N}
dp[S][j]=dp[S][j−1] + i∈S∑(−1)d(S,i)dp[∁Si][j−1]⋅C(N,j−xi)⋅(21)N
最终答案就是
d
p
[
U
]
[
max
x
+
N
]
dp[U][\max x+N]
dp[U][maxx+N] 取模
998244353
998244353
998244353 ,时间复杂度
O
(
2
K
K
N
)
{\rm O}(2^{K}KN)
O(2KKN) 。
如何评价 jiangly 说用
L
G
V
\tt LGV
LGV ,结果代码打的是状压DP ? 其实说的有道理。
CODE
#include<set>
#include<map>
#include<stack>
#include<cmath>
#include<ctime>
#include<queue>
#include<bitset>
#include<cstdio>
#include<cstring>
#include<iostream>
#include<algorithm>
using namespace std;
#define MAXN 1005
#define LL long long
#define DB double
#define ENDL putchar('\n')
#define lowbit(x) (-(x) & (x))
#define FI first
#define SE second
#define SI(x) multiset<x>::iterator
#define MI map<int,int>::iterator
#define eps (1e-4)
#pragma GCC optimize(2)
LL read() {
LL f=1,x=0;char s = getchar();
while(s < '0' || s > '9') {if(s=='-')f = -f;s = getchar();}
while(s >= '0' && s <= '9') {x=x*10+(s-'0');s = getchar();}
return f*x;
}
void putpos(LL x) {
if(!x) return ;
putpos(x/10); putchar('0'+(x%10));
}
void putnum(LL x) {
if(!x) putchar('0');
else if(x < 0) putchar('-'),putpos(-x);
else putpos(x);
}
void AIput(LL x,char c) {putnum(x);putchar(c);}
const int MOD = 998244353;
int n,m,s,o,k;
int dp[1<<10|5][MAXN<<1];
int fac[MAXN<<1],inv[MAXN<<1],invf[MAXN<<1];
int invp[MAXN<<1];
int C(int n,int m) {
if(m < 0 || m > n) return 0;
return fac[n] *1ll* invf[n-m] % MOD *1ll* invf[m] % MOD;
}
int F(int xx) {
return C(m,xx) *1ll* invp[m] % MOD;
}
int xx[MAXN];
int main() {
n = read(); m = read();
for(int i = 1;i <= n;i ++) {
xx[i] = read() + 1;
}
fac[0]=fac[1]=inv[0]=inv[1]=invf[0]=invf[1]=1;
invp[0] = 1;
for(int i = 2;i <= m+3;i ++) {
fac[i] = fac[i-1] *1ll* i % MOD;
inv[i] = (MOD-inv[MOD%i]) *1ll* (MOD/i) % MOD;
invf[i] = invf[i-1] *1ll* inv[i] % MOD;
}
for(int i = 1;i <= m;i ++) invp[i] = invp[i-1] *1ll* inv[2] % MOD;
dp[0][0] = 1;
for(int j = 1;j <= m+1001;j ++) dp[0][j] = 1;
int tp = (1<<n);
for(int i = 1;i < tp;i ++) {
for(int j = 1;j <= m+1001;j ++) {
int kk = 1;
dp[i][j] = dp[i][j-1];
for(int k = n;k > 0;k --) {
if(i & (1<<(k-1))) {
(dp[i][j] += kk*1ll*dp[i-(1<<(k-1))][j-1] % MOD *1ll* F(j-xx[k]) % MOD) %= MOD;
kk = MOD-kk;
}
}
}
}
printf("%d\n",dp[tp-1][m+1001]);
return 0;
}
ABC216H - Random Robots(容斥,状压DP)的更多相关文章
- bzoj2669[cqoi2012]局部极小值 容斥+状压dp
2669: [cqoi2012]局部极小值 Time Limit: 3 Sec Memory Limit: 128 MBSubmit: 774 Solved: 411[Submit][Status ...
- ARC 093 F Dark Horse 容斥 状压dp 组合计数
LINK:Dark Horse 首先考虑1所在位置. 假设1所在位置在1号点 对于此时剩下的其他点的方案来说. 把1移到另外一个点 对于刚才的所有方案来说 相对位置不变是另外的方案. 可以得到 1在任 ...
- [BZOJ2669][CQOI2012]局部最小值(容斥+状压DP)
发现最多有8个限制位置,可以以此为基础DP和容斥. $f_{i,j}=f_{i-1,j}\times (cnt_j-i+1)+\sum_{k\subset j} f_{i-1,k}$ $cnt_j$表 ...
- bzoj3812 主旋律 容斥+状压 DP
题目传送门 https://lydsy.com/JudgeOnline/problem.php?id=3812 题解 考虑对于图的联通性的 DP 的一般套路:总方案 - 不连通的方案. 那么我们只需要 ...
- Comet OJ - Contest #7 C 临时翻出来的题(容斥+状压)
题意 https://www.cometoj.com/contest/52/problem/C?problem_id=2416 思路 这里提供一种容斥的写法(?好像网上没看到这种写法) 题目要求编号为 ...
- HDU5731 Solid Dominoes Tilings 状压dp+状压容斥
题意:给定n,m的矩阵,就是求稳定的骨牌完美覆盖,也就是相邻的两行或者两列都至少有一个骨牌 分析:第一步: 如果是单单求骨牌完美覆盖,请先去学基础的插头dp(其实也是基础的状压dp)骨牌覆盖 hiho ...
- ARC093F Dark Horse 【容斥,状压dp】
题目链接:gfoj 神仙计数题. 可以转化为求\(p_1,p_2,\ldots,p_{2^n}\),使得\(b_i=\min\limits_{j=2^i+1}^{2^{i+1}}p_j\)都不属于\( ...
- BZOJ 2560: 串珠子 (状压DP+枚举子集补集+容斥)
(Noip提高组及以下),有意者请联系Lydsy2012@163.com,仅限教师及家长用户. 2560: 串珠子 Time Limit: 10 Sec Memory Limit: 128 MB Su ...
- UOJ #129 / BZOJ 4197 / 洛谷 P2150 - [NOI2015]寿司晚宴 (状压dp+数论+容斥)
题面传送门 题意: 你有一个集合 \(S={2,3,\dots,n}\) 你要选择两个集合 \(A\) 和 \(B\),满足: \(A \subseteq S\),\(B \subseteq S\), ...
随机推荐
- linux目录结构及定时任务
1. Linux的根目录(最顶层的目录) windows系统有根目录:c盘的根目录就是c:\ d盘的根目录就是d:\ 每个盘(分区)都有自己的根目录 Linux系统, 也支持多个分区 Linux的分区 ...
- bat-静默安装winrar并设置系统级环境变量
@echo off Setlocal enabledelayedexpansion @REM vscode中自动开启延迟环境变量扩展 echo 安装winrar "winrar v.5.71 ...
- jQuery获取市、区县、乡镇、村
效果图: 首先根据自己方法把地区树状结构json字符串拿到 html下拉框和js写法如下: <select class="form-control" style=" ...
- Python递归函数的定义和几个小例子
递归函数 (1)什么是递归函数? 我们都知道,一个函数可以调用其他函数.如果这个函数在内部调用它自己,那么这个函数就叫递归函数. (2)递归函数的作用 举个例子,我们来计算阶乘 n! = 1 * 2 ...
- Nginx+Keepalived+VIP漂移实现HA高可用技术之详细教程
https://www.cnblogs.com/zcc666/p/13141626.html 这个是nginx安装教程地址 https://www.cnblogs.com/zcc666/p/1313 ...
- Python 数据科学手册:读书笔记概论
为防止遗忘,在空闲时间将读书的笔记开始按照章节进行概括总结(2022.1.1): 第二章:NumPy 入门 第三章:Pandas 数据处理 第四章:Matplotlib 数据可视化 第五章:机器学习 ...
- JavaScript基本知识点——带你逐步解开JS的神秘面纱
JavaScript基本知识点--带你逐步解开JS的神秘面纱 在我们前面的文章中已经深入学了HTML和CSS,在网页设计中我们已经有能力完成一个美观的网页框架 但仅仅是网页框架不足以展现出网页的魅力, ...
- Leetcode 1331. 数组序号转换
给你一个整数数组 arr ,请你将数组中的每个元素替换为它们排序后的序号. 序号代表了一个元素有多大.序号编号的规则如下: 序号从 1 开始编号. 一个元素越大,那么序号越大.如果两个元素相等,那么它 ...
- 【Codeforces1706A】 Another String Minimization Problem
官方标签 贪心.字符串 题目描述 输入 输出 样例输入 6 4 5 1 1 3 1 1 5 2 4 1 1 1 1 1 2 4 1 3 2 7 7 5 4 5 5 5 3 5 样例输出 ABABA B ...
- Operating System_via牛客网
题目 链接:https://ac.nowcoder.com/acm/contest/28537/F 来源:牛客网 时间限制:C/C++ 1秒,其他语言2秒 空间限制:C/C++ 131072K,其他语 ...