結果

問題 No.1421 国勢調査 (Hard)
コンテスト
ユーザー chocorusk
提出日時 2021-03-05 22:30:34
言語 C++17
(gcc 13.3.0 + boost 1.87.0)
結果
WA  
実行時間 -
コード長 2,262 bytes
コンパイル時間 3,454 ms
コンパイル使用メモリ 183,972 KB
最終ジャッジ日時 2025-01-19 11:28:44
ジャッジサーバーID
（参考情報） 		judge2 / judge1
このコードへのチャレンジ
（要ログイン）
ファイルパターン 結果
sample AC * 2
other AC * 16 WA * 14
権限があれば一括ダウンロードができます

ソースコード

diff # 

#include <cstdio>
#include <cstring>
#include <iostream>
#include <string>
#include <cmath>
#include <bitset>
#include <vector>
#include <map>
#include <set>
#include <queue>
#include <deque>
#include <algorithm>
#include <complex>
#include <unordered_map>
#include <unordered_set>
#include <random>
#include <cassert>
#include <fstream>
#include <utility>
#include <functional>
#include <time.h>
#include <stack>
#include <array>
#include <list>
#include <atcoder/all>
#define popcount __builtin_popcount
using namespace std;
using namespace atcoder;
typedef long long ll;
typedef pair<int, int> P;
const int MAX=50; //max of width
using BS=bitset<MAX+1>;
vector<BS> Gauss_Jordan(const vector<BS> &v){
    int n=v.size();
    vector<BS> ret=v;
    int d=0;
    for(int k=0; k<MAX; k++){
        int p=-1;
        for(int i=d; i<n; i++){
            if(ret[i][k]){
                p=i;
                break;
            }
        }
        if(p==-1) continue;
        swap(ret[d], ret[p]);
        for(int i=0; i<n; i++){
            if(i==d) continue;
            if(ret[i][k]) ret[i]^=ret[d];
        }
        d++;
    }
    return ret;
}
int main()
{
    int n, m; cin>>n>>m;
    vector<BS> v(m, BS(0));
    int y[10010];
    for(int i=0; i<m; i++){
        int a; cin>>a;
        for(int j=0; j<a; j++){
            int b; cin>>b;
            b--;
            v[i][b]=1;
        }
        cin>>y[i];
    }
    int ans[50]={};
    for(int i=0; i<30; i++){
        for(int j=0; j<m; j++){
            if(y[j]&(1<<i)) v[j][MAX]=1;
            else v[j][MAX]=0;
        }
        auto ret=Gauss_Jordan(v);
        bool x[50]={};
        for(int j=m-1; j>=0; j--){
            int p=-1;
            bool s=0;
            for(int k=0; k<n; k++){
                if(p!=-1){
                    if(ret[j][k]){
                        s^=x[k];
                    }
                }
                if(ret[j][k]){
                    p=k;
                }
            }
            if(p==-1 && ret[j][MAX]){
                cout<<-1<<endl;
                return 0;
            }
            x[p]=(s^ret[j][MAX]);
        }
        for(int j=0; j<n; j++) if(x[j]) ans[j]^=(1<<i);
    }
    for(int i=0; i<n; i++){
        cout<<ans[i]<<endl;
    }
    return 0;
}
0