首页 > 其他分享 >LGV引理

LGV引理

时间:2024-04-26 09:04:56浏览次数:17  
标签:frac 路径 vmatrix cdots 行列式 引理 LGV vdots

LGV引理

行列式引出来的有趣的东西,是与图论的交界处。

LGV 引理大致内容为:

对于一张有向无环图,每条边上都有一个权值 \(w(e)\),记 \(weight(P)\) 表示路径 \(P\) 上所有边的权值乘积,对于一个起点组成的集合 \(A\) 和终点组成的集合 \(B\),满足 \(|A|=|B|\),记 \(e(i,j)\) 表示所有 \(A_i \to B_j\) 的路径的 \(weight\) 之和,那么对于矩阵:

\[M=\begin{bmatrix} e(1,1) & e(1,2) & \cdots & e(1,|B|) \\ e(2,1) & e(2,2) & \cdots & e(2,|B|) \\ \vdots & \vdots & \ddots & \vdots \\ e(|A|,1) & e(|A|,2) & \cdots & e(|A|,|B|) \end{bmatrix} \]

设 \(n=|A|\),引理说明:

\[\det(M)=\sum_{S:A\to B}(-1)^{\tau(p)}\prod_{i=1}^{n}weight(P_i)[P_i:A_i\to B_{p_i}] \]

\(S\) 表示不交路径组,\(p_i\) 表示第 \(i\) 个起点对应的终点。

证明

考虑将行列式展开得到:

\[\begin{aligned} \det(M)&=\sum_p (-1)^{\tau (p)}\prod_{i=1}^{n}e(i,p_i)\\ &=\sum_{p} (-1)^{\tau (p)}\prod_{i=1}^n\sum_{P_i:A_i\to B_{p_i}}weight(P_i) \end{aligned} \]

考虑注意到 \({P_1,P_2,\cdots, P_n}\) 是路径组,考虑枚举路径组。

\[\det(M)=\sum_{T:A\to B}(-1)^{\tau(p)}\prod_{i=1}^{n}weight(T_i)[T_i:A_i\to B_{p_i}] \]

此时的 \(p_i\) 表示为第 \(i\) 个起点的终点为 \(p_i\)。

此时和引理唯一的区别在与 \(T\) 可以是任意路径组。

那么只需要证明相交路径组没贡献为 \(0\)。

那我们只需要对于每个路径组 \(T\) 找到唯一对应的 \(Q\),使得两者相加贡献为 \(0\) 即可。

构造方法为:

找到最小的 \(i\) 使得 \(P_i\) 与别的路径有交点,在选择 \(j\) 使得交点离起点最近,若有多个选最小的 \(j\)。

然后把 \(P_i\) 和 \(P_j\) 交换首个交点到末尾的路径。

那么此时逆序对个数奇偶性改变,后面的乘积不变。

注意:

直接使用 \(\text{LGV}\) 计算的答案只有在平面有向无环图中的意义为不交路径组的方案数。

因为注意到 \(\text{LGV}\) 处理的是点不交路径组,但是非平面图可能边不交,因此这种不合法的方案会被记录。

此时的意义为存在偶数个交点的路径组减去奇数个交点的路径数。

P6657 【模板】LGV 引理

把 \(\text{LGV}\) 矩阵的第 \((i,j)\) 的权值设为 \(\binom{b_j-a_i+n-1}{n-1}\)。

直接跑行列式即可。

Turtles

等价于从 \((1,2)\) 走到 \((n-1,m)\),\((2,1)\) 走到 \((n,m-1)\) 的路径不交方案数,求出方案数后手动行列式即可。

P7736 [NOI2021] 路径交点

对于每一层可以很直接求邻接矩阵的行列式。

就是注意所说的问题,直接暴力求出每个起点到每个终点的路径方案数再套上行列式即可。

[ABC216H] Random Robots

期望转计数,如果这道题只有终点的话那么可以直接使用 \(\text{LGV}\) 定理,其中矩阵中的点 \((i,j)\) 的权值为 \(\binom{m}{B_j-A_i}\),\(A\) 为起点,\(B\) 为终点。

考虑拉普拉斯展开行列式,每次新加入一个一行一列。

具体而言:

\[M=\begin{bmatrix} e(A_1,B_1) & e(A_1,B_2) & \cdots & e(A_1,B_{|S|}) \\ e(A_2,B_1) & e(A_2,B_2) & \cdots & e(A_2,B_{|S|}) \\ \vdots & \vdots & \ddots & \vdots \\ e(A_{|S|},B_1) & e(A_{|S|},B_2) & \cdots & e(A_{|S|},B_{|S|}) \end{bmatrix} \]

