码迷,mamicode.com
首页 > 其他好文 > 详细

codeforces736D. Permutations(线性代数)

时间:2018-08-04 12:05:08      阅读:130      评论:0      收藏:0      [点我收藏+]

标签:技术分享   线性代数   题目   set   alt   stream   ret   pre   image   

题意

技术分享图片

$m \leqslant 500000$,题目打错了

 

Sol

神仙题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

(0)
(0)
   
举报
评论 一句话评论(0
登录后才能评论!
© 2014 mamicode.com 版权所有  联系我们:gaon5@hotmail.com
迷上了代码!