NOIP系列

NOIP2015运输计划

 

唉

真是

这题

卡死我了

tarjan离线lca复杂度O（n）

最后各种卡常，多交几遍才A（洛谷104ms）

%%%zk学长609ms

注意二分的时候左边界要定成0

根据题意，显然先想到二分，然后求交集，在集合里找有没有边删掉之后使得所有路长度<mid

#include<bits/stdc++.h>
using namespace std;
char *TT,*mo,but[(1<<15)+2];
#define getchar() ((TT==mo&&(mo=(TT=but)+fread(but,1,1<<15,stdin)),TT==mo)?0:*TT++)
inline int read(){
    int x=0,c=0,f=1;
    for(;c<'0'||c>'9';c=getchar())f=c!='-';
    for(;c>='0'&&c<='9';c=getchar())x=x*10+c-'0';
    return f?x:-x;
}
int root;
struct Q{
    int a,b,lca,len;
}q[300010];
int head1[300010];
bool comp(Q A,Q B){
    return A.len>B.len;
}
struct T1{
    int next,to,num;
}e1[600010];
int cnt1;
void add1(int u,int v,int c){
    e1[++cnt1].to=v;e1[cnt1].next=head1[u];e1[cnt1].num=c;head1[u]=cnt1;
}
struct T{
    int next,to,cal;
}e[600010];
int cnt;
int head[300010];
int deep[300010];
int p[300010];
void add(int u,int v,int c){
    e[++cnt].to=v;e[cnt].next=head[u];e[cnt].cal=c;head[u]=cnt;
}
int sum[300010];
int n,m;
int fa[300010];
int find(int x){
    return !(fa[x]^x)?x:fa[x]=find(fa[x]);
}
void un(int x,int y){
    x=find(x),y=find(y);
    if(x^y)fa[x]=y;
}
bitset<300010>vis;
void init(int u){
    vis[u]=1;
    for(int i=head[u];i;i=e[i].next){
        if(vis[e[i].to])continue;
        deep[e[i].to]=deep[u]+e[i].cal;
        p[e[i].to]=e[i].cal;
        init(e[i].to);
        un(e[i].to,u);
    }
    for(int i=head1[u];i;i=e1[i].next){
        int v=e1[i].to;
        int ms=e1[i].num;
        if(vis[v]!=0){
            q[ms].lca=find(v);
        }
    } 
}
void dfs(int faf,int u){
    for(int i=head[u];i;i=e[i].next){
        if(!(e[i].to^faf))continue;
        dfs(u,e[i].to);
        sum[u]+=sum[e[i].to];
    }
}
bool check(int mid){
    memset(sum,0,sizeof(sum));
    int ret=0;
    register int i,j;
    for(i=1;i<=m;i++){
        if(q[i].len<=mid)break;
        sum[q[i].a]++;
        sum[q[i].b]++;
        sum[q[i].lca]-=2;
        ret=max(ret,q[i].len-mid);
    }
    i--;
    dfs(0,root);
    for(j=1;j<=n;j++){
        if(p[j]>=ret&&sum[j]==i)return 1;
    }
    return 0;
}
int main(){
    srand(time(NULL));
    n=read(),m=read();
    for(int i=1;i<=n;i++)fa[i]=i;
    int x,y,z;
    for(int i=1;i<n;i++){
        x=read(),y=read(),z=read();
        add(x,y,z);
        add(y,x,z);
    }
    root=rand()%n+1;
    deep[root]=0;
    for(int i=1;i<=m;i++){
        q[i].a=read(),q[i].b=read();
        add1(q[i].a,q[i].b,i);
        add1(q[i].b,q[i].a,i);
    }
    init(root);
    for(int i=1;i<=m;i++){
        q[i].len=deep[q[i].a]+deep[q[i].b]-2*deep[q[i].lca];
    }
    sort(q+1,q+m+1,comp);
    int l=0,r=300000000;
    while(r>l){
        int mid=l+r>>1;
        //cout<<mid<<endl;
        if(check(mid))r=mid;
        else l=mid+1;
    }
    cout<<l;
    return 0;
}