设 \(f(S,mx)\) 表示当前考虑了集合 \(S\),最大值小于等于 \(mx\) 的答案。

首先如果 \(B_{|S|}<mx\) 直接有 \(f(S,mx-1)\to f(S,mx)\)。

否则直接对最后一列做拉普拉斯展开:

\[\begin{aligned} \det(M)&=\sum_{i=1}^{|S|}c_{i,|S|}(-1)^{i+|S|}A_{i,|S|} \end{aligned} \]

\(A_{i,j}\) 为余子式。\(c_{i,|S|}\) 表示矩阵中点的权值,那么在本题中:

\[f(S,mx)=\sum_{i=1}^{|S|}(-1)^{i+|S|}\binom{m}{B_{|S|}-A_i}f(S\backslash i,mx-1) \]

复杂度为 \(O(nm 2^n)\)

神秘搬的题目

给定一个 \(n\times n\) 的矩阵主对角线上点的权值,求有多少种填完剩下位置的方案,使得每行每列的数都单调不降。

\(n\le 5000\)

首先注意到左下角填数和右上角填数是对称的,不妨只考虑左下角。

一共只有对行和列分别的限制,我们观察一个维度限制一个维度。

我们对每一行观察,将每一行 \((j,a_{i,j})\) 看成一个坐标,那么每一行对应一条路径 \(P_i\),即从 \((0,0)\) 走到 \((i,a_{i,i})\) 的方案数,每次能从 \((x,y)\) 走到 \((x+1,\ge y)\)。

在考察对列的限制,即对于每一行的路径,设 \(d_{i,j}\) 表示路径 \(P_i\) 上横坐标为 \(j\) 的点的纵坐标,如果不存在则忽略,那么 \(\forall j,d_{i,j}\le d_{i+1,j}\)。

考虑小于等于并不优美,将第 \(i\) 行所有值加上 \(i-1\),不妨设 \(b_i=a_i+i-1\),那么 \(\forall j,d_{i,j}<d_{i+1,j}\),其实换句话说,就是 \(n\) 条路径,每条路径从 \((0,0)\) 走到 \((i,b_i)\),每次只能从 \((x,y)\) 走到 \((x+1,\ge y)\) ,求所有路径都不相交的方案数。

路径不相交自然想到 \(\text{LGV}\) 引理,但是 \(\text{LGV}\) 引理只能计数平面图 \(\text{DAG}\),而这个图并非平面图,考虑转化为网格图计算。

转化一下对 \(d\) 的限制,把每条路径看成一个函数,设 \(g_{i,j}\) 表示函数 \(i\) 纵坐标取到 \(j\) 的时候对应的横坐标。

那么有 \(\forall j,g_{i,j}>g_{i+1,j}\), 那么考虑如下转化:

\(n\) 个起点,第 \(i\) 个起点位于 \((i-1,0)\),第 \(i\) 个终点位于 \((0,a_i+i-1)\),每次只能从 \((x,y)\) 走到 \((x-1,y)\) 或者 \((x,y+1)\),求有多少种不交路径数。

怎么对应原图的呢?考虑当前是第 \(i\) 次走水平步 (向左走),此时纵坐标为 \(j\),那么对于对应原图的坐标 \((i,j)\),最后一步对应 \((i,b_i)\)。

好,现在是我们熟知的网格图了,套上模板可知矩阵上 \((i,j)\) 为 \(\binom{b_j+i-1}{i-1}\),答案为行列式的值。

\[\begin{vmatrix} \binom{b_1+0}{0} & \binom{b_2+0}{0} & \cdots & \binom{b_n+0}{0} \\ \binom{b_1+1}{1} & \binom{b_2+1}{1} & \cdots & \binom{b_n+2}{2} \\ \vdots & \vdots & \ddots & \vdots \\ \binom{b_1+n-1}{n-1} & \binom{b_2+n-1}{n-1} & \cdots &\binom{b_n+n-1}{n-1} \end{vmatrix} \]

打开组合数:

\[\begin{vmatrix} \frac{(b_1+0)!}{b_1!0!} & \frac{(b_2+0)!}{b_2!0!} & \cdots &\frac{(b_n+0)!}{b_n!0!} \\ \frac{(b_1+1)!}{b_1!1!} & \frac{(b_2+1)!}{b_2!1!} & \cdots &\frac{(b_n+1)!}{b_n!1!} \\ \vdots & \vdots & \ddots & \vdots \\ \frac{(b_1+n-1)!}{b_1!(n-1)!} & \frac{(b_2+n-1)!}{b_2!(n-1)!} & \cdots &\frac{(b_n+n-1)!}{b_n!(n-1)!} \end{vmatrix} \]

把每一行的 \((i-1)!\) 提出来

