「2020-2021 集训队作业」Yet Another Linear Algebra Problem（行列式，Binet-Cauchy 公式）

最新推荐文章于 2021-09-26 11:55:30 发布

DD(XYX)

最新推荐文章于 2021-09-26 11:55:30 发布

阅读量399

点赞数 2

分类专栏：数学

本文链接：https://blog.csdn.net/weixin_43960414/article/details/119008209

版权

数学专栏收录该内容

118 篇文章 1 订阅

订阅专栏

题面

出题人： $\tt TLY$ 太阳神：Tiw_Air_OAO

${\tt「2020—2021~集训队作业」Yet~Another~Linear~Algebra~Problem}\\\\ {\tt_{传统~~~~~1000\,ms~~~~~512\,MiB}}$

题目描述

您需要解决两个独立（但类似）的子问题：

问题一：给定 $n$ 个在 $G F (3)$ （数字都模 $3$ ）上的 $m$ 维向量，记它们张成的线性空间为 $V$ 。求从 $n$ 个向量中选出一组向量，使得它们是 $V$ 的基的方案数。对 $3$ 取模。

问题二：给定 $n$ 个在 $G F (2)$ 上的 $m$ 维向量，记它们张成的线性空间为 $V$ 。其中第 $i$ 个向量有颜色 $c_i$ 。求从每种颜色中恰好选出一个向量，使得它们是 $V$ 的基的方案数。对 $2$ 取模。

注：为了凸显主要矛盾而忽略次要矛盾，保证 $V$ 的维数为 $m$ 。

输入格式

第一行包含一个正整数 $t a s k i d$ ，表示需要解决的问题编号。

第二行包含两个正整数 $n, m$ ，含义见上。

接下来输入 $n$ 行：

如果 $t a s k i d = 1$ ，则第 $i$ 行包含 $m$ 个非负整数 $v_{i,1},v_{i,2},\ldots,v_{i,m}$ ，描述了第 $i$ 个向量。

如果 $t a s k i d = 2$ ，则第 $i$ 行包含 $m + 1$ 个非负整数 $v_{i,1},v_{i,2},\ldots,v_{i,m},c_i$ ，描述了第 $i$ 个向量与它的颜色。

输出格式

输出一行一个正整数表示答案。

样例 1

输入

输出

样例 2

输入

输出

样例 3

输入

输出

样例 4

输入

输出

样例 5

输入

输出

数据范围与提示

对于 $100\%$ 的数据， $taskid\in\{1,2\},1\leq n,m\leq500$ 。

当 $t a s k i d = 1$ 时，则 $v_{i,j}\in\{0,1,2\}$ 。

当 $t a s k i d = 2$ 时，则 $v_{i,j}\in\{0,1\},c_i\in[1,m]$ 。

${\rm subtask~1(50~pts)}:taskid=1$ 。

${\rm subtask~2(50~pts)}:taskid=2$ 。

$_{{\rm LibreOJ}~{\tt Powered~by}~{\rm SYZOJ~NG}}$

题解

用到的知识很多。

首先，我们知道一个 $n\times n$ 的矩阵是满秩的，等价于它的行列式非 $0$ ，那么选 $m$ 个向量，是否能成为 $V$ 的一组基底，就取决于它们拼成的 $m\times m$ 矩阵是否行列式非 $0$ 。

对于 $\rm subtask~1$ ，我们可以把答案表示为：
$\sum_{1\leq i_1<i_2<...<i_m\leq n} \left[~\left|\begin{matrix} —\overset{\rightarrow}{v_{i_1}}—\\ —\overset{\rightarrow}{v_{i_2}}—\\ \vdots\\ —\overset{\rightarrow}{v_{i_m}}—\\ \end{matrix}\right|\not\equiv0~\right]\mod 3$

中间那个条件符号不好处理啊！但是我们知道 $1^2,2^2$ 在模 $3$ 意义下都是等于 $1$ 的，也就是说，把行列式平方，刚好同余于该条件表达式！那么我们就平方试试，顺便把其中一个行列式矩阵转置一下：
$\sum_{1\leq i_1<i_2<...<i_m\leq n} \left|\begin{matrix} |&|&|&|\\ \overset{\rightarrow}{v_{i_1}}&\overset{\rightarrow}{v_{i_2}}&\dots&\overset{\rightarrow}{v_{i_m}}\\ |&|&|&|\\ \end{matrix}\right| \cdot \left|\begin{matrix} —\overset{\rightarrow}{v_{i_1}}—\\ —\overset{\rightarrow}{v_{i_2}}—\\ \vdots\\ —\overset{\rightarrow}{v_{i_m}}—\\ \end{matrix}\right| \mod 3$

再运用 $\tt Binet\_Cauchy$ 公式，可以知道上式就等于
$\det\left(\left(\begin{matrix} |&|&|&|\\ \overset{\rightarrow}{v_{1}}&\overset{\rightarrow}{v_{2}}&\dots&\overset{\rightarrow}{v_{n}}\\ |&|&|&|\\ \end{matrix}\right) \cdot \left(\begin{matrix} —\overset{\rightarrow}{v_{1}}—\\ —\overset{\rightarrow}{v_{2}}—\\ \vdots\\ —\overset{\rightarrow}{v_{n}}—\\ \end{matrix}\right)\right)$

