反射容斥

反射容斥

一类关于格路计数的容斥手法。

考虑一种卡特兰数通项公式的组合证明。

一种卡特兰数的等价描述是从 (0,0) 走到 (n,n) 的路径数,每次只能走简单格路(每一维的变化量的绝对值的和为 1),不能经过 y=x+1 这条直线的路径方案数。

首先如果没有不经过直线的限制的方案数为 (2nn),现在考虑如果不满足限制会发生什么。

对于一个经过 y=x+1 的路径,找到路径与直线的第一个交点 A,将之前的路径沿着 y=x+1 翻折,那么起点翻折后为 (1,1)

那每一条从 (1,1) 走到 (n,n) 的路径都恰好对应一条从 (0,0) 走到 (n,n) 经过 y=x+1 的路径。

所以有公式:Catlan(n)=(2nn)(2nn+1)

运用这个方法可以很好的运用到计数到达 T(n,m)(保证 nm) 且不经过 y=x+1 的直线的方案数。

方案数为:f(n,m)=(n+mn)(n+mn+1),证明与上面类似。

再推广一波,从 (0,0)(n,m) 不经过直线 y=x+b 的方案数为:(n+mn)(n+mn+b)

继续思考:求从 (0,0)(n,m) 不经过直线 y=x+ly=x+r 的方案数 (l<0<r)

先想:反射容斥能解决的是:可以找到路径第一次和某条直线的交点,然后翻折过去。

定义由 L,R 构成的字符串例如 RLR 表示依次先与 y=x+r 相交再和 y=x+l 相交再和 y=x+r 相交。(如果连续多次与一条直线相交只记录一次)

顺次处理字符串中的字符,然后前面的部分关于字符对应的直线翻折,假设最后起点为 x,y,答案为 (n+mn)(nx+mynx)

使用反射容斥可以计算的是所有包含当前串为子序列的字符串所代表的路径总数。(画图方便理解)

要求的答案为字符串为空集。

f(S) 表示处理 S 串得到的答案(即计算。使用容斥可得到。

ans=f()f(L)f(R)+f(LR)+F(RL)F(LRL)F(RLR)

证明:一个串 LRLRL (R 开头同理),空集有 1 的贡献。能成为子序列的有 S 的前缀以及 S 去掉第一位的所有前缀,那么刚好会做 1 的贡献,第一个串的贡献是 1,1,1,1(第一位开始),第二个串的贡献是 1,1,1(第二位开始),两柿相加和恰好为 1,再加上空集的贡献恰好为 0。所以只有当 S 为空集才会贡献。

化简可得:

ans=kZ(n+mnk(rl))(n+mnk(rl)+r)

注意:k 可以取到所有整数,但是有值的位置只有 O(nrl) 个。

【集训队作业2018】count

题目要求统计长度为 n,值域在 [1,m] 且每种值至少出现一次的不同构的笛卡尔树个数。(如果有多个最大值取编号最小的)

首先 n>m 必然无解。

那么对于每个点,左儿子的取值要严格小于它,右儿子的取值小于等于它。

f(rt) 为括号序列, f(rt)=(f(lson))+f(rson)

定义向左链长即为 f(rt) 的最大嵌套层数。

观察到向左链长只要小于等于 m 为笛卡尔树合法的充要条件。

必要性:如果向左链长大于了 m,那么至少有 m+1 个不同的数,不合法。

充分性:考虑对每个节点分配合法的值域范围 [1,x],当前节点直接取 x,左儿子分配 [1,x1],右儿子分配 [1,x] 即可。

注意到一个括号序列恰好对应一棵二叉树,转为对括号序列的计数。

将左括号看成 1,右括号看成 1,转化为计数从 (0,0) 走到 (2n,0),每次只能向右上/右下走,且不能经过 y=1y=m+1 的路径方案数。

发现这与我们介绍的双边界反射容斥不大相同。

考虑一种常见的手法,通过对坐标轴的旋转/对称/平移将已知推未知。

将坐标轴旋转 45 度,变为计数从 (0,0) 走到 (n,n),每次只能向左或者向上走,且不经过 y=x+1y=xm1 的方案数。

套用公式可得:

ans=kZ(n+nnk(2+m))(n+nnk(2+m)+1)

P3266 [JLOI2015] 骗我呢

首先注意到每行 m 个数且互不相同,一共有 m+1 个数,那么中间只会有一个数没被用。

fi,j 表示第 i 行不用值 j 的方案数。又因为 xi,j<xi1,j+1

fi,j=k=0j+1fi1,k=fi1,j+1+k=0jfi1,k=fi1,j+1+fi,j1

每个点是由它右下角的点以及左侧的点转移而来,这并不符合我们的需求。

考虑将底 i 行的所有点集体向右平移 i 个为位置。那么有:

fi,j=fi1,j+fi,j1

这是我们所会的,即从 (0,0) 走到 (n+m+1,n) 不经过 y=x+1y=xm2 的方案数,带入柿子即可。(我知道你想要最终的柿子)

ans=kZ(n+n+m+1n+m+1k(3+m))(n+n+m+1n+m+1k(3+m)+1)

posted @   Hanghang007  阅读(909)  评论(2编辑  收藏  举报
相关博文:
阅读排行:
· 阿里最新开源QwQ-32B,效果媲美deepseek-r1满血版,部署成本又又又降低了!
· Manus重磅发布:全球首款通用AI代理技术深度解析与实战指南
· 开源Multi-agent AI智能体框架aevatar.ai,欢迎大家贡献代码
· 被坑几百块钱后,我竟然真的恢复了删除的微信聊天记录!
· AI技术革命,工作效率10个最佳AI工具
点击右上角即可分享
微信分享提示