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+;

 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*],son[N],d[N],id[N],ne,nb,n,m,cnt;

 int read(){
     static int f,x;static char ch;
     x=f=;ch=getchar();
     while(!isdigit(ch)){f=(ch=='-');ch=getchar();}
     while(isdigit(ch)){x=x*+ch-'';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]=;
     for(int i=head[k];i;i=e[i].nxt){
         int v=e[i].to;
         if(v==fa)continue;
         d[v]=d[k]+;
         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+,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+,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(,,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(,,cnt,id[b]+,id[a]));
     return ret;
 }

 int main(){
     scanf("%d%d",&n,&m);
     foru(i,,m){re(a[i].f);re(a[i].t);re(a[i].w);F[i]=i;}
     sort(a+,a++m,cmp);
     foru(i,,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-;
     dfs(,);
     build(,);
     foru(i,,cnt)upd(,,cnt,i,inf);
     foru(i,,nb){
         if(d[b[i].f]<d[b[i].t])swap(b[i].f,b[i].t);
         upd(,,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\n",find(u,v));
     }
     return ;
 }