\[\prod_{i=1}^{n}\dfrac{1}{(i-1)!}\begin{vmatrix} \frac{(b_1+0)!}{b_1} & \frac{(b_2+0)!}{b_2!} & \cdots &\frac{(b_n+0)!}{b_n!} \\ \frac{(b_1+1)!}{b_1} & \frac{(b_2+1)!}{b_2!} & \cdots &\frac{(b_n+1)!}{b_n!} \\ \vdots & \vdots & \ddots & \vdots \\ \frac{(b_1+n-1)!}{b_1!} & \frac{(b_2+n-1)!}{b_2!} & \cdots &\frac{(b_n+n-1)!}{b_n!} \end{vmatrix} \]

接下来你需要一些线性代数的基本功:初等矩阵变幻去化简行列式,其实你会发现,行列式的值等价于范德蒙德行列式。

拿 \(3\times 3\) 的矩阵而言

\[\begin{vmatrix} 1 & 1 & 1 \\ (b_1+1) & (b_2+1) & (b_3+1) \\ (b_1+1)(b_1+2) & (b_2+1)(b_2+2) & (b_3+1)(b_3+2) \end{vmatrix} \]

先将第三行减去第二行的两倍,再用第二行减去第一行的一倍可以得到:

\[\begin{vmatrix} 1 & 1 & 1 \\ b_1 & b_2 & b_3 \\ (b_1+1)b_1 & (b_2+1)b_2 & (b_3+1)b_3 \end{vmatrix} \]

再将第三行减去第二行的第一行得到:

\[\begin{vmatrix} 1 & 1 & 1 \\ b_1 & b_2 & b_3 \\ b_1^2 & b_2^2 & b_3^2 \end{vmatrix} \]

这个手法可以很好的推广到高维,具体地式子不再展开说明。

现在转化为求范德蒙德行列式的值,即:

\[\begin{vmatrix} 1 & 1 & \cdots & 1\\ b_1 & b_2 & \cdots &b_n \\ \vdots & \vdots & \ddots &\vdots\\ b_1^{n-1} & b_2^{n-1} & \cdots &b_n^{n-1} \end{vmatrix} \]

考虑将除了第一行的每一行同时减去上一行的 \(b_1\) 得到

\[\begin{vmatrix} 1 & 1 & \cdots & 1\\ 0 & b_2-b_1 & \cdots &b_n-b_1 \\ \vdots & \vdots & \ddots &\vdots\\ 0 & b_2^{n-2}(b_2-b_1) & \cdots &b_n^{n-2}(b_n-b_1) \end{vmatrix} \]

注意到此时第一列只有 \((1,1)\) 非零,那么对第一列行列式展开那么得到值其实就为除去第一行第一列的行列式值。

此时除去第一行第一列之后,把每一列 \((b_i-b_1)\) 提取出来再不断递归即可。

那么得到:

\[\begin{vmatrix} 1 & 1 & \cdots & 1\\ b_1 & b_2 & \cdots &b_n \\ \vdots & \vdots & \ddots &\vdots\\ b_1^{n-1} & b_2^{n-1} & \cdots &b_n^{n-1} \end{vmatrix} =\prod_{1\le i<j \le n}(b_j-b_i) \]

那么最终的答案也就为(右上角的方案和左下角的是对称的,方案数相同)

\[ans=\left( \prod_{i=1}^{n}\dfrac{1}{(i-1)!}\prod_{1\le i<j\le n}(b_j-b_i)\right)^2 \]

直接计算即可得到通过。

后记:使用一些高端的多项式科技可以加速计算。

P9041 [PA2021] Fiolki 2

现在考虑如果一个区间满足区间的最大不交路径组大小恰好为 \(i\),那么说明转到 \(\text{LGV}\) 矩阵上之后存在恰好 \(i\) 行 \(i\) 列的子式的行列式值不为 \(0\),换句话说,\(\text{LGV}\) 矩阵的秩大小为 \(i\)。

考虑现在对每条边赋予一个非零权值,利用拓扑序求出每个起点到当前点的权值乘积。(这里本来应该对每条边设一个形式元 \(x\),然后准确表示,但是维度成本太高,随机一个权值使得秩变小的概率为 \(\dfrac{1}{p}\),证明同 P10102)

设 \(f_i,j\) 表示第 \(i\) 个点从 \(j\) 走到这的权值乘积。

现在只需对每个区间求出 \(f_{i,j}\) 形成的矩阵的秩,从而得到最大不交路径大小。

考虑扫描线,维护前缀线性基(不会的看幸运数字),将线性基中的编号 \(id_i\) 全部拉出来从大到小排序,那么左端点在 \([id_{i+1},id_i)\) 的秩的大小为 \(i\)。

