結果

問題 No.922 東北きりきざむたん
コンテスト
ユーザー tute7627
提出日時 2019-11-08 22:20:31
言語 C++14
(gcc 13.3.0 + boost 1.87.0)
結果
AC  
実行時間 237 ms / 2,000 ms
コード長 5,131 bytes
コンパイル時間 2,520 ms
コンパイル使用メモリ 192,004 KB
実行使用メモリ 51,252 KB
最終ジャッジ日時 2024-09-15 01:35:49
合計ジャッジ時間 6,523 ms
ジャッジサーバーID
（参考情報） 		judge6 / judge4
このコードへのチャレンジ
（要ログイン）
ファイルパターン 結果
sample AC * 4
other AC * 26
権限があれば一括ダウンロードができます

ソースコード

diff # 

#include<bits/stdc++.h>
using namespace std;

#define endl '\n'
#define ALL(a)  (a).begin(),(a).end()
#define ALLR(a)  (a).rbegin(),(a).rend()
#define spa << " " <<
#define test cout<<"test"<<endl;
#define fi first
#define se second
#define MP make_pair
#define PB push_back
#define EB emplace_back
#define rep(i,n,m) for(ll i = n; i < (ll)(m); i++)
#define rrep(i,n,m) for(ll i = m - 1; i >= (ll)(n); i--)
using ll = long long;
using ld = long double;
const ll MOD = 1e9+7;
//const ll MOD = 998244353;
const ll INF = 1e18;
using P = pair<ll, ll>;
template<typename T>
void chmin(T &a,T b){if(a>b)a=b;}
template<typename T>
void chmax(T &a,T b){if(a<b)a=b;}
void pmod(ll &a,ll b){a=(a+b)%MOD;}
void pmod(ll &a,ll b,ll c){a=(b+c)%MOD;}
void qmod(ll &a,ll b){a=(a*b)%MOD;}
void qmod(ll &a,ll b,ll c){a=(b*c)%MOD;}
ll median(ll a,ll b, ll c){return a+b+c-max({a,b,c})-min({a,b,c});}
void ans1(bool x){if(x) cout<<"Yes"<<endl;else cout<<"No"<<endl;}
void ans2(bool x){if(x) cout<<"YES"<<endl;else cout<<"NO"<<endl;}
void ans3(bool x){if(x) cout<<"Yay!"<<endl;else cout<<":("<<endl;}
template<typename T1,typename T2>
void ans(bool x,T1 y,T2 z){if(x)cout<<y<<endl;else cout<<z<<endl;}  
template<typename T>
void debug(vector<vector<T>>&v,ll h,ll w){for(ll i=0;i<h;i++)
{cout<<v[i][0];for(ll j=1;j<w;j++)cout spa v[i][j];cout<<endl;}};
void debug(vector<string>&v,ll h,ll w){for(ll i=0;i<h;i++)
{for(ll j=0;j<w;j++)cout<<v[i][j];cout<<endl;}};
template<typename T>
void debug(vector<T>&v,ll n){if(n!=0)cout<<v[0];
for(ll i=1;i<n;i++)cout spa v[i];cout<<endl;};
template<typename T>
vector<vector<T>>vec(ll x, ll y, T w){
  vector<vector<T>>v(x,vector<T>(y,w));return v;}
ll gcd(ll x,ll y){ll r;while(y!=0&&(r=x%y)!=0){x=y;y=r;}return y==0?x:y;}
vector<ll>dx={1,0,-1,0,1,1,-1,-1};
vector<ll>dy={0,1,0,-1,1,-1,1,-1};
template<typename T>
vector<T> make_v(size_t a,T b){return vector<T>(a,b);}
template<typename... Ts>
auto make_v(size_t a,Ts... ts){
  return vector<decltype(make_v(ts...))>(a,make_v(ts...));
}

