疫情控制

H 国有 n 个城市，这 n 个城市用 n-1 条双向道路相互连通构成一棵树，1号城市是首都，也是树中的根节点。
H 国的首都爆发了一种危害性极高的传染病。
当局为了控制疫情，不让疫情扩散到边境城市（叶子节点所表示的城市），决定动用军队在一些城市建立检查点，使得从首都到边境城市的每一条路径上都至少有一个检查点，边境城市也可以建立检查点。
但要注意的是，首都是不能建立检查点的。
现在，在H国的一些城市中已经驻扎有军队，且一个城市可以驻扎多个军队。
军队总数为 m 支。
一支军队可以在有道路连接的城市间移动，并在除首都以外的任意一个城市建立检查点，且只能在一个城市建立检查点。
一支军队经过一条道路从一个城市移动到另一个城市所需要的时间等于道路的长度（单位：小时）。
请问：最少需要多少个小时才能控制疫情？
注意：不同的军队可以同时移动。

输入格式
第一行一个整数n，表示城市个数。
接下来的n-1行，每行3个整数，u、v、w，每两个整数之间用一个空格隔开，表示从城市u到城市v有一条长为w的道路，数据保证输入的是一棵树，且根节点编号为1。
接下来一行一个整数m，表示军队个数。
接下来一行m个整数，每两个整数之间用一个空格隔开，分别表示这m个军队所驻扎的城市的编号。

输出格式
共一行，包含一个整数，表示控制疫情所需要的最少时间。如果无法控制疫情则输出-1。

数据范围
2≤m≤n≤50000,
0<w<109

输入样例：
4
1 2 1
1 3 2
3 4 3
2
2 2
输出样例：
3

#pragma GCC optimize(2)

#include<bits/stdc++.h>

#define ll long long
using namespace std;
const int N = 5e4 + 10;
int head[N], ver[N * 2], edge[N * 2], Next[N * 2], tot;
int n, m;
int t, f[N][20], dist[N][20], d[N];
ll cnt, r, l, mid;
int troop[N];
bool blo[N];

int tim[N], ned[N], atot, btot, ctot;
pair<int, int> h[N];
bool need[N];

void add(int x, int y, int z) {
    ver[++tot] = y;
    edge[tot] = z;
    Next[tot] = head[x];
    head[x] = tot;
}

void read() {
    cin >> n;
    for (int i = 1; i <= n - 1; i++) {
        int x, y, z;
        scanf("%d%d%d", &x, &y, &z);
        r += z;
        add(x, y, z);
        add(y, x, z);
    }
    cin >> m;
    for (int i = 1; i <= m; i++)
        scanf("%d", &troop[i]);
}

void bfs() {
    t = (int) log2(n) + 1;
    queue<int> q;
    q.push(1);
    d[1] = 1;
    while (!q.empty()) {
        int x = q.front();
        q.pop();
        for (int i = head[x]; i; i = Next[i]) {
            int y = ver[i];
            if (d[y])continue;
            d[y] = d[x] + 1;
            dist[y][0] = edge[i];
            f[y][0] = x;
            for (int j = 1; j <= t; j++) {
                f[y][j] = f[f[y][j - 1]][j - 1];
                dist[y][j] = dist[y][j - 1] + dist[f[y][j - 1]][j - 1];
            }
            q.push(y);
        }
    }
}

bool dfs(int x) {
    bool is_leaf = true;
    if (blo[x])return true;
    for (int i = head[x]; i; i = Next[i]) {
        int y = ver[i];
        if (y == f[x][0])continue;
        is_leaf = false;
        if (!dfs(y))return false;
    }
    return !is_leaf;
}

void init() {
    for (int i = 1; i <= n; i++) {
        need[i] = blo[i] = false;
    }
    atot = btot = ctot = 0;
}

bool check(ll lim) {
    init();
    for (int i = 1; i <= m; i++) {
        int x = troop[i];
        cnt = 0;
        for (int j = t; j >= 0; j--)
            if (f[x][j] > 1 && lim - dist[x][j] - cnt >= 0) {
                cnt += dist[x][j];
                x = f[x][j];
            }
        if (f[x][0] == 1 && lim - dist[x][0] - cnt > 0)
            h[++ctot] = make_pair(lim - dist[x][0] - cnt, x);
        else blo[x] = true;
    }
    sort(h + 1, h + ctot + 1);
    for (int i = head[1]; i; i = Next[i])
        if (!dfs(ver[i]))
            need[ver[i]] = true;
    for (int i = 1; i <= ctot; i++)
        if (need[h[i].second] && h[i].first < dist[h[i].second][0])
            need[h[i].second] = false;
        else
            tim[++atot] = h[i].first;
    for (int i = head[1]; i; i = Next[i])
        if (need[ver[i]])ned[++btot] = edge[i];
    if (atot < btot)
        return false;
    sort(tim + 1, tim + atot + 1);
    sort(ned + 1, ned + btot + 1);
    int i = 1, j = 1;
    while (i <= atot && j <= btot)
        tim[i] >= ned[j] ? (i++, j++) : (i++);
    return j > btot;

}

void solve() {
    bool flag = false;
    while (l < r) {
        mid = (l + r) >> 1;
        check(mid) ? (r = mid, flag = true) : (l = mid + 1);
    }
    cout << (!flag ? -1 : l);
}

int main() {
    read();
    bfs();
    solve();
    return 0;
}