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