解析
二分的单调性较为明显,一路推导下去的性质都不算太难想,正解的思路还是不难想到的
但从头到尾都实现很考验思维的严密性和代码能力
然后我就双重被考验挂了qwq
第一交的时候一个地方把倍增的dis数组写成pl,判断封口也有问题…
但竟然有70
真实考试的时候绝对不能出打错数组这种错!提交前还是要谨慎一些!
另外那个判封口…只能说路漫漫其修远兮了…
我对于最后根下方儿子(暂且叫它次根好了)与军队贪心匹配的地方写的和题解不太一样
题解的实现是:爬到根下方的军队不考虑,然后把所有军队提出来一起贪心
我的做法是先给有军队的次根一个剩余最少的军队,再维护一个堆看是否需要用外面的军队把里面的军队替换出来
这两个实现我觉得差不多吧,谈不上孰优孰劣
然而我就写挂了qwq
代码
#include<bits/stdc++.h>
using namespace std;
#define ll long long
const int N=1e5+100;
ll read() {ll x=0,f=1;char c=getchar();while(!isdigit(c)) {if(c=='-')f=-1;c=getchar();}while(isdigit(c)) {x=x*10+(c^48);c=getchar();}return x*f;
}int n,m;
int x[N];
struct node{int to,nxt,w;
}p[N<<1];
int fi[N],ecnt;
inline void addline(int x,int y,int w){p[++ecnt]=(node){y,fi[x],w};fi[x]=ecnt;return;
}
ll dis[N][20];
int pl[N][20];
int que[N],num,jd[N];
void dfs(int x,int f){for(int k=1;k<=16;k++){pl[x][k]=pl[pl[x][k-1]][k-1];dis[x][k]=dis[x][k-1]+dis[pl[x][k-1]][k-1];}jd[x]=1;for(int i=fi[x];~i;i=p[i].nxt){int to=p[i].to;if(to==f) continue;jd[x]=0;if(x==1) que[++num]=to;pl[to][0]=x;dis[to][0]=p[i].w;dfs(to,x);}//printf("x=%d jd=%d\n",x,jd[x]);return;
}bool vis[N],tag[N];
vector<ll>v[N];
int solve(int x,int f){//printf("solve:x=%d vis=%d jd=%d\n",x,vis[x],jd[x]);if(vis[x]) return tag[x]=1;if(jd[x]) return tag[x]=0;tag[x]=1;for(int i=fi[x];~i;i=p[i].nxt){int to=p[i].to;if(to==f) continue;tag[x]&=solve(to,x);}//printf("x=%d tag=%d\n",x,tag[x]);return tag[x];
}
ll w[N],cnt,add[N],o;
priority_queue<ll>q;
struct aaa{ll ned,val;bool operator < (const aaa oth)const{return ned<oth.ned;}
}a[N];;
int numa;
bool cmp(ll a,ll b){return a>b;}
void init(){memset(vis,0,sizeof(vis));memset(tag,0,sizeof(tag));for(int i=1;i<=num;i++){v[que[i]].clear();v[que[i]].shrink_to_fit();}cnt=o=numa=0;while(!q.empty()) q.pop();
}
bool check(ll tot){init();for(int o=1;o<=m;o++){int pos=x[o];ll lft=tot;for(int k=16;k>=0;k--){if(dis[pos][k]>lft||pl[pos][k]==0||pl[pos][k]==1) continue;lft-=dis[pos][k];pos=pl[pos][k];} //printf("pos=%d lft=%lld\n",pos,lft);if(pl[pos][0]==1){v[pos].push_back(lft);}else vis[pos]=1;}if(solve(1,0)) return true;for(int i=1;i<=num;i++){int now=que[i];if(v[now].size()==0&&!tag[now]){w[++cnt]=dis[now][0];continue;}sort(v[now].begin(),v[now].end());if(!tag[now]){if(v[now][0]-dis[now][0]>0){//printf("now=%d v=%lld dis=%lld\n",now,v[now][0],dis[now][0]);a[++numa]=(aaa){dis[now][0],v[now][0]-dis[now][0]};}for(int j=1,tp=v[now].size();j<tp;j++){if(v[now][j]>dis[now][0]) add[++o]=v[now][j]-dis[now][0];}}else{for(int j=0,tp=v[now].size();j<tp;j++){if(v[now][j]>dis[now][0]) add[++o]=v[now][j]-dis[now][0];}}}if(o<cnt) return false;sort(a+1,a+1+numa);sort(add+1,add+1+o);sort(w+1,w+1+cnt,cmp);//printf("w: ");for(int i=1;i<=cnt;i++) printf("%lld ",w[i]);putchar('\n');//printf("add: ");for(int i=1;i<=o;i++) printf("%lld ",add[i]);putchar('\n');//printf("aaa: ");for(int i=1;i<=numa;i++) printf("(%lld %lld) ",a[i].ned,a[i].val);putchar('\n');int pos=0;for(int k=1;k<=o;k++){while(pos<numa&&add[k]>=a[pos+1].ned){pos++;q.push(a[pos].val);}if(q.empty()) continue;ll now=q.top();q.pop();if(add[k]<now) swap(add[k],now);q.push(now);}sort(add+1,add+1+o,cmp);//printf("add: ");for(int i=1;i<=o;i++) printf("%lld ",add[i]);putchar('\n');for(int i=1;i<=cnt;i++){if(w[i]>add[i]) return false;}return true;
}int main() {
#ifndef ONLINE_JUDGE//freopen("a.in","r",stdin);//freopen("a.out","w",stdout);
#endifmemset(fi,-1,sizeof(fi));ecnt=-1;n=read();for(int i=1;i<n;i++){int x=read(),y=read(),w=read();addline(x,y,w);addline(y,x,w);}m=read();for(int i=1;i<=m;i++) x[i]=read();dfs(1,0);//printf("check=%d\n",check(8));//return 0;ll st=0,ed=1e14;while(st<ed){ll mid=(st+ed)>>1;if(check(mid)) ed=mid;else st=mid+1;}printf("%lld\n",st<1e14?st:-1);return 0;
}
/*
4
1 2 10
1 3 1
4 1 5
3
3 4 4
*/