(其实是贺的:https://www.luogu.com.cn/paste/whl2joo4

LGV 引理

LGV 引理,即 Lindström–Gessel–Viennot lemma .

一个带权 DAG \(G\) 中有起点集 \(A=\{a_1, a_2, \cdots, a_n\}\),终点集 \(B=\{b_1, b_2, \cdots, b_n\}\) .

不相交路径定义为:一组 \(u\to v\) 的不相交路径 \(P=(P_1,P_2,\cdots,P_n)\),满足 \(P_i\) 是一条 \(a_i\) 到 \(b_{\sigma(P)_i}\) 的路径,其中 \(\sigma(P)\) 是 \(P\) 对应的排列 .

且对于任意 \(i\neq j\),\(P_i, P_j\) 无公共点(即不相交)

令 \(\omega(P)\) 为路径 \(P\) 上所有边权的积 .

对顶点 \(u,v\),定义

\[e(u,v) = \sum_{P:a \to b}{\omega(P)}
\]

令 \(\tau(p)\) 表示排列 \(p\) 的逆序对个数 .


令矩阵

\[M=\begin{bmatrix}e(a_1, b_1) & e(a_1, b_2) & \cdots & e(a_1, b_n) \\ e(a_2, b_1) & e(a_2, b_2) & \cdots & e(a_2, b_n) \\ \vdots & \vdots & \ddots & \vdots \\ e(a_n, b_1) & e(a_n, b_2) & \cdots & e(a_n, b_n)\\\end{bmatrix}
\]

则有

\[\boxed{\det(M)=\sum_{P:A\to B}{(-1)^{\tau(P)}\prod_{i=1}^{n}\omega(P_i)}}
\]

证明:

展开行列式的定义,得 LGV 引理等价于

\[\sum_p (-1)^{\tau(p)}\prod_{i=1}^ne(a_i, b_{p(i)})=\sum_{P:A\to B}{(-1)^{\tau(P)}\prod_{i=1}^{n}\omega(P_i)}
\]

其中 \(p\) 是一个 \(1\) 到 \(n\) 的排列 .

妈呀左右两边长得这么像!

观察组合意义!我们可以发现:

  • LHS 其实就是选 \(n\) 条起点和终点分别互不相同的路径,然后把他们的带符号和加起来 .
  • RHS 几乎一样,只不过 \(\displaystyle \sum_p\) 变成了 \(\displaystyle \sum_{P:A\to B}\),也就是要求路径不相交 .

故我们考虑构造双射 .

如果两个人的路径相交了,我们可以将相交后的部分取反,就相当于这两人分别走到了对方的终点去 . 这样是不是就不交了!!

是不是很平凡! 证明真牛逼 orz .


AJH 大佬说边权只要是交换环就行,也就是说可以是 GF .

太强了 orz

不相交路径计数

LGV 引理的直接应用 .

有一个图 \(G\) 是 平面图 且是 DAG .

从每个 \(a_i\) 到 \(b_i\),在满足路径不交的前提下,所有方案中路径边权乘积之和 .

是不是就是板子啊 .

你用一种方法算出 \(e\),然后是不是构造出 \(M\) 高斯消元求行列式就完了 .


为什么要是平面图?


一般的 LGV 题符合要求的匹配只有一组(e.g. 网格图),因此按照匹配的顺序列矩阵求出的就是不交路径数量 .

然而如果不是平面图,直接求 \(\det(M)\) 其实并不一定是不交路径数量 .

具体可以看 NOI2021 路径交点 .

例题

Luogu6657. 【模板】LGV 引理

link: https://www.luogu.com.cn/problem/P6657


一个 \(n\times n\),每步只能往右下走 .

\(m\) 个棋子,初始在 \((a_i,1)\),要到 \((b_i,n)\) .

求不交路径方案数,对 \(998244353\) 取模 .

不相交路径计数 弱化版 .

显然对于任何一个路径,\(\omega(P)=1\),这样我们才能计数嘛 .

因为是网格图,所以拿出我们老生常谈的网格图两点路径计数:

\[e(u,v)=\dbinom{v-u+n-1}{n-1}
\]

(因为从 \((x_1,y_1)\) 到 \((x_2,y_2)\) 的方案数是 \(\dbinom{x_2-x_1+y_2-y_1}{x_2-x_1}\))

然后预处理阶乘及其逆元算 binom,然后高斯消元求行列式就完了 .

时间复杂度 \(O(n+m^3)\) .

CF348D Turtle

link: https://codeforces.com/contest/348/problem/D


网格图去掉一些点,两个点从 \((1,1)\) 走到 \((n,m)\),求不交路径数取模 \(10^9+7\) .

大力 LGV 引理 .

把两个相同起始点拆开,起点拆成 \((1,2)\),\((2,1)\),终点拆成 \((n-1,m)\),\((n,m-1)\) .

然后 LGV 引理可得一个二阶行列式,对角线法则乘开即可 .

然后变成带限制路径计数,做四边 DP 即可 .

时间复杂度 \(O(nm)\) .

Monotonic Matrix

link: https://www.nowcoder.com/acm/contest/139/A

求满足以下条件 \(n\times m\) 矩阵的数量:

  • \(A_{i,j}\in\{0, 1, 2\}\) .
  • \(A_{i,j}\le A_{i,j+1}\),\(A_{i,j}\le A_{i+1,j}\) .

对 \(10^9+7\) 取模 .

从起点 \((0,0)\) 终点 \((n,m)\) 画一条非降路径,然后 \(0,1\) 是沿着网格走的过程 .

从起点 \((-1,1)\) 终点 \((n-1,m+1)\) 画一条非降路径,然后 \(1,2\) 是沿着网格走的过程 .

然后用 LGV 引理算出不交的路径的方案数即可 .

然而这玩意也是二阶行列式,就是几个组合数乘一下减一下 .

时间复杂度在于求组合数 .

习题

Reference

LGV 引理的更多相关文章

  1. 2021牛客暑期多校训练营9C-Cells【LGV引理,范德蒙德行列式】

    正题 题目链接:https://ac.nowcoder.com/acm/contest/11260/C 题目大意 一个平面上,\(n\)个起点\((0,a_i)\)分别对应终点\((i,0)\),每次 ...

  2. P7736-[NOI2021]路径交点【LGV引理】

    正题 题目链接:https://www.luogu.com.cn/problem/P7736 题目大意 有\(k\)层的图,第\(i\)层有\(n_i\)个点,每层的点从上到下排列,层从左到右排列.再 ...

  3. LGV 引理小记

    讲个笑话,NOI 之前某场模拟赛让我知道了这个神奇的科技,于是准备 NOI 之前学完,结果鸽着鸽着就鸽掉了,考 day1 之前一天本来准备花一天时间学的,然后我就开玩笑般地跟自己说,这么 trivia ...

  4. LGV 引理——二维DAG上 n 点对不相交路径方案数

    文章目录 引入 简介 定义 引理 证明 例题 释疑 扩展 引入 有这样一个问题: 甲和乙在一张网格图上,初始位置 ( x 1 , y 1 ) , ( x 2 , y 2 ) (x_1,y_1),(x_ ...

  5. P6657-[模板]LGV 引理

    正题 题目链接:https://www.luogu.com.cn/problem/P6657 题目大意 给出$n\times n$的棋盘,$m$个起点第$i$个为$(1,a_i)$,对应$m$个终点第 ...

  6. ZROI 暑期高端峰会 A班 Day1 组合计数

    AGC036F Square Constriants 一定有 \(l_i<p_i\le r_i\). 考虑朴素容斥,枚举每个数是 \(\le l_i\) 还是 \(\le r_i\).对于 \( ...

  7. ACM模板_axiomofchoice

    目录 语法 c++ java 动态规划 多重背包 最长不下降子序列 计算几何 向量(结构体) 平面集合基本操作 二维凸包 旋转卡壳 最大空矩形 | 扫描法 平面最近点对 | 分治 最小圆覆盖 | 随机 ...

  8. 【置顶】Trotyl's OI tree

    \(\rm thx\):@\(\rm UntilMadow\) ! \(\color{Green}{\rm Pupil}\) :只会一点点 \(\color{blue}{\text{Expert}}\ ...

  9. NOI2021游记

    NOI2021游记 前言 写于 2021.7.28,成绩榜刚出后几个小时.总分 345 拿到银牌 183 名. 我的高中 OI 生活在这里画上句号.结局对我而言虽然不够完美,但是无论怎样都是我人生道路 ...

随机推荐

  1. mysql二进制日志和mysql备份工具介绍以及日志恢复

    mysql备份: 三种备份方式 冷备:数据库停机,在进行备份 热备:lock table锁表,read   数据库只可以读不能写,在备份 温备:备份时数据库正常运行 备份类型:完整备份:全部备份,部分 ...

  2. 满满干货!手把手教你实现基于eTS的分布式计算器

    最近收到很多小伙伴反馈,想基于扩展的TS语言(eTS)进行HarmonyOS应用开发,但是不知道代码该从何处写起,从0到1的过程让新手们抓狂. 本期我们将带来"分布式计算器"的开发 ...

  3. boost::bind 不能处理函数重载 (error: no matching function for call to 'bind')

    前言 最近任务多.工期紧,没有时间更新博客,就水一期吧.虽然是水,也不能太失水准,刚好最近工作中遇到一个 boost::bind 的问题,花费了半天时间来定位解决,就说说它吧. 问题背景 项目中使用了 ...

  4. 《Mybatis 手撸专栏》第10章:使用策略模式,调用参数处理器

    作者:小傅哥 博客:https://bugstack.cn 沉淀.分享.成长,让自己和他人都能有所收获! 一.前言 你这代码写的,咋这么轴呢! 说到轴,让我想起初中上学时老师说的话:"你那脑 ...

  5. 架构师必备:HBase行键设计与应用

    首先要回答一个问题,为何要使用HBase? 随着业务不断发展.数据量不断增大,MySQL数据库存在这些问题: MySQL支持的数据量为TB级,不能一直保留历史数据.而HBase支持的数据量为PB级,适 ...

  6. JavaScript数据类型BigInt实践之id数值太大,导致前后端交互异常

    项目开发中前后端数据交互常会使用id作为主键索引,通常id数值都不大,使用number类型就可以表示处理,但对于一些分布式id或其他情况,id数值太大且超过了JS的最大处理数(Math.pow(2,  ...

  7. AsList()方法详解

    AsList()方法详解 在Java中,我们应该如何将一个数组array转换为一个List列表并赋初始值?首先想到的肯定是利用List自带的add()方法,先new一个List对象,再使用add()方 ...

  8. [自制操作系统] 第04回 完善MBR

    目录 一.前景回顾 二.改写MBR 三.实现loader 一.前景回顾 在之前我们说到,MBR的作用便是加载操作系统内核到指定位置.而MBR需要通过读取硬盘来获得操作系统内核.在上一回我们已经讲解了硬 ...

  9. 打字练习-编程语言关键字系列-html

    以下是小编整理的部分html关键字,专门给有需要的朋友进行打字练习用,通过打字练习的方式,既提高了打字速度,又可以熟悉html关键字~~~ www, url, http, W3C, html, htm ...

  10. 整理orcal常用sql语句

    1.表插入列 alter table XMJ_ONE add column1 NUMBER(38) default 0;comment on column XMJ_ONE.column1 is '字段 ...