Submission #385915

#TimeUsernameProblemLanguageResultExecution timeMemory
385915ogibogi2004Toll (APIO13_toll)C++14
78 / 100
2598 ms20996 KiB
#include<bits/stdc++.h> using namespace std; #define ll long long const ll INF=1e15; const ll MAXN=3e5+6; const ll MAXM=3e5+6; struct edge { int u,v; ll w; int idx; bool operator<(edge const& other)const { return w<other.w; } }; ll ans,sum; int n,m,k; bool alwaysin[MAXN]; vector<pair<int,int> >g[MAXN]; ll pari[MAXN]; ll pari1[MAXN]; vector<edge>sometimes_in; edge old[MAXM]; edge nw[20]; vector<int>comps; struct DSU { int par[MAXN],sz[MAXN]; ll h[MAXN]; void init() { for(int i=1;i<MAXN;i++) { par[i]=i;sz[i]=1; h[i]=pari[i]; } } void init(vector<int>v) { for(auto xd:v) { par[xd]=xd; sz[xd]=1; } } int getRoot(int u) { if(u==par[u])return u; return par[u]=getRoot(par[u]); } void join(int p,int q) { if(sz[p]>sz[q]) { par[q]=p; sz[p]+=sz[q]; h[p]+=h[q]; } else { par[p]=q; sz[q]+=sz[p]; h[q]+=h[p]; } } }dsu1,dsu2; int lei[MAXN]; ll atmost[MAXN]; int par[MAXN]; int depth[MAXN]; ll subtree[MAXN]; void dfs1(int u,int p,int last_edge_idx=-69) { atmost[u]=INF; depth[u]=depth[p]+1; subtree[u]=pari1[u]; lei[u]=last_edge_idx; par[u]=p; for(auto xd:g[u]) { if(xd.first==p)continue; dfs1(xd.first,u,xd.second); subtree[u]+=subtree[xd.first]; } } void takovai(edge e) { int x=e.u,y=e.v; if(depth[x]<depth[y])swap(x,y); while(depth[x]>depth[y]) { atmost[x]=min(atmost[x],e.w); x=par[x]; } while(x!=y) { atmost[x]=min(atmost[x],e.w); atmost[y]=min(atmost[y],e.w); x=par[x];y=par[y]; } } void dfs(int u,int p) { for(auto xd:g[u]) { if(xd.first==p)continue; if(xd.second>m) { //cout<<"%%%\n"; sum+=atmost[xd.first]*subtree[xd.first]; } dfs(xd.first,u); } } int main() { ios_base::sync_with_stdio(false); cin.tie(NULL); cout.tie(NULL); cin>>n>>m>>k; for(int i=1;i<=m;++i) { cin>>old[i].u>>old[i].v>>old[i].w; old[i].idx=i; } sort(old+1,old+m+1); for(int i=0;i<k;++i) { cin>>nw[i].u; cin>>nw[i].v; nw[i].idx=m+i+1; } for(int i=1;i<=n;++i)cin>>pari[i]; dsu1.init(); for(int i=0;i<k;++i) { int p=dsu1.getRoot(nw[i].u); int q=dsu1.getRoot(nw[i].v); if(p!=q)dsu1.join(p,q); } for(int i=1;i<=m;++i) { int p=dsu1.getRoot(old[i].u); int q=dsu1.getRoot(old[i].v); if(p!=q) { alwaysin[i]=1; dsu1.join(p,q); } } dsu1.init(); for(int i=1;i<=m;++i) { int p=dsu1.getRoot(old[i].u); int q=dsu1.getRoot(old[i].v); if(p!=q) { if(alwaysin[i]!=1)sometimes_in.push_back(old[i]); dsu1.join(p,q); } } dsu1.init(); for(int i=1;i<=m;++i) { if(alwaysin[i]) { dsu1.join(dsu1.getRoot(old[i].u),dsu1.getRoot(old[i].v)); } } for(int i=1;i<=n;++i) { if(dsu1.getRoot(i)==i)comps.push_back(i); } for(int i=0;i<k;++i) { nw[i].u=dsu1.getRoot(nw[i].u); nw[i].v=dsu1.getRoot(nw[i].v); } for(int i=0;i<sometimes_in.size();++i) { sometimes_in[i].u=dsu1.getRoot(sometimes_in[i].u); sometimes_in[i].v=dsu1.getRoot(sometimes_in[i].v); } int root=dsu1.getRoot(1); for(int i=0;i<comps.size();++i) { pari1[comps[i]]=dsu1.h[comps[i]]; } for(int mask=1;mask<(1<<k);++mask) { for(int i=0;i<comps.size();++i) { dsu2.par[comps[i]]=comps[i]; dsu2.sz[comps[i]]=1; dsu2.h[comps[i]]=dsu1.h[comps[i]]; } bool ne_stava_brat=0; for(int i=0;i<comps.size();++i)g[comps[i]].clear(); for(int i=0;i<k;++i) { if(mask&(1<<i)) { int p=dsu2.getRoot(nw[i].u); int q=dsu2.getRoot(nw[i].v); if(p==q) { ne_stava_brat=1; break; } dsu2.join(p,q); g[nw[i].u].push_back({nw[i].v,nw[i].idx}); g[nw[i].v].push_back({nw[i].u,nw[i].idx}); } } if(ne_stava_brat)continue; vector<edge>ivan; for(int i=0;i<sometimes_in.size();++i) { int p=dsu2.getRoot(sometimes_in[i].u); int q=dsu2.getRoot(sometimes_in[i].v); if(p==q){ivan.push_back(sometimes_in[i]);continue;} dsu2.join(p,q); g[sometimes_in[i].u].push_back({sometimes_in[i].v,i}); g[sometimes_in[i].v].push_back({sometimes_in[i].u,i}); //ivan.push_back(sometimes_in[i]); } dfs1(root,0); for(auto xd:ivan)takovai(xd); sum=0; dfs(root,0); ans=max(ans,sum); } cout<<ans<<endl; return 0; } /* 5 5 1 3 5 2 1 2 3 2 3 5 2 4 4 4 3 6 1 3 10 20 30 40 50 */

