1. 程式人生 > >洛谷P2447 [SDOI2010]外星千足蟲(異或方程組)

洛谷P2447 [SDOI2010]外星千足蟲(異或方程組)

找到 art pre pro code main getc earth getchar()

題意

題目鏈接

Sol

異或高斯消元的板子題。

bitset優化一下,復雜度\(O(\frac{nm}{32})\)

找最優解可以考慮高斯消元的過程,因為異或的特殊性質,每次向下找的時候找到第一個1然後交換就行,這樣顯然是最優的

#include<bits/stdc++.h>
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> b[MAXN];
void Gauss() {
    int ans = 0;
    for(int i = 1; i <= N; i++) {
        int j = i;
        while(!b[j][i] && j < M + 1) 
            j++;
        if(j == M + 1) {puts("Cannot Determine"); return ;}
        ans = max(ans, j);
        swap(b[i], b[j]);
        for(int j = 1; j <= M; j++) {
            if(i == j || !b[j][i]) continue;
            b[j] ^= b[i];
        }
    }
    printf("%d\n", ans);
    for(int i = 1; i <= N; i++)
        puts(!b[i][N + 1] ? "Earth" : "?y7M#");
}
int main() {
    N = read(); M = read();
    for(int i = 1; i <= M; i++) {
        string s; cin >> s; 
        b[i][N + 1] = read();
        for(int j = 1; j <= N; j++) b[i][j] = (s[j - 1] == '0' ? 0 : 1);
    }
    Gauss();
    return 0;
}

洛谷P2447 [SDOI2010]外星千足蟲(異或方程組)