Luogu 1967(货车运输)(最大生成森林+LCA)

传送门

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

题解:由于要装尽量多的货物,所以先用kruskal跑一个最大生成森林,然后再因为存在短板效应,所以在这些生成的树上跑倍增(可参考倍增求LCA),求出每条询问的路径上的最小边权。

 

注意:生成森林的father数组最后一定要再find一遍把应该连的都连好。为什么单纯做并查集时不用再练一波呢?因为并查集的query都是直接调用的find函数,没连好的在递归时自己就连上了。

当然如果在询问时调用find的话也OK,就不用在提前find这一遍。

#include<cstdio>
#include<cstring>
#include<iostream>
#include<algorithm>
using namespace std;
const int N=1e4+4,M=5e4+4,INF=0x3f3f3f3f;
int n,m,fa[N],etot=0;
int head[N],dep[N];
int mn[15][N],f[15][N];
struct A {
	int u,v,w;
	friend bool operator <(const A &p,const A &q) {
		return p.w>q.w;
	}
}a[M];
struct EDGE {
	int v,nxt,w;
}e[M<<1];
inline int read() {
	int x=0;char c=getchar();
	while (c<'0'||c>'9') c=getchar();
	while (c>='0'&&c<='9') x=x*10+c-'0',c=getchar();
	return x;
}
inline void adde(int u,int v,int w) {
	e[++etot].nxt=head[u],e[etot].v=v,e[etot].w=w,head[u]=etot;
	e[++etot].nxt=head[v],e[etot].v=u,e[etot].w=w,head[v]=etot;
}
int find(int x) {
	return x^fa[x]?fa[x]=find(fa[x]):x;
}
void dfs(int p,int father) {
	for (int j=1;(1<<j)<=dep[p];++j) {
		f[j][p]=f[j-1][f[j-1][p]];
		mn[j][p]=min(mn[j-1][p],mn[j-1][f[j-1][p]]);
	}
	for (int i=head[p];~i;i=e[i].nxt) {
		int v=e[i].v;
		if (v^father) {
			dep[v]=dep[p]+1;
			f[0][v]=p;
			mn[0][v]=e[i].w;
			dfs(v,p);
		}
	}
}
inline int query(int x,int y) {
	if (fa[x]^fa[y]) return -1;
	int ret=INF;
	if (dep[x]<dep[y]) x^=y^=x^=y;
	int t=dep[x]-dep[y];
	for (int i=14;~i;--i)
		if (t&(1<<i))
			ret=min(ret,mn[i][x]),x=f[i][x];
	if (x==y) return ret;
	for (int i=14;~i;--i)
		if (f[i][x]^f[i][y])
			ret=min(ret,min(mn[i][x],mn[i][y])),x=f[i][x],y=f[i][y];
	return ret=min(ret,min(mn[0][x],mn[0][y]));
}
int main() {
//	freopen("1.txt","r",stdin);
	memset(head,-1,sizeof(head));
	n=read(),m=read();
	for (register int i=1;i<=m;++i)
		a[i].u=read(),a[i].v=read(),a[i].w=read();
	sort(a+1,a+m+1);
	int blo=n;
	for (register int i=1;i<=n;++i) fa[i]=i;
	for (register int i=1;i<=m&&blo>1;++i) {
		int u=find(a[i].u),v=find(a[i].v);
		if (u^v) adde(a[i].u,a[i].v,a[i].w),fa[v]=u,--blo;
	}
	memset(mn,INF,sizeof(mn));
	for (register int i=1;i<=n;++i) {
		find(i);//!!!Must merge again!
		if (!dep[i]) dep[i]=1,dfs(i,0);
	}
	int q=read();
	for (register int i=0;i<q;++i) {
		int u=read(),v=read();
		printf("%d\n",query(u,v));
	}
	return 0;
}

 

  • 0
    点赞
  • 0
    收藏
    觉得还不错? 一键收藏
  • 0
    评论
评论
添加红包

请填写红包祝福语或标题

红包个数最小为10个

红包金额最低5元

当前余额3.43前往充值 >
需支付:10.00
成就一亿技术人!
领取后你会自动成为博主和红包主的粉丝 规则
hope_wisdom
发出的红包
实付
使用余额支付
点击重新获取
扫码支付
钱包余额 0

抵扣说明:

1.余额是钱包充值的虚拟货币,按照1:1的比例进行支付金额的抵扣。
2.余额无法直接购买下载,可以购买VIP、付费专栏及课程。

余额充值