洛谷P2680 运输计划

https://www.luogu.org/problem/P2680
题意：給定一棵树，和m个路径起点终点，即m个运输任务，把任意一条边的边权变为0，使得m个任务同时开始，最大完成时间最小。
思路：要删的那条边或者说边权变为0的那条边，一定在最长路径上。
但是这个并没有什么用处。
要最大值最小，那么二分答案，转化为判定性问题，即判断能不能删一条边使得所有任务完成时间都小于等于mid。
check方法：枚举所有路径，只考虑长度大于mid的，如果存在一条【所有考虑的路径都经过】并且【所有考虑的路径减去这个边的权值都小于等于mid】的边，那么mid就是合法的，否则不合法。
【所有考虑的路径减去这个边的权值都小于等于mid】等价于【最长的路径减去这个边的权值小于等于mid】
怎么知道是不是所有考虑的路径都经过呢？用树上差分维护。
复杂度$O(log(n\ast t_{max})\ast (n+m))$
有一个点卡常数卡的厉害，需要快读+链式前向星才卡过。用刘汝佳版vector实现的邻接表过不了。

#include<bits/stdc++.h>
using namespace std;
const int maxn=300000+100;

int n,m,max_sub,max_path,num;
int fa[maxn][20],deep[maxn],w[maxn],d[maxn],dist[maxn];
int head[maxn],tot;
struct Edge{
    int to,dist,next;
}e[maxn*2];
struct Query{
    int u,v,lca,path;
}q[maxn];

void Read(int &x)
{
	int f=1;x=0;char s=getchar();
	while(s>'9' or s<'0'){if(s=='-')f=-1;s=getchar();}
	while(s>='0' and s<='9'){x=x*10+s-'0';s=getchar();}
	x*=f;
}

void AddEdge(int f,int t,int d)
{
    e[++tot]=(Edge){t,d,head[f]};
    head[f]=tot;
    e[++tot]=(Edge){f,d,head[t]};
    head[t]=tot;
}

void dfs(int u,int f)
{
    for(int i=head[u];i;i=e[i].next)
    {
        int v=e[i].to;
        if(v==f)continue;
        fa[v][0]=u;
        deep[v]=deep[u]+1;
        w[v]=w[u]+e[i].dist;
        dist[v]=e[i].dist;
        dfs(v,u);
    }
}

void lca_init()
{
    for(int i=1;i<=n;i++)for(int j=1;(1<<j)<n;j++)fa[i][j]=-1;
    for(int j=1;(1<<j)<n;j++)
        for(int i=1;i<=n;i++)
            if(fa[i][j-1]!=-1)
            {
                int a=fa[i][j-1];
                fa[i][j]=fa[a][j-1];
            }
}

int lca(int p,int q)
{
    int log;
    if(deep[p]<deep[q])swap(p,q);
    for(log=1;(1<<log)<=deep[p];log++);
    log--;
    int ans=-(1<<30);
    for(int i=log;i>=0;i--)
        if(deep[p]-(1<<i)>=deep[q])
        {
            p=fa[p][i];
        }
    if(p==q)return p; //lca为p或q
    for(int i=log;i>=0;i--)
        if(fa[p][i]!=-1 && fa[p][i]!=fa[q][i])
        {
            p=fa[p][i];
            q=fa[q][i];
        }
    return fa[p][0];//lca为fa[p]或fa[q]
}

int dfs2(int u,int f)
{   
    int cnt=d[u];
    for(int i=head[u];i;i=e[i].next)
    {
        int v=e[i].to;
        if(v==f)continue;
        cnt+=dfs2(v,u);
    }
    if(cnt==num)max_sub=max(max_sub,dist[u]);
    return cnt;
}

bool check(int mid)
{
    num=0;
    memset(d,0,sizeof(d));
    for(int i=1;i<=m;i++)
    {
        if(q[i].path>mid)
        {
            int x=q[i].lca;
            num++;
            d[q[i].u]++;d[q[i].v]++;
            d[x]-=2;
        }
    }  
    if(num==0)return true; 
    max_sub=0;
    dfs2(1,-1);
    if(max_path-max_sub<=mid)return true;
    else return false;
}

int main()
{
    //freopen("input.in","r",stdin);
    Read(n);Read(m);
    int a,b,c;
    for(int i=1;i<n;i++)
    {
        Read(a);Read(b);Read(c);
        AddEdge(a,b,c);
    }
    fa[1][0]=-1;
    deep[1]=0;
    w[1]=0;
    dfs(1,-1);
    lca_init();
    for(int i=1;i<=m;i++)
    {
        Read(q[i].u);Read(q[i].v);
        q[i].lca=lca(q[i].u,q[i].v);
        q[i].path=w[q[i].u]+w[q[i].v]-2*w[q[i].lca];
        max_path=max(max_path,q[i].path);
    }
    int l=0,r=max_path,mid,ans;      
    while(l<=r)
    {
        mid=(l+r)/2;
        if(check(mid))ans=mid,r=mid-1;
        else l=mid+1;                
    }
    cout<<ans<<endl;
    return 0;
}