結果
問題 | No.1582 Vertexes vs Edges |
ユーザー | 沙耶花 |
提出日時 | 2021-07-02 21:47:01 |
言語 | C++17 (gcc 13.3.0 + boost 1.87.0) |
結果 |
AC
|
実行時間 | 156 ms / 2,000 ms |
コード長 | 1,924 bytes |
コンパイル時間 | 3,872 ms |
コンパイル使用メモリ | 264,304 KB |
最終ジャッジ日時 | 2025-01-22 15:52:02 |
ジャッジサーバーID (参考情報) |
judge2 / judge3 |
(要ログイン)
ファイルパターン | 結果 |
---|---|
sample | AC * 3 |
other | AC * 36 |
ソースコード
#include <stdio.h> #include <bits/stdc++.h> #include <atcoder/all> using namespace atcoder; using mint = modint1000000007; using namespace std; #define rep(i,n) for (int i = 0; i < (n); ++i) #define Inf 1000000001 struct lca{ vector<int> depth; vector<vector<int>> parents; int max_j=18; lca(int n,vector<vector<int>> &E){ depth.assign(E.size(),-1); parents.assign(E.size(),vector<int>(max_j,-1)); rep(i,100){ if((1<<i)>E.size()){ max_j = i; break; } } stack<int> s; s.push(n); depth[n] = 0; while(s.size()!=0){ int k = s.top(); s.pop(); for(int i=0;i<E[k].size();i++){ int x = E[k][i]; if(depth[x]!=-1)continue; s.push(x); depth[x] = depth[k]+1; for(int j=0;j<max_j;j++){ if(j==0){ parents[x][j] = k; } else{ parents[x][j] = parents[parents[x][j-1]][j-1]; } if(parents[x][j] == -1)break; } } } } int kth_parent(int u,int k){ for(int i=0;i<max_j;i++){ if(k==0)break; if(u==-1)break; if(k%2){ u = parents[u][i]; } k/=2; } return u; } int query(int u,int v){ if(depth[u]<depth[v])swap(u,v); u = kth_parent(u,depth[u]-depth[v]); if(u==v){ return u; } for(int i=max_j-1;i>=0;i--){ if(parents[u][i]!=parents[v][i]){ u = parents[u][i]; v = parents[v][i]; } } return parents[u][0]; } int get_distance(int u,int v){ return depth[u]+depth[v]-2*depth[query(u,v)]; } }; int main(){ int N; cin>>N; vector<vector<int>> E(N); rep(i,N-1){ int u,v; cin>>u>>v; u--;v--; E[u].push_back(v); E[v].push_back(u); } lca L(0,E); mf_graph<int> G(N+2); int S = N,T = N+1; rep(i,N){ if(L.get_distance(0,i)%2==0){ rep(j,E[i].size()){ int u = i,v = E[i][j]; G.add_edge(u,v,1); } G.add_edge(S,i,1); } else{ G.add_edge(i,T,1); } } cout<<G.flow(S,T)<<endl; return 0; }