結果
問題 | No.3009 Union-Find でつながろう! |
ユーザー |
![]() |
提出日時 | 2025-01-17 22:18:24 |
言語 | C++17 (gcc 13.3.0 + boost 1.87.0) |
結果 |
AC
|
実行時間 | 303 ms / 2,000 ms |
コード長 | 2,212 bytes |
コンパイル時間 | 4,914 ms |
コンパイル使用メモリ | 259,176 KB |
実行使用メモリ | 6,144 KB |
最終ジャッジ日時 | 2025-01-17 22:18:49 |
合計ジャッジ時間 | 8,747 ms |
ジャッジサーバーID (参考情報) |
judge4 / judge2 |
(要ログイン)
ファイルパターン | 結果 |
---|---|
sample | AC * 3 |
other | AC * 27 |
ソースコード
/** author: shobonvip created: 2025.01.17 22:11:07 **/ #include<bits/stdc++.h> using namespace std; //* ATCODER #include<atcoder/all> using namespace atcoder; typedef modint998244353 mint; //*/ /* BOOST MULTIPRECISION #include<boost/multiprecision/cpp_int.hpp> using namespace boost::multiprecision; //*/ typedef long long ll; #define rep(i, s, n) for (int i = (int)(s); i < (int)(n); i++) #define rrep(i, s, n) for (int i = (int)(n)-1; i >= (int)(s); i--) #define all(v) v.begin(), v.end() template <typename T> bool chmin(T &a, const T &b) { if (a <= b) return false; a = b; return true; } template <typename T> bool chmax(T &a, const T &b) { if (a >= b) return false; a = b; return true; } template <typename T> T max(vector<T> &a){ assert(!a.empty()); T ret = a[0]; for (int i=0; i<(int)a.size(); i++) chmax(ret, a[i]); return ret; } template <typename T> T min(vector<T> &a){ assert(!a.empty()); T ret = a[0]; for (int i=0; i<(int)a.size(); i++) chmin(ret, a[i]); return ret; } template <typename T> T sum(vector<T> &a){ T ret = 0; for (int i=0; i<(int)a.size(); i++) ret += a[i]; return ret; } typedef bitset<4096> B; int main(){ ios_base::sync_with_stdio(false); cin.tie(NULL); int m,n; cin >> m >> n; vector<B> a(m); rep(i,0,m) { rep(j,0,n){ char x;cin >> x; a[i][j] = (x=='1'); } } vector<B> e(n); rep(i,0,n){ rep(j,0,n){ e[i][j]=1; } } rep(i,0,n) { rep(j,0,m){ if(a[j][i]){ e[i]&=a[j]; } } } // scc atcoder::scc_graph scc(n); for (int i=0; i<n; i++){ for (int j=0; j<n; j++){ if (e[i][j]){ scc.add_edge(i, j); } } } vector<vector<int>> new_graph = scc.scc(); vector<int> cores(n); int k = new_graph.size(); for (int i=0; i<k; i++){ for (int j: new_graph[i]){ cores[j] = i; } } dsu uf(k); for (int i=0; i<n; i++){ for (int j=0; j<n; j++){ if (cores[i] == cores[j]){ continue; } if (e[i][j]){ uf.merge(cores[i],cores[j]); } } } int cnt=0; rep(i,0,k){ if(uf.leader(i)==i)cnt++; } cout<<cnt<<endl; }