結果
問題 | No.1421 国勢調査 (Hard) |
ユーザー |
![]() |
提出日時 | 2021-03-05 22:34:35 |
言語 | C++17 (gcc 13.3.0 + boost 1.87.0) |
結果 |
AC
|
実行時間 | 173 ms / 2,000 ms |
コード長 | 2,271 bytes |
コンパイル時間 | 3,038 ms |
コンパイル使用メモリ | 183,228 KB |
最終ジャッジ日時 | 2025-01-19 11:30:35 |
ジャッジサーバーID (参考情報) |
judge1 / judge1 |
(要ログイン)
ファイルパターン | 結果 |
---|---|
sample | AC * 2 |
other | AC * 30 |
ソースコード
#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_popcountusing namespace std;using namespace atcoder;typedef long long ll;typedef pair<int, int> P;const int MAX=50; //max of widthusing 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(p==-1 && 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;}