| # | Time | Username | Problem | Language | Result | Execution time | Memory |
|---|---|---|---|---|---|---|---|
| 76283 | nxteru | Factories (JOI14_factories) | C++14 | 0 ms | 0 KiB |
This submission is migrated from previous version of oj.uz, which used different machine for grading. This submission may have different result if resubmitted.
#include <iostream>
#include <algorithm>
#include <vector>
#include <set>
#include <map>
#include <queue>
#include <stack>
#include <cstdio>
#include <cstring>
#include <string>
#include <math.h>
using namespace std;
typedef long long ll;
typedef double D;
typedef pair<int,ll> P;
#define M 1000000007
#define F first
#define S second
#define PB push_back
#define INF 100000000000000000
int n,q,k,par[500005][20],dp[500005],pr[500005],vs[500005],id[500005],seg[1<<20],qx[500005],qy[500005];
ll x[500005],y[500005],ans,dps[500005];
vector<P>g[500005];
vector<P>G[500005];
vector<int>cm;
vector<int>z;
void dfs(int v,int p,int d,ll da){
dps[v]=da;
vs[k]=v;
pr[v]=k++;
par[v][0]=p;
dp[v]=d;
for(int i=0;i<g[v].size();i++)if(g[v][i].F!=p)dfs(g[v][i].F,v,d+1,da+g[v][i].S);
}
int rmin(int a,int b){
if(dps[a]>dps[b])return b;
return a;
}
void up(int a,int b){
a+=(1<<19)-1;
seg[a]=b;
while(a>0){
a=(a-1)/2;
seg[a]=rmin(seg[a*2+1],seg[a*2+2]);
}
}
int que(int a,int b,int l,int r,int o){
if(r<a||b<l)return n;
if(a<=l&&r<=b)return seg[o];
return rmin(que(a,b,l,(l+r-1)/2,o*2+1),que(a,b,(l+r+1)/2,r,o*2+2));
}
int lca(int v,int u){
if(dp[v]>dp[u])swap(v,u);
for(int i=0;i<20;i++){
if((dp[u]-dp[v])>>i&1)u=par[u][i];
}
if(u==v)return v;
for(int i=19;i>=0;i--){
if(par[v][i]!=par[u][i]){
v=par[v][i];
u=par[u][i];
}
}
return par[v][0];
}
void s_dfs(int v,int lm,int rm){
if(lm+1<=v-1){
int l=que(lm+1,v-1,0,(1<<19)-1,0);
if(l!=n){
G[id[vs[v]]].PB(P(id[vs[l]],dps[vs[l]]-dps[vs[v]]));
G[id[vs[l]]].PB(P(id[vs[v]],dps[vs[l]]-dps[vs[v]]));
s_dfs(l,lm,v);
}
}
if(v+1<=rm-1){
int r=que(v+1,rm-1,0,(1<<19)-1,0);
if(r!=n){
G[id[vs[v]]].PB(P(id[vs[r]],dps[vs[r]]-dps[vs[v]]));
G[id[vs[r]]].PB(P(id[vs[v]],dps[vs[r]]-dps[vs[v]]));
s_dfs(r,v,rm);
}
}
}
void a_dfs(int v,int p){
for(int i=0;i<G[v].size();i++){
int u=G[v][i].F;
ll c=G[v][i].S;
if(u!=p){
a_dfs(u,v);
x[v]=min(x[v],x[u]+c);
y[v]=min(y[v],y[u]+c);
}
}
ans=min(ans,x[v]+y[v]);
}
int main(void){
scanf("%d%d",&n,&q);
for(int i=0;i<n;i++)id[i]=-1;
for(int i=0;i<1<<20;i++)seg[i]=n;
for(int i=0;i<n-1;i++){
int a,b;
ll c;
scanf("%d%d%lld",&a,&b,&c);
g[a].PB(P(b,c));
g[b].PB(P(a,c));
}
dfs(0,-1,0,0);
dps[n]=INF;
for(int i=0;i<19;i++){
for(int v=0;v<n;v++){
if(par[v][i]==-1)par[v][i+1]=-1;
else par[v][i+1]=par[par[v][i]][i];
}
}
while(q){
k=0;
cm.clear();
z.clear();
int s,t;
scanf("%d%d",&s,&t);
for(int i=0;i<s;i++){
scanf("%d",&qx[i]);
if(id[qx[i]]==-1){
cm.PB(pr[qx[i]]);
up(pr[qx[i]],pr[qx[i]]);
id[qx[i]]=k++;
}
}
for(int i=0;i<t;i++){
scanf("%d",&qy[i]);
if(id[qy[i]]==-1){
cm.PB(pr[qy[i]]);
up(pr[qy[i]],pr[qy[i]]);
id[qy[i]]=k++;
}
}
sort(cm.begin(),cm.end());
for(int i=1;i<cm.size();i++){
int a=lca(vs[cm[i]],vs[cm[i-1]]);
if(id[a]==-1){
id[a]=k++;
up(pr[a],pr[a]);
}
z.PB(a);
}
int r=que(cm[0],cm.back(),0,(1<<19)-1,0);
s_dfs(r,-1,n);
ans=INF;
for(int i=0;i<k;i++){
x[i]=INF;
y[i]=INF;
}
for(int i=0;i<s;i++)x[id[qx[i]]]=0;
for(int i=0;i<t;i++)y[id[qy[i]]]=0;
a_dfs(id[vs[r]],-1);
printf("%lld\n",ans);
for(int i=0;i<z.size();i++){
up(pr[z[i]],n);
id[z[i]]=-1;
}
for(int i=0;i<s;i++){
up(pr[qx[i]],n);
id[qx[i]]=-1;
}
for(int i=0;i<t;i++){
up(pr[qy[i]],n);
id[qy[i]]=-1;
}
for(int i=0;i<k;i++)G[i].clear();
q--;
}
}