【NOIP2015 Day2】运输计划

问题描述

公元 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行包含两个正整数ui和vi，表示第j个运输计划是从uj号星球飞往vj号星球。

输出格式

/输出文件名为transport.out。/
共1行，包含1个整数，表示小P的物流公司完成阶段性工作所需要的最短时间。

样例输入：

6 3
1 2 3
1 6 4
3 1 7
4 3 6
3 5 5
3 6
2 5
4 5

样例输出：

11

说明

输入输出样例1说明：
将第1条航道改造成虫洞，则三个计划耗时分别为11、12、11，故需要花费的时间为12；
将第2条航道改造成虫洞，则三个计划耗时分别为7、15、11，故需要花费的时间为15；
将第3条航道改造成虫洞，则三个计划耗时分别为4、8、11，故需要花费的时间为11；
将第4条航道改造成虫洞，则三个计划耗时分别为11、15、5，故需要花费的时间为15；
将第5条航道改造成虫洞，则三个计划耗时分别为11、10、6，故需要花费的时间为11。
故将第3条或第5条航道改造成虫洞均可使得阶段性工作的耗时最短，需要花费的时间为11。

数据范围

---

题解

由题意，使一条边长度变为零，要使最长距离尽可能短。
暴力枚举显然不行。求“最大的最小”，可用二分答案。
每次二分一个长度x，判断是否能让所有路径长度不超过x。
如何判断？
枚举所有长度大于x的路径，将路径上每条边权值加一。假设开始有k 条路径长度大于x， 显然只有将权值为k 的一条边改为虫洞，才有可能让所有路径长度不超过x 。
至于如何快速修改边权，可用树上差分。
具体见代码。

---

代码

#include <cstdio>
#include <iostream>
#include <cmath>
#include <algorithm>
#include <cstring>
using namespace std;
const int Q=666666,P=333333;
int nn[Q],e[Q],v[Q],last[P]={0},dis[P][21],n,s,f[P][21],dep[P],tot=0,bbq[P],m;
bool ch[P]={0};
struct dt{
    int a,b,l,p;
}q[P];
void add(int x,int y,int z)
{
    v[++tot]=z;
    e[tot]=y;
    nn[tot]=last[x];
    last[x]=tot;
}   
void search(int x)
{
    int y,t;
    ch[x]=true;
    for(t=last[x];t;t=nn[t])
    {
        y=e[t];
        if(ch[y])continue;
        dep[y]=dep[x]+1;
        dis[y][0]=v[t];
        f[y][0]=x;
        search(y);
    }
}
void lca(int x,int y,int &temp,int &dian)
{
    int i,det;temp=0;
    if(dep[x]<dep[y])swap(x,y);
    det=dep[x]-dep[y];
    for(i=s;i>=0;--i)
        if((det>>i)&1)temp+=dis[x][i],x=f[x][i];
    if(x==y){dian=x;
        return;
    }
    for(i=s;i>=0;--i)
        if(f[x][i]!=f[y][i])
            temp+=dis[x][i]+dis[y][i],x=f[x][i],y=f[y][i];
    temp+=dis[x][0]+dis[y][0];
    dian=f[x][0];
    return;
}
void Mod_Root(int x)
{
    ch[x]=true;
    int t,y;
    for(t=last[x];t;t=nn[t])
    {
        y=e[t];
        if(ch[y])continue;
        Mod_Root(y);
        bbq[x]+=bbq[y];
    }
}
bool hac(int x)
{
    int i,temp=0,ha=0;
    memset(ch,0,sizeof(ch));
    memset(bbq,0,sizeof(bbq));
    for(i=1;i<=m;i++)
    {
        if(q[i].l>x)++bbq[q[i].a],++bbq[q[i].b],bbq[q[i].p]-=2,temp++,ha=max(ha,q[i].l-x);
    }
    Mod_Root(1);
    for(i=1;i<=n;i++)
        if(bbq[i]==temp&&dis[i][0]>=ha)return true;
    return false;
}
int find(int l,int r)
{
    while(l<=r)
    {
        int mid=(l+r)>>1;
        if(hac(mid))r=mid-1;
        else l=mid+1;
    }
    return l;
}
int tmain()
{
    int i,x,y,j;
    scanf("%d%d",&n,&m);
    s=int(floor(log2(n)))+1;
    for(i=1;i<n;i++)
    {
        scanf("%d%d%d",&x,&y,&j);
        add(x,y,j),add(y,x,j);
    }
    dep[1]=1,f[1][0]=1;
    search(1);
    for(j=1;j<=s;j++)
        for(i=1;i<=n;i++)
            f[i][j]=f[f[i][j-1]][j-1],dis[i][j]=dis[i][j-1]+dis[f[i][j-1]][j-1];
    for(i=1;i<=m;i++)
    {
        scanf("%d%d",&q[i].a,&q[i].b);
        lca(q[i].a,q[i].b,q[i].l,q[i].p);
    }
    printf("%d",find(0,300000000));
    return 0;
}