原题的旅行

题解

先从全都是文明城市情况考虑:显然就是 x x x y y y 的路径上最大的边的最小值,直接重构树即可。

考虑非文明城市的加入:

因为从 x x x y y y 的路径一定是经过若干个文明城市,所以这时路径的权值就定义为:两个文明城市直接的最短路。

然后同上,问题是这样的边是 n 2 n^2 n2 的,考虑从性质入手优化:

对于 M S T MST MST 上的一条边,对应原图上两个文明城市 i , j i,j i,j 的一条最短路的权值和,因为边 i → j i \to j ij M S T MST MST 上,所以有更优的选择使得不需要连 i → j i \to j ij 而让 i , j i,j i,j 连通,所以记录每个点到所有文明城市的最短距离以及是哪个城市, i , j i,j i,j 就一定会被路径上的一条边枚举到。

注意,该题卡spfa

题解

#include<bits/stdc++.h>
#define LL long long
using namespace std;
const int N=2e5+10;
char c[N];
LL d[N],val[N<<1];
int n,v[N],m,be[N],cnt,fa[N<<1],f[N<<1][22],dep[N<<1],now,Q;
int tot,ver[N<<1],fst[N],nxt[N<<1],w[N<<1];
inline void add(int x,int y,int z){ver[++tot]=y;nxt[tot]=fst[x];fst[x]=tot;w[tot]=z;}
vector<int> to[N<<1];
priority_queue<pair<int,int> >q;
struct Node{
	int x,y;
	LL w;
}e[N<<2];
bool cmp(const Node x,const Node y){return x.w<y.w;}
inline int get(int x){return (fa[x]==x)?x:fa[x]=get(fa[x]);}
void dfs1(int x){
	for(auto y:to[x]){
		dep[y]=dep[x]+1;
		f[y][0]=x;
		dfs1(y);
	}
}
inline int LCA(int x,int y){
	if(dep[x]<dep[y])swap(x,y);
	int cha=dep[x]-dep[y];
	for(int i=20;i>=0;--i)if(cha&(1<<i))x=f[x][i];
	if(x==y)return x;
	for(int i=20;i>=0;--i)if(f[x][i]!=f[y][i])x=f[x][i],y=f[y][i];
	return f[x][0];
}
int main(){
	scanf("%d%d",&n,&m);
	scanf("%s",c+1);
	for(int i=1;i<=m;++i){
		int x,y,z;
		scanf("%d%d%d",&x,&y,&z);
		add(x,y,z);add(y,x,z);
	}
	for(int i=1;i<=n;++i)d[i]=1e17;
	for(int i=1;i<=n;++i)if(c[i]=='1')q.push(make_pair(0,i)),v[i]=1,d[i]=0,be[i]=i;
	while(q.size()){
		int x=q.top().second,y=-q.top().first;q.pop();
		if(d[x]!=y)continue;
	//	cout<<"FAQ "<<x<<" "<<d[x]<<" "<<be[x]<<endl;
		v[x]=1;
		for(int i=fst[x];i;i=nxt[i]){
			int y=ver[i];
			if(v[y])continue;
			if(d[y]>d[x]+w[i]){
				d[y]=d[x]+w[i];
				be[y]=be[x];
				q.push(make_pair(-d[y],y));
			}
		}
	}
//	for(int i=1;i<=n;++i)cout<<"FAQ "<<i<<" "<<d[i]<<" "<<be[i]<<endl;
	for(int i=1;i<=n;++i){
		for(int j=fst[i];j;j=nxt[j]){
			int y=ver[j];
			if(be[i]!=be[y])e[++cnt]=(Node){be[i],be[y],d[i]+d[y]+w[j]};
		}
	}
	sort(e+1,e+cnt+1,cmp);
	now=n;
	for(int i=1;i<=n;++i)fa[i]=i;
	for(int i=1;i<=cnt;++i){
		int x=get(e[i].x),y=get(e[i].y);
		if(x==y)continue;
		now++;
		fa[now]=now;
		fa[x]=fa[y]=now;
		to[now].push_back(x);
		to[now].push_back(y);
		val[now]=e[i].w;
	}
	dfs1(now);
	for(int j=1;j<=20;++j)for(int i=1;i<=now;++i)f[i][j]=f[f[i][j-1]][j-1];
	scanf("%d",&Q);
	while(Q--){
		int x,y;
		scanf("%d%d",&x,&y);
		printf("%lld\n",val[LCA(x,y)]);
	}
	return 0;
}
  • 0
    点赞
  • 0
    收藏
    觉得还不错? 一键收藏
  • 0
    评论
评论
添加红包

请填写红包祝福语或标题

红包个数最小为10个

红包金额最低5元

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

抵扣说明:

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

余额充值