【傳送門：BZOJ4326】

簡要題意：

　　給出一棵n個點的樹，每條邊有邊權，給出m條路徑

　　在可以將一條邊的邊權變成0的情況下，求出m條路徑的最大值最小

題解：

　　樹上差分+二分

　　首先把原來的圖構建出來，然後求出原圖的m條路徑的長度

　　然後二分答案，如果有路徑的長度大於二分的答案，那麽說明這條路徑是需要把一條邊的邊權變為0的

　　那麽如果有多條路徑的長度大於二分的答案，那麽說明這條邊必須在所有的路徑當中

　　用樹上差分來處理這些路徑上的邊出現的次數

　　然後枚舉每條邊，如果當前邊出現的次數為路徑的個數並且這些路徑中的最長路徑-這條邊的權值<=二分答案的話，那麽說明當前的二分答案是正確的

參考代碼：

#include<cstdio>
#include<cstring>
#include<cstdlib>
#include<algorithm>
#include<cmath>
using namespace std;
struct node
{
    int x,y,d,next;
}a[610000];int len,last[310000];
void ins(int x,int y,int d)
{
    len++;
    a[len].x=x;a[len].y=y;a[len].d=d;
    a[len].next
 
=last[x];last[x]=len;
}
int bin[21],f[310000][21],dep[310000],top[310000];
void dfs(int x)
{
    for(int i=1;bin[i]<=dep[x];i++) f[x][i]=f[f[x][i-1]][i-1];
    for(int k=last[x];k;k=a[k].next)
    {
        int y=a[k].y;
        if(y!=f[x][0])
        {
            dep[y]=dep[x]+1;
            f[y][
 
0]=x;
            top[y]=top[x]+a[k].d;
            dfs(y);
        }
    }
}
int LCA(int x,int y)
{
    if(dep[x]<dep[y]) swap(x,y);
    for(int i=20;i>=0;i--)
    {
        if(dep[x]-dep[y]>=bin[i])
        {
            x=f[x][i];
        }
    }
    if(x==y) return x;
    for(int i=20;i>=0;i--)
    {
        if(dep[x]>=bin[i]&&f[x][i]!=f[y][i])
        {
            x=f[x][i];y=f[y][i];
        }
    }
    return f[x][0];
}
int cf[310000];
int X[310000],Y[310000],dis[310000];
int fe[310000];
void qz(int x)
{
    for(int k=last[x];k;k=a[k].next)
    {
        int y=a[k].y;
        if(y!=f[x][0])
        {
            fe[y]=a[k].d;
            qz(y);
            cf[x]+=cf[y];
        }
    }
}
int lca[310000];
int main()
{
    bin[0]=1;for(int i=1;i<=20;i++) bin[i]=bin[i-1]<<1;
    int n,m;
    scanf("%d%d",&n,&m);
    for(int i=1;i<n;i++)
    {
        int x,y,d;
        scanf("%d%d%d",&x,&y,&d);
        ins(x,y,d);ins(y,x,d);
    }
    f[1][0]=0;dep[1]=1;top[1]=0;fe[1]=0;dfs(1);
    memset(cf,0,sizeof(cf));
    for(int i=1;i<=m;i++)
    {
        scanf("%d%d",&X[i],&Y[i]);
        lca[i]=LCA(X[i],Y[i]);
        dis[i]=top[X[i]]+top[Y[i]]-2*top[lca[i]];
    }
    int l=0,r=1<<31-1,ans;
    while(l<=r)
    {
        int mid=(l+r)/2;bool bk=false;
        int tot=0,mmax=0;
        memset(cf,0,sizeof(cf));
        for(int i=1;i<=m;i++)
        {
            if(dis[i]>mid)
            {
                tot++;mmax=max(mmax,dis[i]);
                cf[X[i]]++;cf[Y[i]]++;cf[lca[i]]-=2;
            }
        }
        qz(1);
        for(int i=1;i<=n;i++)
        {
            if(cf[i]==tot&&mmax-fe[i]<=mid)
            {
                bk=true;break;
            }
        }
        if(bk==true)
        {
            ans=mid;
            r=mid-1;
        }
        else l=mid+1;
    }
    printf("%d\n",ans);
    return 0;
}

BZOJ4326: NOIP2015 運輸計劃