Compilation message (stderr)

toll.cpp: In function 'int main()':
toll.cpp:180:15: warning: comparison of integer expressions of different signedness: 'int' and 'std::vector<edge>::size_type' {aka 'long unsigned int'} [-Wsign-compare]
  180 |  for(int i=0;i<sometimes_in.size();++i)
      |              ~^~~~~~~~~~~~~~~~~~~~
toll.cpp:186:15: warning: comparison of integer expressions of different signedness: 'int' and 'std::vector<int>::size_type' {aka 'long unsigned int'} [-Wsign-compare]
  186 |  for(int i=0;i<comps.size();++i)
      |              ~^~~~~~~~~~~~~
toll.cpp:192:16: warning: comparison of integer expressions of different signedness: 'int' and 'std::vector<int>::size_type' {aka 'long unsigned int'} [-Wsign-compare]
  192 |   for(int i=0;i<comps.size();++i)
      |               ~^~~~~~~~~~~~~
toll.cpp:199:16: warning: comparison of integer expressions of different signedness: 'int' and 'std::vector<int>::size_type' {aka 'long unsigned int'} [-Wsign-compare]
  199 |   for(int i=0;i<comps.size();++i)g[comps[i]].clear();
      |               ~^~~~~~~~~~~~~
toll.cpp:218:16: warning: comparison of integer expressions of different signedness: 'int' and 'std::vector<edge>::size_type' {aka 'long unsigned int'} [-Wsign-compare]
  218 |   for(int i=0;i<sometimes_in.size();++i)
      |               ~^~~~~~~~~~~~~~~~~~~~
#Verdict Execution timeMemoryGrader output
Fetching results...
#Verdict Execution timeMemoryGrader output
Fetching results...
#Verdict Execution timeMemoryGrader output
Fetching results...
#Verdict Execution timeMemoryGrader output
Fetching results...
#Verdict Execution timeMemoryGrader output
Fetching results...