luogu2680 运输计划

http://www.elijahqi.win/archives/1371
题目背景
公元 2044 年,人类进入了宇宙纪元。

题目描述
L 国有 n 个星球,还有 n-1 条双向航道,每条航道建立在两个星球之间,这 n-1 条航道连通了 L 国的所有星球。

小 P 掌管一家物流公司,该公司有很多个运输计划,每个运输计划形如:有一艘物

流飞船需要从 ui 号星球沿最快的宇航路径飞行到 vi 号星球去。显然,飞船驶过一条航道 是需要时间的,对于航道 j,任意飞船驶过它所花费的时间为 tj,并且任意两艘飞船之 间不会产生任何干扰。

为了鼓励科技创新,L 国国王同意小 P 的物流公司参与 L 国的航道建设,即允许小 P 把某一条航道改造成虫洞,飞船驶过虫洞不消耗时间。

在虫洞的建设完成前小 P 的物流公司就预接了 m 个运输计划。在虫洞建设完成后, 这 m 个运输计划会同时开始,所有飞船一起出发。当这 m 个运输计划都完成时,小 P 的 物流公司的阶段性工作就完成了。

如果小 P 可以自由选择将哪一条航道改造成虫洞,试求出小 P 的物流公司完成阶段 性工作所需要的最短时间是多少?

输入输出格式
输入格式：

输入文件名为 transport.in。

第一行包括两个正整数 n、m,表示 L 国中星球的数量及小 P 公司预接的运输计划的数量,星球从 1 到 n 编号。

接下来 n-1 行描述航道的建设情况,其中第 i 行包含三个整数 ai, bi 和 ti,表示第

i 条双向航道修建在 ai 与 bi 两个星球之间,任意飞船驶过它所花费的时间为 ti。

接下来 m 行描述运输计划的情况,其中第 j 行包含两个正整数 uj 和 vj,表示第 j个 运输计划是从 uj 号星球飞往 vj 号星球。

输出格式：

输出 共1行,包含1个整数,表示小P的物流公司完成阶段性工作所需要的最短时间。

输入输出样例
输入样例#1： 复制

6 3
1 2 3
1 6 4
3 1 7
4 3 6
3 5 5
3 6
2 5
4 5
输出样例#1： 复制

11
说明
所有测试数据的范围和特点如下表所示

请注意常数因子带来的程序效率上的影响。

真tmd一个卡常题 卡了半天没卡过去

二分答案 然后每次对所有的线路扫描一遍 然后 求出我交集 然后用线路中最大值–交集中的边的最大值 看能否小于我二分的答案

#include<cstdio>
#define N 300010
#include<algorithm>
#include<cstring>
using namespace std;
inline char gc(){
    static char now[1<<16],*T,*S;
    if (T==S){T=(S=now)+fread(now,1,1<<16,stdin);if (T==S) return EOF;}
    return *S++;
}
inline int read(){
    int x=0;char ch=getchar();
    while (ch<'0'||ch>'9') ch=getchar();
    while (ch<='9'&&ch>='0'){x=x*10+ch-'0';ch=getchar();}
    return x;
}
struct node{
    int y,next,z;
}data[N<<1];
struct node1{
    int st,ed,value,lca;
}edge[N];
int num,n,m,h[N],Log[N],fa[N][20],sum[N][20],dep[N],flag[N];
void dfs(int x){
    for (int i=h[x];i;i=data[i].next){
        int y=data[i].y,z=data[i].z;
        if (fa[x][0]==y) continue;fa[y][0]=x;dep[y]=dep[x]+1;sum[y][0]=z;
        for (int j=1;j<=Log[dep[y]];++j){
            fa[y][j]=fa[fa[y][j-1]][j-1];sum[y][j]=sum[y][j-1]+sum[fa[y][j-1]][j-1];
        }dfs(y);
    }
}
inline int calc(int x,int y,int id){
    int sum1=0;
    if (dep[x]<dep[y]) swap(x,y);
    int dis=dep[x]-dep[y];
    for (int i=0;i<=Log[dis];++i){
        if ((1<<i)&dis) sum1+=sum[x][i],x=fa[x][i];
    }edge[id].lca=x;
    if (x==y)return sum1;
    for (int i=Log[dep[y]];i>=0;--i){
        if (fa[x][i]!=fa[y][i]){
            sum1+=sum[x][i]+sum[y][i];x=fa[x][i];y=fa[y][i];
        }
    }edge[id].lca=fa[x][0];
    sum1+=sum[x][0]+sum[y][0];return sum1;
}
void dfs1(int x){
    for (int i=h[x];i;i=data[i].next){
        int y=data[i].y;if (fa[x][0]==y) continue;
        dfs1(y);flag[x]+=flag[y];
    }
}
inline bool cmp(node1 a,node1 b){
    return a.value>b.value;
}
inline bool check(int x){
    memset(flag,0,sizeof(flag));
    int max_dis=0,k=0;
    for (int i=1;i<=m;++i){
        if (edge[i].value>x){
            flag[edge[i].st]++;flag[edge[i].ed]++;flag[edge[i].lca]-=2;
            max_dis=max(max_dis,edge[i].value);++k;
        }else break;
    }dfs1(1);int maxans=0;
    for (int i=2;i<=n;++i){
        if (flag[i]==k){
            maxans=max(maxans,sum[i][0]);
        }
    }return max_dis-maxans<=x;
}
int main(){
    //freopen("2680.in","r",stdin);
    n=read();m=read();
    for (int i=1;i<n;++i){
        int x=read(),y=read(),z=read();
        data[++num].y=y;data[num].next=h[x];data[num].z=z;h[x]=num;
        data[++num].y=x;data[num].next=h[y];data[num].z=z;h[y]=num;
    }Log[0]=-1;num=0;
    for (int i=1;i<=n;++i) Log[i]=Log[i>>1]+1;
    dfs(1);
    for (int i=1;i<=m;++i){
        int x=read(),y=read();
        edge[i].st=x;edge[i].ed=y;edge[i].value=calc(x,y,i);
    }sort(edge+1,edge+m+1,cmp);
    int l=0,r=edge[1].value;
    while (l<=r){
        int mid=l+r>>1;
        if (check(mid)) r=mid-1;else l=mid+1;   
    }printf("%d",l);
    return 0;
}