CF736D Permutations

CF736D Permutations
首先领接矩阵表示哪个数可以放哪个位置，设为1，然后它所构成的$n$阶积和式就是每个条件都成立的合法方案数
然后考虑一个条件失效，那么通过那个位置的代数余子式，就可以得到这一行失效后的合法方案数，又已知伴随矩阵各个位置代表各个代数余子式，$A$*=$|A|\cdot A^{-1}$ 由于$|A|$为奇数，所以$A$*给个位置上奇偶性只取决于$A^{-1}$，于是高斯消元求$A^{-1}$即可

由于每个位置都为0/1，所以可以用bitset优化

const int N=2005;
int n,m;
std::bitset<2*N> a[N];
void gauss(){
    for (int i=1;i<=n;i++){
        int now=0;
        for (int j=i;j<=n;j++){
            if (a[j][i]){
                now=j;
                break;
            }
        }
        if (now!=i){
            std::swap(a[now],a[i]);
        }
        for (int j=1;j<=n;j++){
            if (j!=i&&a[j][i]){
                a[j]^=a[i];
            }
        }
    }
}
void solve(){
    std::cin>>n>>m;
    std::vector<int> x(m+1),y(m+1);
    for (int i=1;i<=m;i++){
        std::cin>>x[i]>>y[i];
        a[x[i]][y[i]]=1;
    }
    for (int i=1;i<=n;i++){
        a[i][i+n]=1;
    }

    gauss();
    for (int i=1;i<=m;i++){
        if (a[y[i]][x[i]+n]){
            std::cout<<"NO\n";
        }else{
            std::cout<<"YES\n";
        }
    }
}