#二进制分组+dij# [校测 普通快乐]

本文介绍了一种针对特殊节点的最短路径问题解决方案。通过Dijkstra算法变体,结合优先队列实现高效寻径,旨在寻找从特殊节点出发到达任意非特殊节点的次短路径与最大权值边最小值。

摘要生成于 C知道 ,由 DeepSeek-R1 满血版支持, 前往体验 >

Title

在这里插入图片描述

给旅行者很像,但是这道题目是单向边。


Solution

非正解: 把所有特殊点压进队列里,然后跑dij,记录次小值和最大值,可以枚举任一个非特殊点求解;
正解: 在这里插入图片描述


Code(非正解)

#include<cstdio> 
#include<algorithm>
#include<cstring>
#include<iostream>
#include<queue>
#define mp(x,y) make_pair((x),(y))
#define rep(i,x,y) for(int i=x;i<=y;i++)
using namespace std; 
const int N=2e5+20; 
int read(){
	int p=0,f=1; char c=getchar(); 
	while (!isdigit(c)) {if (c=='-') f=-1; c=getchar();}
	while (isdigit(c)) p=(p<<3)+(p<<1)+c-48,c=getchar(); 
	return p*f; c 
}
int n,m,K,b[N],d[N],tot,head[N],ans=2147483647,D[N]; 
bool vis[N],B[N]; 
int X[N],Y[N],Z[N]; 
priority_queue<pair<int,int> >q; 
struct node{
	int y,z,nt; 
}a[N*2];
void add(int x,int y,int z){
	a[++tot]=(node){y,z,head[x]}; head[x]=tot; 
}
void dij(){
	while(q.size()){
		int x=q.top().second; q.pop(); 
		if (!vis[x]) vis[x]=1; else continue; 
		for(int i=head[x];i;i=a[i].nt){
			int y=a[i].y; 
			if (d[y]>=d[x]+a[i].z){
				D[y]=d[y],d[y]=d[x]+a[i].z; 
				q.push(mp(-d[y],y)); 
			}
		}
	}
}
int main(){
	memset(d,0x3f,sizeof(d)); 
	memset(D,0x3f,sizeof(D)); 
	memset(vis,0,sizeof(vis)); 
	n=read(),m=read(),K=read(); 
	rep(i,1,m){
		int x=read(),y=read(),z=read(); 
		add(x,y,z),add(y,x,z); 
		X[i]=x,Y[i]=y,Z[i]=z; 
	}
	rep(i,1,K) b[i]=read(),q.push(mp(0,b[i])),d[b[i]]=0,D[b[i]]=0,B[b[i]]=1;  
	dij();
	rep(i,1,m) if (B[X[i]]&&B[Y[i]]) ans=min(Z[i],ans); 
	rep(i,1,n) if (!B[i]) ans=min(ans,D[i]+d[i]); 
	printf("%d",ans); 
	return 0; 
} 

Code(std)

#include<cstdio>
#include<cstdlib>
#include<cstring>
#include<iostream>
#include<queue>
using namespace std;

int k,n;
int m;
struct edge{
	int y,next,c;
}s[400010];
struct node{
	int x,c;
	bool operator<(const node a)const{
		return c>a.c;
	}
};
int first[100010],len=0;
int g[100010],dis[100010];
priority_queue<node> f;
bool tf[100010],we[100010];
int ans=1e9;

void read(int&x){
	char ch=getchar();x=0;
	while(ch<'0' || ch>'9') ch=getchar();
	while(ch>='0' && ch<='9') x=x*10+ch-'0',ch=getchar();
}

void ins(int x,int y,int c){
	s[++len]=(edge){y,first[x],c};first[x]=len;
}

void solve(int x){
	for(int i=1;i<=n;i++) dis[i]=1e9,tf[i]=false;
	for(int i=1;i<=k;i++) if((g[i]>>x)&1) f.push((node){g[i],0}),dis[g[i]]=0;
	while(!f.empty()){
		node x=f.top();f.pop();
		if(tf[x.x]) continue;
		tf[x.x]=true;
		for(int i=first[x.x];i!=0;i=s[i].next)if(dis[x.x]+s[i].c<dis[s[i].y]) {
			dis[s[i].y]=dis[x.x]+s[i].c;
			f.push((node){s[i].y,dis[s[i].y]});
		}
	}
	for(int i=1;i<=k;i++) if(!((g[i]>>x)&1)) ans=min(ans,dis[g[i]]);
}

int main(){
	read(n);read(m);read(k); 
	int x,y,c;
	for(int i=1;i<=m;i++){
		read(x);read(y);read(c);
		ins(x,y,c);ins(y,x,c);
	}
	for(int i=1;i<=k;i++) read(g[i]);
	for(int i=0;i<=16;i++) solve(i);
	printf("%d\n",ans);
}
评论
添加红包

请填写红包祝福语或标题

红包个数最小为10个

红包金额最低5元

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

抵扣说明:

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

余额充值