【Luogu P1967】货车运输
Luogu P1967
题目大意:给定一张图和q个询问,询问x节点和y节点的路径之间最小边权最大可以是多少。
可以发现对于一条边\(E(x,y)\),如果x到y有另一条路径且最小边权大于\(E(x,y)\),那么这条边完全可以不考虑了。
综上所述,我们可以考虑对这张图使用Kruskal构建最大生成树。
那么两个点之间的路径就可以使用倍增LCA来求解了。
注意:原图可能并不连通,构成的可能是一个森林,所以要加入特判。
#include<cstdio>
#include<algorithm>
#include<cstring>
using namespace std;
const int maxm=50005,maxn=10005;
struct data
{
int sta,val,to,next;
bool operator <(const data&x) const
{
return val>x.val;
}
}e[2*maxm],edge[2*maxm];
int cnt,head[maxn],n,q,m,x,y,ans,z,fa[maxn],f[maxn][32],val[maxn][32],d[maxn],lg[maxn];
void add(int u,int v,int w)
{
e[++cnt].sta=u;
e[cnt].to=v;
e[cnt].val=w;
}
int getf(int v)
{
if (fa[v]!=v) return fa[v]=getf(fa[v]);
return fa[v];
}
inline void merge(int x,int y)
{
fa[getf(x)]=getf(y);
}
inline bool check(int x,int y)
{
return fa[getf(x)]==fa[getf(y)];
}
void dfs(int now,int fat)
{
f[now][0]=fat;d[now]=d[fat]+1;
for (int i=1;i<=lg[d[now]];i++)
{
f[now][i]=f[f[now][i-1]][i-1];
val[now][i]=min(val[now][i-1],val[f[now][i-1]][i-1]);
//记录路径上的最小边权,注意理解。
}
for (int i=head[now];i;i=edge[i].next)
{
if (edge[i].to==fat) continue;
val[edge[i].to][0]=edge[i].val;
dfs(edge[i].to,now);
}
}
int lca(int x,int y)
{
if (!check(x,y)) return -1;
if (d[x]<d[y]) swap(x,y);
ans=2147483647;
/*注意这里不能写ans=min(val[x][0],val[y][0]);这两条边不一定是必选的。
因为y可能本身就是lca(x,y)*/
while (d[x]>d[y])
{
ans=min(ans,val[x][lg[d[x]-d[y]]-1]);//注意顺序,更新x之前先更新最小边权。
x=f[x][lg[d[x]-d[y]]-1];
}
if (x==y) return ans;
for (int i=lg[d[x]]-1;i>=0;i--)
if (f[x][i]!=f[y][i])
{
ans=min(ans,min(val[x][i],val[y][i]));
x=f[x][i];y=f[y][i];
//同样注意顺序
}
return ans=min(ans,min(val[x][0],val[y][0]));
}
int main()
{
memset(val,0x3f,sizeof(val));
scanf("%d%d",&n,&m);
for (int i=1;i<=m;i++)
{
scanf("%d%d%d",&x,&y,&z);
add(x,y,z);
add(y,x,z);
}
for (int i=1;i<=n;i++) fa[i]=i;
sort(e+1,e+1+cnt);
int tmpcnt=cnt;
cnt=0;
for (int i=1;i<=tmpcnt;i++)
if (!check(e[i].sta,e[i].to))
{
merge(e[i].sta,e[i].to);
edge[++cnt].to=e[i].to;
edge[cnt].val=e[i].val;
edge[cnt].next=head[e[i].sta];
head[e[i].sta]=cnt;
edge[++cnt].to=e[i].sta;
edge[cnt].val=e[i].val;
edge[cnt].next=head[e[i].to];
head[e[i].to]=cnt;
}
for (int i=1;i<=n;i++) lg[i]=lg[i-1]+(1<<lg[i-1]==i);
for (int i=1;i<=n;i++)
if (!d[i]) dfs(i,0);
scanf("%d",&q);
for (int i=1;i<=q;i++)
{
scanf("%d%d",&x,&y);
printf("%d\n",lca(x,y));
}
return 0;
}