[BZOJ5481]矩阵

题目

传送门 to BZOJ（可是已经凉透了）

题目概要
对于每个 $i\in [2,n]$，求出有多少个 $i\times i$ 的 $01$ 矩阵满足：每一行每一列的和都是 $2$ 。

数据范围与约定
$n⩽10^7$ 。

思路

将第 $i$ 行都看做节点 $r_i$，第 $i$ 列都看做节点 $c_i$，若 $a_{ij}=1$ 则将 $r_i$ 和 $c_i$ 之间连一条无向边。可以发现我们就是做 二分图环剖分。事实上 $\mathrm{o}\mathrm{e}\mathrm{i}\mathrm{s}$ 已收录此数列。

无论如何，需要先求出 $g_n$ 为，大小为 $2n$ 的环有多少种。显然 $\mathfrak{X}$ 部和 $\mathfrak{Y}$ 部可以分别排列，避免重复而规定 $\mathfrak{X}$ 的某个元素为首，又考虑到环是不分方向的，可知
$$g_n=\frac{n!(n-1)!}{2}(n⩾2)$$

思路壹

枚举其中的一个环，要求必须包含某个特定点。记答案为 $f_n$ 则
$$\begin{array}{rl}{f}_n& =\sum _{x=2}^n\left(\genfrac{}{}{0ex}{}{n}{x}\right)\left(\genfrac{}{}{0ex}{}{n-1}{x-1}\right)g_x{f}_{n-x}\\ & =\frac{n!(n-1)!}{2}\sum _{x=2}^n\frac{x!(x-1)!}{x![(n-x)!{]}^2(x-1)!}{f}_{n-x}\\ \Rightarrow \frac{f_n}{(n!{)}^2}& =\frac{1}{2n}\sum _{x=2}^n\frac{f_{n-x}}{[(n-x)!{]}^2}\end{array}$$

简单前缀和可做到 $\mathcal{O}(n)$ 。

思路贰

有标号连通块划分？记 $F(x)={\sum }_{n⩾2}\frac{n!(n-1)!}{2}\cdot \frac{x^n}{(n!{)}^2}={\sum }_{n⩾2}\frac{x^n}{2n}$，则只需求出
$$G(x)=\exp [F(x)]$$

注意这里的 $F(x)$ 不是常见的指数型生成函数的形式；因为它的组合数系数实际上是 ${\left(\genfrac{}{}{0ex}{}{n}{x}\right)}^2$，所以阶乘需平方。

这是可以求解的。只需注意到
$$F^{\prime }(x)=\sum _{n⩾2}\frac{x^{n-1}}{2}=\frac{x}{2-2x}$$

而后有
$$G^{\prime }(x)=G(x)F^{\prime }(x)\Rightarrow (2-2x)\cdot G^{\prime }(x)=xG(x)$$

设 $G(x)={\sum }_{n⩾0}{f}_n{x}^n$，对比系数有
$$2(n+1)f_{n+1}-2n{f}_n=f_{n-1}$$

整理后可以 $\mathcal{O}(n)$ 递推。

代码

此处给出 思路壹 的代码。

#include <cstdio>
const int mod = 998244353;
const int MAXN = 10000005;
int n, inv[MAXN];
int main () {
	scanf("%d", &n);
	inv[0] = inv[1] = 1;
	for(int i = 2; i <= n; ++i) // 逆元
		inv[i] = 1ll * (mod - mod/i) * inv[mod%i] % mod;
	int fac = 1, sumF = 1, F, ans = 0, lstF = 0;
	// lstF = F_{i-1}, F = F_{i}, sumF = \sum_{i=0}^{n-2} F_i
	for(int i = 2; i <= n; ++i) {
		fac = 1ll * fac * i % mod; // fac = i!
		F = 1ll * sumF * inv[i] % mod * inv[2] % mod;
		ans = (ans + 1ll * F * fac % mod * fac % mod) % mod;
		sumF = (sumF + lstF) % mod;
		lstF = F;// F_{i-1}
	}
	printf("%d\n", ans); // ans = \sum_{i=1}^{n} f_i
}

代码源自 $\text{Ark}$ 。这个链接只是个人主页罢了。