UOJ#551【UNR #4】校园闲逛 【多项式好题（雾）】

题目描述

Link

$n$ 个点 $m$ 条边的无向带权图，$1\le$ 边权 $\le V$，$Q$ 次询问，从 $x$ 到 $y$ 有多少条权值和为 $w$ 的路径，路径可以经过重复边和重复点。

$1\le n\le 8,0\le m\le 3\ast 10^5,1\le V\le 65000,0\le Q\le 10000$

题目分析

暴力 DP：$dp[i][j][v]=\sum _{k=1}^{n}dp[i][k][v^{\prime }]\ast num[k][j][v-v^{\prime }]$

显然是个卷积形式。

• 可以直接分治NTT，看做 $F_v=F_{v^{\prime }}\ast W_{v-v^{\prime }}$，$F$ 是个矩阵，比如 Code1 ，卡着时限过，复杂度 $O(n^{3}V\log V+n^{2}V{\log }^{2}V)$。
  
  

• 可以用一些黑科技，因为分治时若 $l\ne 0$ 则 $2l>r$，所以 $F_{0,r-l}$ 都已经算好了，而此时 $F_{l,r}$ 都是由 $F_{0,l-1}$ 贡献来的，所以用 $F_{l,r}\ast F_{0,r-l}$ 当作新的 $F_{l,r}$，这样不会算重也不会算漏。
  
  所以分治时当 $l\ne 0$ 时就不用继续递归了，直接在这一层处理完。这样复杂度就变成了 $T(V)=T(\frac{V}{2})+n^{3}V+n^{2}V\log V=O(n^{3}V+n^{2}V\log V)$
  比如 Code2 。
  
  

• 可以用一些高科技
  令 $g_{x,y}$ 为 $x\to y$ 路径权值和对应方案数的生成函数，$G$ 为 $g$ 对应的 $n\times n$ 的矩阵，$W$ 为边权多项式对应的矩阵，那么有 $G=GW+I$，于是 $G=\frac{I}{I-W}$。
  显然 $W$ 矩阵里的每个多项式的常数项都不为 1，因此 $I-W$ 有逆，对这个元素为多项式的矩阵求逆，用类似高斯消元的初等变换法，需要 $n$ 次多项式求逆，$O(n^2)$ 次 DFT，$n^3$ 次点值乘法和多项式加法，复杂度 $O(n^{3}V+n^{2}V\log V)$，略有一些(大)常数。

  或者转换一下思路，令 $A(z)$ 表示权值为 $i$ 的路径的矩阵生成函数，即多项式的每项系数都是一个矩阵，$x^i$ 表示权值和为 $i$，矩阵内的元素表示方案数，并令 $W$ 为对应的边权矩阵多项式，那么有 $A=AW+I$，于是 $A=\frac{I}{I-W}$，但此时 $W$ 是个多项式，$I$ 相当于是这个多项式的常数项，要对多项式求逆，用牛顿迭代法，有 $x_1=x_0(2-A{x}_0)$
  实现上有些需要注意的地方，虽然多项式系数为矩阵，但是DFT的时候因为只有数乘，所以可以对矩阵的每个元素分别 DFT。然后每个位置的点值现在是一个矩阵，做乘法的时候需要注意 $A{x}_0$ 的顺序不能变，因为矩阵不满足交换律。最开始的牛顿迭代的式子是 $A{x}_0-I=0$，那么平方之后是 $A{x}_{0}A{x}_0-2A{x}_0+I=0$，即 $A{x}_0(2-A{x}_0)=I$。因为左逆等于右逆，所以式子也可以写为 $x_{0}A-I=0$，推出 $x_1=(2-x_{0}A)x_0$，可以发现两个形式是一样的。
  复杂度同样是 $O(n^{3}V+n^{2}V\log V)$

  比如 Code3 。