codeforces736D. Permutations(线性代数)

题意

题目打错了

Sol

神仙题Orz

构造矩阵B,使得B[b[i]][a[i]] = 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
 
*/

本文参与腾讯云自媒体分享计划,欢迎正在阅读的你也加入,一起分享。

发表于

我来说两句

0 条评论
登录 后参与评论

相关文章

来自专栏数据结构与算法

BZOJ2118: 墨墨的等式(最短路 数论)

这里有一篇讲的非常好的博客https://blog.csdn.net/w4149/article/details/66476606?locationNum=3&...

661
来自专栏hanlp学习笔记

hanlp源码解析之中文分词算法

词图指的是句子中所有词可能构成的图。如果一个词A的下一个词可能是B的话,那么A和B之间具有一条路径E(A,B)。一个词可能有多个后续,同时也可能有多个前驱,它们...

1572
来自专栏xingoo, 一个梦想做发明家的程序员

[翻译]CURAND Libaray--Host API--(2)

2.3 返回值 所有的CURAND host端的函数返回值都是curandStatus_t.如果调用没有错误,则返回成功,即返回值为CURAND_STATUS_...

22410
来自专栏GreenLeaves

SQL学习之汇总数据之聚集函数

一、 1、我们经常需要汇总数据而不用把他们实际检索出来,为此SQL提供了专门的函数,以便于分析数据和报表生成,这些函数的功能有: (1)确定表中行数(或者满足单...

1995
来自专栏owent

POJ PKU 1990 MooFest 解题报告

561
来自专栏图形学与OpenGL

机械版CG 实验4 裁剪

了解二维图形裁剪的原理(点的裁剪、直线的裁剪、多边形的裁剪),利用VC+OpenGL实现直线的裁剪算法。

1721
来自专栏linux、Python学习

必备神技能 | MySQL 查找删除重复行

本文讲述如何查找数据库里重复的行。这是初学者十分普遍遇到的问题。方法也很简单。这个问题还可以有其他演变,例如,如何查找“两字段重复的行”(#mysql IRC ...

710
来自专栏hadoop学习笔记

hanlp源码解析之中文分词算法详解

词图指的是句子中所有词可能构成的图。如果一个词A的下一个词可能是B的话,那么A和B之间具有一条路径E(A,B)。一个词可能有多个后续,同时也可能有多个前驱,它们...

1043
来自专栏每日一篇技术文章

OpenGLES_实战04_教你绘制球体

写这篇文章主要给初学者一个绘制球体的思路,苹果给我们封装的类,帮助我们简化了不少代码,如果纯OpenGL 做这样一个练习代码量还是挺多的。

1791
来自专栏有趣的Python和你

MySQL学习之联结表内联结左联结右联结

1165

扫码关注云+社区