NOIP2015 day2t3 运输计划

统计路径经过次数的差分+二分答案+倍增lca

---

BZOJ 4326（10080ms）
codevs 4632（最后一个点过不了）
写了俩小时还是没法过codevs最后一个点感觉自己萌萌哒，主函数里第一个乱糟糟的for循环其实就是建边，只不过被me拆掉写了······

但是CSDN id：_傲寒 的代码过了，可以去看一看他的
http://blog.csdn.net/qq_30401759/article/details/49777821

---

贴代码：

#include<cstdio>#include<cstring>#include<algorithm>using namespace std;int n,m,head[300001],tw=0;struct path{    int pre,to,len;}w[600001];int anc[20][300001],lca[300001],dep[300001],maxdist=-1,cnt,now,mid;int st[300001],ed[300001],pass[300001],dis[300001],dist[300005],ok;void build_tree(int u,int f){    anc[0][u]=f;    dep[u]=dep[f]+1;    for(int i=head[u];i;i=w[i].pre)        if(w[i].to!=f){            dis[w[i].to]=w[i].len+dis[u];            build_tree(w[i].to,u);        }}void getst(){    int i,j;    for(j=1;j<=19;j++)        for(i=1;i<=n;i++)            anc[j][i]=anc[j-1][anc[j-1][i]];}int getlca(int u,int v){    if( dep[u] < dep[v] ) swap(u,v);    int t = dep[u] - dep[v];    for( int p = 0; t; t>>=1,p++ )        if( t & 1 ) u = anc[p][u];    if( u == v ) return u;    for( int p = 19; p >= 0; p-- )         if( anc[p][u] != anc[p][v] )            u = anc[p][u], v = anc[p][v];    return anc[0][u];}void updata(int u,int f){    if(ok==1)   return;    int i;    for(i=head[u];i;i=w[i].pre)        if(w[i].to!=f){            updata(w[i].to,u);            pass[u]+=pass[w[i].to];        }    if(pass[u]==cnt&&maxdist-dis[u]+dis[anc[0][u]]<mid){        ok=1;now=u;    }}int solve(int L,int R){    if(L>R) return R;    memset(pass,0,sizeof(pass));    mid=(L+R)/2;cnt=0;    for(register int i=1;i<=m;i++){        if(dist[i]>=mid){            cnt++;            pass[st[i]]+=1;            pass[ed[i]]+=1;            pass[lca[i]]+=-2;        }    }    ok=0;    updata(1,0);    if(ok==1) return solve(L,mid-1);    return solve(mid+1,R);}char ch;int data;inline int read(){    data=0;ch=0;    while(ch<'0' || ch>'9') ch=getchar();    while(ch>='0' && ch<='9') data=data*10+ch-'0',ch=getchar();    return data;}int main(){    int t1;    scanf("%d%d",&n,&m);    for(register int i=1;i<n;i++){        t1=read();        w[++tw].pre=head[t1];        w[tw].to=read();        w[tw].len=read();        head[t1]=tw;        w[++tw].pre=head[w[tw-1].to];        w[tw].to=t1;        w[tw].len=w[tw-1].len;        head[w[tw-1].to]=tw;    }    build_tree(1,0);    getst();    for(register int i=1;i<=m;i++){        st[i]=read();ed[i]=read();        lca[i]=getlca(st[i],ed[i]);        dist[i]=dis[st[i]]+dis[ed[i]]-2*dis[lca[i]];        maxdist=max(maxdist,dist[i]);    }    printf("%d",solve(1,maxdist));}