LUOGU 1967 货车运输 [noip 2013]

*

1967 货车运输

*
题目描述
AA 国有 n n 座城市，编号从 1 1 到 nn ，城市之间有 mm 条双向道路。每一条道路对车辆都有重量限制，简称限重。现在有 qq 辆货车在运输货物， 司机们想知道每辆车在不超过车辆限重的情况下，最多能运多重的货物。

输入输出格式
输入格式：
第一行有两个用一个空格隔开的整数 n,mn,m ，表示 AA 国有 nn 座城市和 mm 条道路。

接下来 mm 行每行 3 3 个整数 x, y, zx,y,z ，每两个整数之间用一个空格隔开，表示从 x x 号城市到 y y 号城市有一条限重为 zz 的道路。注意： xx 不等于 yy ，两座城市之间可能有多条道路 。

接下来一行有一个整数 q，表示有 q 辆货车需要运货。

接下来 q 行，每行两个整数 x、y，之间用一个空格隔开，表示一辆货车需要从 x 城市运输货物到 y 城市，注意： x 不等于 y 。

输出格式：
共有 qq 行，每行一个整数，表示对于每一辆货车，它的最大载重是多少。如果货车不能到达目的地，输出 -1−1

---

并查集+最大生成树+LCA倍增

ps：打完以后发现可以树链剖分，但是懒得重新写了

洛谷一直卡在95分，后来发现是忘了考虑图的联通性问题（只dfs（）了一次）

---

AC代码如下

#include<cstdio>
#include<cstring>
#include<algorithm>
using namespace std;
struct edge
{
    int x,y,z;
 }e[500005];
struct edge2
{
    int length,to,next;
}E[20005];
bool v[10005];
int n,tot,a,b,Q,m,f[10005][30],len[10005][30],fa[10005],d[10005],head[10005],lg[100005];
void add(int x,int y,int z)
{
    E[++tot].next=head[x];E[tot].length=z;E[tot].to=y;head[x]=tot;
}
bool cmp(edge a,edge b){return a.z>b.z;}
int find(int x)
{
    if(fa[x]==x)return x;
    else return fa[x]=find(fa[x]);
}
void dfs(int father,int x,int value)
{
    d[x]=d[father]+1;
    f[x][0]=father;len[x][0]=value;
    for(int i=1;(1<<i)<d[x];i++)
    {
        f[x][i]=f[f[x][i-1]][i-1];
        len[x][i]=min(len[x][i-1],len[f[x][i-1]][i-1]);
    }
    for(int i=head[x];i;i=E[i].next)
      if(E[i].to!=father)
        dfs(x,E[i].to,E[i].length);
}
int ask(int x,int y)
{
    if(d[x]<d[y]){int tmp=x;x=y;y=tmp;}
    int minn=1e9;
    while(d[x]>d[y])
    {
        minn=min(minn,len[x][lg[d[x]-d[y]]-1]);
        x=f[x][lg[d[x]-d[y]]-1];
    }
    if(x==y)return minn;
    for(int  k=lg[d[x]];k>=0;k--)
      if(f[x][k]!=f[y][k])
      {
        minn=min(minn,len[x][k]);x=f[x][k];
        minn=min(minn,len[y][k]);y=f[y][k];
      }
    minn=min(minn,len[x][0]);
    minn=min(minn,len[y][0]);
    return minn;
}
int main()
{
    memset(len,0x7f,sizeof(len));
    scanf("%d%d",&n,&m);
    for(int i=1;i<=m;i++)scanf("%d%d%d",&e[i].x,&e[i].y,&e[i].z);
    sort(e+1,e+m+1,cmp);
    for(int i=1;i<=n;i++)fa[i]=i,lg[i]=lg[i-1]+(1<<lg[i-1]==i);
    for(int i=1;i<=m;i++)
    {
        int p=find(e[i].x),q=find(e[i].y);
        if(p==q)continue;
        fa[p]=q;
        add(e[i].x,e[i].y,e[i].z);add(e[i].y,e[i].x,e[i].z);
    }
    for(int i=1;i<=n;i++)
      if(!v[find(i)])
        dfs(0,find(i),0),v[find(i)]=1;
    scanf("%d",&Q);
    while(Q--)
    {
        scanf("%d%d",&a,&b);
        if(find(a)!=find(b))printf("-1\n");
        else printf("%d\n",ask(a,b));
    }
    return 0;

}