NOIP2013D1T3货车运输 (生成树+树链剖分)

给出一个图，询问图上两点间路径上最小边权的最大值。

先跑一次最大生成树。

树剖维护路径最小边权。

树剖又双叒叕写挂了。

#include<cstring>
#include<cstdio>
#include<cmath>
#include<algorithm>
#include<cctype>
#define foru(i,x,y) for(int i=x;i<=y;i++)
#define mm(a) memset(a,0,sizeof(a))
#define ford(i,x,y) for(int i=x;i>=y;i--)
#define re(x) x=read()
using namespace std;
typedef long long LL;
typedef double db;
const int inf=1e9;
const int N=2e5+10;

struct dat{int f,t,w;}a[N],b[N];
struct edge{int to,nxt,w;}e[N];
int head[N],top[N],f[N],F[N],siz[N],t[N*20],son[N],d[N],id[N],ne,nb,n,m,cnt;

int read(){
    static int f,x;static char ch;
    x=f=0;ch=getchar();
    while(!isdigit(ch)){f=(ch=='-');ch=getchar();}
    while(isdigit(ch)){x=x*10+ch-'0';ch=getchar();}
    return f?-x:x;
}

bool cmp(dat a,dat b){
    return a.w>b.w;
}

int gf(int a){
    if(a==F[a])return a;
    return F[a]=gf(F[a]);
}

void add(int a,int b,int c){
    e[++ne]=(edge){b,head[a],c};head[a]=ne;
    e[++ne]=(edge){a,head[b],c};head[b]=ne;
}

void dfs(int k,int fa){
    f[k]=fa;siz[k]=1;
    for(int i=head[k];i;i=e[i].nxt){
        int v=e[i].to;
        if(v==fa)continue;
        d[v]=d[k]+1;
        dfs(v,k);
        siz[k]+=siz[v];
        if(siz[v]>siz[son[k]])son[k]=v;
    }
}

void build(int k,int tp){
    id[k]=++cnt;top[k]=tp;
    if(son[k])build(son[k],tp);
    for(int i=head[k];i;i=e[i].nxt){
        int v=e[i].to;
        if(v==f[k])continue;
        if(v!=son[k])build(v,v);
    }
}

#define mid ((L+R)>>1)
#define ls (k<<1)
#define rs (k<<1|1)

void upd(int k,int L,int R,int p,int x){
    if(p<L||p>R)return;
    if(L==R){t[k]=x;return;}
    upd(ls,L,mid,p,x);upd(rs,mid+1,R,p,x);
    t[k]=min(t[ls],t[rs]);
}

int quiry(int k,int L,int R,int l,int r){
    if(r<L||l>R)return inf;
    if(l<=L&&R<=r)return t[k];
    return min(quiry(ls,L,mid,l,r),quiry(rs,mid+1,R,l,r));
}

int find(int a,int b){
    int ret=inf;
    while(top[a]!=top[b]){
        if(d[top[a]]<d[top[b]])swap(a,b);
        //又写成了(d[a]<d[b]) 
        int tmp=quiry(1,1,cnt,id[top[a]],id[a]);
        ret=min(ret,tmp);
        a=f[top[a]];
    }
    if(d[a]<d[b])swap(a,b);
    if(a!=b)ret=min(ret,quiry(1,1,cnt,id[b]+1,id[a]));
    return ret;
}

int main(){
    scanf("%d%d",&n,&m);
    foru(i,1,m){re(a[i].f);re(a[i].t);re(a[i].w);F[i]=i;}
    sort(a+1,a+1+m,cmp);
    foru(i,1,m){
        int f1=gf(a[i].f),f2=gf(a[i].t);
        if(f1!=f2){
            F[f2]=f1;
            add(a[i].f,a[i].t,a[i].w);
            b[++nb]=a[i];
        }
    }
    mm(f);m=n-1;
    dfs(1,0);
    build(1,1);
    foru(i,1,cnt)upd(1,1,cnt,i,inf);
    foru(i,1,nb){
        if(d[b[i].f]<d[b[i].t])swap(b[i].f,b[i].t);
        upd(1,1,cnt,id[b[i].f],b[i].w);
    }
    int q=read(),u,v;
    while(q--){
        re(u);re(v);
        if(gf(u)!=gf(v)){puts("-1");continue;}
        printf("%d
",find(u,v));
    }
    return 0;
}