因此求这个行列式就行了。

对于 $\rm subtask~2$ ，在模 $2$ 意义下，并没有变得更简单，而是多加了颜色条件。我们可以这样翻译这个条件：所选的 $m$ 个向量的颜色必须互不相同，每种颜色都要有向量。

这让我们又想到了一个关于行列式的性质：若行列式的某一行或某一列全为 $0$ ，那么行列式为 $0$ 。

因此，我们可以用一个矩阵 $A$ 表示颜色， $A_{i,j}$ 定义为向量 $i$ 的颜色是( $1$ )否( $0$ )为 $j$ 。每种颜色都要有向量，也就等价于 $A$ 中对应的 $m\times m$ 的子矩阵每行每列都有恰好一个 $1$ ，行列式不为 $0$ 。

由于是取模 $2$ ，向量组成的行列式本身就等价于不等于 $0$ 的条件表达式，可以不用平方。我们就把 $A$ 的行列式乘过来表示答案：
$\sum_{1\leq i_1<i_2<...<i_m\leq n} \left|\begin{matrix} |&|&|&|\\ \overset{\rightarrow}{v_{i_1}}&\overset{\rightarrow}{v_{i_2}}&\dots&\overset{\rightarrow}{v_{i_m}}\\ |&|&|&|\\ \end{matrix}\right| \cdot \left|\begin{matrix} A_{i_1,1}&A_{i_1,2}&\dots&A_{i_1,m}\\ A_{i_2,1}&A_{i_2,2}&\dots&A_{i_2,m}\\ \vdots&\vdots&\ddots&\vdots\\ A_{i_m,1}&A_{i_m,2}&\dots&A_{i_m,m}\\ \end{matrix}\right| \mod 2$

又可以使用 $\tt Binet\_Cauchy$ 公式了，把它变为
$\det\left(\left(\begin{matrix} |&|&|&|\\ \overset{\rightarrow}{v_{1}}&\overset{\rightarrow}{v_{2}}&\dots&\overset{\rightarrow}{v_{n}}\\ |&|&|&|\\ \end{matrix}\right) \cdot A\right)$

直接求行列式。

复杂度都是 $O(n^3)$ 。

CODE

$^{_{代码并不一定是博客终结的地方}}$

#include<vector>
#include<cstdio>
#include<cstring>
#include<iostream>
#include<algorithm>
using namespace std;
#define MAXN 505
#define ENDL putchar('\n')
#define LL long long
#define DB double
#define lowbit(x) ((-x) & (x))
LL read() {
	LL f = 1,x = 0;char s = getchar();
	while(s < '0' || s > '9') {if(s=='-')f = -f;s = getchar();}
	while(s >= '0' && s <= '9') {x=x*10+(s-'0');s = getchar();}
	return f * x;
}
int n,m,i,j,s,o,k;
int MOD = 1;
struct mat{
	int n,m;
	int s[500][500];
	mat(){n=m=0;memset(s,0,sizeof(s));}
	inline void set(int N,int M) {n=N;m=M;}
}X,A;
inline mat operator * (mat a,mat b) {
	mat c;c.set(a.n,b.m);
	for(int i = 0;i < a.n;i ++) {
		for(int k = 0;k < a.m;k ++) {
			if(a.s[i][k])
			for(int j = 0;j < b.m;j ++) {
				c.s[i][j] += a.s[i][k]*b.s[k][j];
			}
		}
	}
	for(int i = 0;i < c.n;i ++) {
		for(int j = 0;j < c.m;j ++) c.s[i][j] %= MOD;
	}return c;
}
mat T_mat(mat a) {
	for(int i = 0;i < a.n;i ++) {
		for(int j = i;j < a.m;j ++) {
			swap(a.s[i][j],a.s[j][i]);
		}
	}
	swap(a.n,a.m); return a;
}
int det(mat a) {
	if(a.n != a.m) return 0;
	int as = 1;
	for(int i = 0;i < a.n;i ++) {
		if(!(a.s[i][i]%MOD))
		for(int j = i;j < a.n;j ++) {
			if(a.s[j][i]%MOD) {swap(a.s[i],a.s[j]);as = MOD-as;break;}
		}
		if(!(a.s[i][i]%MOD)) return 0;
		for(int j = i+1;j < a.n;j ++) {
			if(a.s[j][i]%MOD) {
				for(int k = a.m-1;k >= i;k --) {
					(a.s[j][k] += MOD-a.s[i][k]*a.s[i][i]*a.s[j][i]%MOD) %= MOD;
				}
			}
		}
	}
	for(int i = 0;i < a.n;i ++) as = as * a.s[i][i] % MOD;
	return (as+MOD)%MOD;
}
int main() {
	int op = read();
	MOD = 4-op;
	n = read(); m = read();
	X.set(m,n);
	A.set(n,m);
	for(int i = 0;i < n;i ++) {
		for(int j = 0;j < m;j ++) {
			X.s[j][i] = read()%MOD;
		}
		if(op == 2) A.s[i][read()-1] = 1;
	}
	printf("%d\n",op == 1 ? det(X * T_mat(X)) : det(X * A));
	return 0;
}