标签:技术分享 线性代数 题目 set alt stream ret pre image
$m \leqslant 500000$,题目打错了
神仙题Orz
构造矩阵$B$,使得$B[b[i]][a[i]] = 1$
那么他的行列式的奇偶性也就对应了生成排列数列数量的奇偶性(定义)
删除一个位置相当于去掉对答案的贡献,也就是代数余子式的值
代数余子式可以由伴随矩阵求出$A^{*} = |A| A^{-1}$
这里只需要奇偶性,因此不需要求出实际行列式的值。
矩阵可以用bitset加速,可以过掉这个题
#include<cstdio> #include<bitset> #include<iostream> using namespace std; const int MAXN = 2001; inline int read() { char c = getchar(); int x = 0, f = 1; while(c < ‘0‘ || c > ‘9‘) {if(c == ‘-‘) f = -1; c = getchar();} while(c >= ‘0‘ && c <= ‘9‘) x = x * 10 + c - ‘0‘, c = getchar(); return x * f; } int N, M; bitset<MAXN * 2 + 10> b[MAXN]; int x[500001], y[500001]; int main() { N = read(); M = read(); for(int i = 1; i <= N; i++) b[i][i + N] = 1; for(int i = 1; i <= M; i++) { x[i] = read(), y[i] = read(); b[x[i]][y[i]] = 1; } for(int i = 1, j; i <= N; i++) { for(j = i; j <= N; j++) if(b[j][i]) {swap(b[i], b[j]); break;} for(int k = 1; k <= N; k++) if(b[k][i] && (k != i)) b[k] ^= b[i]; } for(int i = 1; i <= N; i++, puts("")) for(int j = 1; j <= 2 * N; j++) cout << b[i][j] << " "; return 0; } /* 3 7 1 1 1 3 2 2 2 3 3 1 3 2 3 3 */
codeforces736D. Permutations(线性代数)
标签:技术分享 线性代数 题目 set alt stream ret pre image
原文地址:https://www.cnblogs.com/zwfymqz/p/9416394.html