vector<ll>val(110000);
struct LCA{
  ll datasize = 21;
  ll n,root;//rootを根とする頂点数nの根付き木
  vector<ll>depth;//0から
  vector<vector<ll>>G,data;
  vector<ll>group;
  vector<P>dp;
  ll minbuf,minsum;
  LCA(vector<vector<ll>>g,ll r):G(g),n(g.size()),root(r){
    data.assign(n,vector<ll>(datasize,-1));
    depth.assign(n,-1LL);
    group.assign(n,-1LL);
    dp.assign(n,MP(0LL,0LL));
    minsum=0;
    build();
  }
  void build(){
    ll cnt=0;
    rep(i,0,n){
      if(depth[i]==-1){
        minbuf=INF;
        dfs(i,0,cnt);
        cnt++;
        dfsmin(i,0,MP(dp[i].fi,dp[i].fi+dp[i].se));
        minsum+=minbuf;
        //cout<<i spa minbuf<<endl;
      }
    }
    for(ll i=1;i<datasize;i++){
      for(ll j=0;j<n;j++){
        if(data[j][i-1]==-1)data[j][i]=-1;
        else data[j][i]=data[data[j][i-1]][i-1];
      }
    }
  }
  pair<ll,ll> dfs(ll k,ll d, ll x){
    depth[k]=d;
    group[k]=x;
    ll num=0;
    ll sumbuf=0;
    for(ll i=0;i<G[k].size();i++){
      if(depth[G[k][i]]==-1){
        data[G[k][i]][0]=k;
        auto p=dfs(G[k][i],d+1,x);
        num+=p.fi,sumbuf+=p.se;
      }
    }
    dp[k]=MP(num+val[k],sumbuf+num);
    return MP(num+val[k],sumbuf+num);
  }
  void dfsmin(ll k,ll d,P par){
    auto tmp=par;
    tmp.fi-=dp[k].fi;
    tmp.se-=dp[k].se+dp[k].fi;
    tmp.se+=tmp.fi;
    chmin(minbuf,tmp.se+dp[k].se);
    //cout<<k spa tmp.fi spa tmp.se spa par.fi spa par.se<<endl;
    //cout<<dp[k].fi spa dp[k].se spa minbuf<<endl;
    rep(i,0,G[k].size()){
      if(depth[G[k][i]]==d+1){
        dfsmin(G[k][i],d+1,MP(dp[k].fi+tmp.fi,dp[k].se+tmp.se));
      }
    }
  }
  ll deep(ll x){
    return x>=0?depth[x]:-1;
  }
  ll query(ll x,ll y){
    if(depth[x]<depth[y])swap(x,y);
    ll tmp = datasize - 1;
    while(depth[x]!=depth[y]){
      while(tmp>=1&&depth[y]>=deep(data[x][tmp]))tmp--;
      x = data[x][tmp];
    }
    tmp = datasize - 1;
    while(x!=y){
      while(tmp>=1&&data[x][tmp]==data[y][tmp])tmp--;
      x = data[x][tmp];
      y = data[y][tmp];
    }
    return x;
  }
};

struct UnionFind {
  vector<ll> data;
  UnionFind(ll size) : data(size, -1) { }
  bool unite(ll x, ll y) {
    x = root(x); y = root(y);
    if (x != y) {
      if (data[y] < data[x]) swap(x, y);
      data[x] += data[y]; data[y] = x;
    }
    return x != y;
  }
  bool find(ll x, ll y) {
    return root(x) == root(y);
  }
  ll root(ll x) {
    return data[x] < 0 ? x : data[x] = root(data[x]);
  }
  ll size(ll x) {
    return -data[root(x)];
  }
};
int main(){
  cin.tie(nullptr);
  ios_base::sync_with_stdio(false);
  //cout<<fixed<<setprecision(10);
  ll res=0,buf=0;
  bool judge = true;
  ll n,m,q;cin>>n>>m>>q;
  vector<vector<ll>>g(n);
  UnionFind uf(n);
  rep(i,0,m){
    ll u,v;cin>>u>>v;
    u--;v--;
    g[u].PB(v);
    g[v].PB(u);
    uf.unite(u,v);
  }
  vector<P>p;
  rep(i,0,q){
    ll a,b;cin>>a>>b;
    a--;b--;
    if(uf.find(a,b))p.EB(a,b);
    else val[a]++,val[b]++;
  }
  LCA lca(g,0);
  rep(i,0,p.size()){
    lca.minsum+=lca.deep(p[i].fi)+lca.deep(p[i].se)
      -2*lca.deep(lca.query(p[i].fi,p[i].se));
  }
  cout<<lca.minsum<<endl;
  return 0;
}
0