复杂度 \(O(nk^2+mk)\)。

标签:frac,路径,vmatrix,cdots,行列式,引理,LGV,vdots
From: https://www.cnblogs.com/Hanghang007/p/18159160

相关文章

  • LGV引理
    LGV引理行列式引出来的有趣的东西,是与图论的交界处。LGV引理大致内容为:对于一张有向无环图,每条边上都有一个权值\(w(e)\),记\(weight(P)\)表示路径\(P\)上所有边的权值乘积,对于一个起点组成的集合\(A\)和终点组成的集合\(B\),满足\(|A|=|B|\),记\(e(i,j)\)表示所有\(......
  • 置换群 / Polya 原理 / Burnside 引理 学习笔记
    置换群/Polya原理/Burnside引理学习笔记在GJOI上做手链强化,经过长达三小时的OEIS和手推无果后开摆,喜提rnk12,故开始学习置换群相关内容。笔记主要以Polya原理和Burnside引理的应用为主,所以会非常简单,很大一部分的群论概念和证明不会写,因为我不会。基础群论定......
  • LGC引理
    LGV引理内容(不会证明)\(\omega(P)\)表示P这条路径上所有边的边权之积。(路径计数时,可以将边权都设为1)(事实上,边权可以为生成函数)e(u,v)表示u到v的每一条路径P的$\omega(P)$之和,即$e(u,v)=\sum\limits_{P:u\rightarrowv}\omega(P)$。起点集合\(A\),是有向无......
  • LGV引理学习笔记
    ReferenceOI-wiki介绍LGV引理(Lindström–Gessel–Viennotlemma)用来解决有向无环图上不相交路径计数,注意仅适用于有向无环图。给定\(n\)个起点构成的集合\(S\)和\(n\)个终点构成的集合\(T\),定义\(\omega(P)\)表示路径\(P\)上所有边权的乘积(计数时设边权为\(1\)......
  • Burnside 引理 与 Pólya 定理 学习笔记
    为了防止明天就把好不容易听完的东西都还给rabbit_lb了,还是记一点吧。1.群论基础1.1群(group)的定义给定集合\(G\)和\(G\)上的二元运算\(\cdot\),满足下列条件称之为群:封闭性:若\(a,b\inG\),则\(a\cdotb\inG\)。结合律:对于任意\(a,b,c\inG\),有\((a\cdotb)\cd......
  • LGV 引理
    这个东西一般用来求DAG中从初始点集合\(S\)到终止点集合\(T\)的有符号不相交路径方案数(不相交指的是点不会同时出现在两个路径中),\(n=|S|=|T|\)。设\(P(w)\)表示路径\(w\)上的边权的乘积,\(e(s,t)\)表示\(\sum_{w:s\tot}P(w)\)。\[A=\begin{bmatrix}&e(A_1,B_1)&e......
  • 【数学】LGV 引理
    题目描述这是一道模板题。有一个\(n\timesn\)的棋盘,左下角为\((1,1)\),右上角为\((n,n)\),若一个棋子在点\((x,y)\),那么走一步只能走到\((x+1,y)\)或\((x,y+1)\)。现在有\(m\)个棋子,第\(i\)个棋子一开始放在\((a_i,1)\),最终要走到\((b_i,n)\)。问有多少种方案,使得......
  • 对acwing355异象石引理的证明
    首先我们抽象一下这道题的模型,然后把引理记住模型:对于一棵树上选定的一些点,把他们连通起来的最小边数我们先考虑一种朴素做法,对于任何一种方案,任取其中两个点,那么这个方案一定包含这两个点之间的路径就是说,我们依次添加每个点,对于每一个新添加进来的点,让这个点与其已经添加的点......
  • 周期引理的代数证明
    翻译自https://zhuanlan.zhihu.com/p/85169630字符串是0-index.周期引理:对于长为\(n\)的字符串\(s\),如果\(p,q\)均为\(s\)的周期,并且\(p+q-\gcd(p,q)\leqn\),那么\(\gcd(p,q)\)也是\(s\)的周期。定义\(s_p(i)=s_{i\bmodp},s_q(i)=s_{i\bmodp}\),用生成函数来描......
  • 【模板】图的计数相关:行列式及求值、Matrix-Tree 定理、BEST 定理、LGV 引理
    归类为线性代数、图论。证明都是神仙,特别是名字带“理”的,不证了。行列式定义行列式(Determinant)是对\(n\)阶方阵\(A\)定义的,是一个标量。\(A\)的\(n\)阶行列式\(det(A)\)或\(|A|\)定义如下:\[det(A)=\sum_p(-1)^{\mu(p)}\prod_{i}A[i][p_i].\]这里将排列的逆序对数......