【Luogu】P2680运输计划（树上差分+二分）

　　题目链接

　　总体思路……怎么说呢……是个暴力吧……

　　首先用倍增预处理出每条路径的长度。　　然后按长度把路径排序。

　　然后二分答案。对于当前答案mid检验，怎么检验呢？

　　首先差分把所有长度比mid大的链上除了LCA之外的所有点权+1。dfs求出每个点的点权，顺便记下有多少点是被所有路径经过的。对于这些点，如果有一个点使得它到他的父亲那条边被删掉后，长度最长的路径都能卡过mid，那么返回true.

　　那么检查所有点发现长度最长的路径卡不过mid，那么返回false。

　　然后这题搞完了。

　　

#include<cstdio>
#include<cstdlib>
#include<cctype>
#include<cstring>
#include<algorithm>

using namespace std;
int n,m;

inline int read(){
    int num=,f=;
    char ch=getchar();
    while(!isdigit(ch)){
        if(ch=='-')    f=-;
        ch=getchar();
    }
    while(isdigit(ch)){
        num=num*+ch-'';
        ch=getchar();
    }
    return num*f;
}

struct Edge{
    int next,to,dis;
}edge[];
int head[],num;
inline void add(int from,int to,int dis){
    edge[++num]=(Edge){head[from],to,dis};
    head[from]=num;
}
struct Node{
    int from,to,dis,lca;
    bool operator <(const Node &a)const{
        return dis>a.dis;
    }
}q[];
int s[][],d[][];
int deep[];
int vis[],e[];
int pos[],tot;
int Max,ans;

void find(int x,int fa){
    deep[x]=deep[fa]+;
    for(int i=head[x];i;i=edge[i].next){
        int to=edge[i].to;
        if(to==fa)    continue;
        s[to][]=x;d[to][]=edge[i].dis;find(to,x);
    }
}

void dfs(int x,int fa,int size){
    e[x]=vis[x];
    for(int i=head[x];i;i=edge[i].next){
        int to=edge[i].to;
        if(to==fa)    continue;
        dfs(to,x,size);
        e[x]+=e[to];
    }
    if(e[x]==size)    pos[++tot]=x;
}

bool check(int limit){
    memset(vis,,sizeof(vis));
    int size=;tot=;
    for(int i=;i<=m;++i){
        if(q[i].dis<=limit)    break;
        vis[q[i].from]++;vis[q[i].to]++;vis[q[i].lca]-=;
        size++;
    }
    dfs(,,size);
    for(int i=;i<=tot;++i)
        if(q[].dis-d[pos[i]][]<=limit)    return ;
    return ;
}

int main(){
    n=read(),m=read();
    for(int i=;i<n;++i){
        int from=read(),to=read(),dis=read();
        add(from,to,dis);
        add(to,from,dis);
    }
    for(int i=;i<=m;++i)    q[i]=(Node){read(),read(),,};
    find(,);
    for(int j=;j<;++j)
        for(int i=;i<=n;++i){
            s[i][j]=s[s[i][j-]][j-];
            d[i][j]=d[i][j-]+d[s[i][j-]][j-];
        }
    for(int i=;i<=m;++i){
        int from=q[i].from,to=q[i].to;    int &dis=q[i].dis,&lca=q[i].lca;
        if(deep[from]<deep[to])    swap(from,to);
        int f=deep[from]-deep[to];
        for(int j=;(<<j)<=f;++j)
            if(f&(<<j)){
                dis+=d[from][j];
                from=s[from][j];
            }
        if(from==to){
            lca=from;
            Max=max(Max,dis);
            continue;
        }
        for(int j=;j>=;--j){
            if(s[from][j]==s[to][j])    continue;
            dis+=d[from][j]+d[to][j];
            from=s[from][j];to=s[to][j];
        }
        dis+=d[from][]+d[to][];
        lca=s[from][];
        Max=max(Max,dis);
    }
    sort(q+,q+m+);
    int l=,r=Max;
    while(l<=r){
        int mid=(l+r)>>;
        if(check(mid)){    ans=mid;r=mid-;    }
        else        l=mid+;
    }
    printf("%d",ans);
    return ;
}