codevs3731 luogu2296 寻找道路

22 篇文章 0 订阅
3 篇文章 0 订阅
博客探讨了在有向图中寻找满足特定条件的最短路径问题,即路径上的所有点都能直接或间接到达终点。文章通过实例解释了当起点与终点不连通时,路径不存在的情况。博主尝试使用BFS和Dijkstra算法解决,通过建立反向图,从终点开始找到无法到达的点,并在起点进行Dijkstra搜索,避免选择特定标记的点。在调试过程中,博主强调了初始化、判断条件和数据结构的重要性,并提到可能需要堆优化来提高效率。
摘要由CSDN通过智能技术生成

题目描述

在有向图G 中,每条边的长度均为1 ,现给定起点和终点,请你在图中找一条从起点到终点的路径,该路径满足以下条件:

1 .路径上的所有点的出边所指向的点都直接或间接与终点连通。

2 .在满足条件1 的情况下使路径最短。

注意:图G 中可能存在重边和自环,题目保证终点没有出边。

请你输出符合条件的路径的长度。

输入输出格式

输入格式:

输入文件名为road .in。

第一行有两个用一个空格隔开的整数n 和m ,表示图有n 个点和m 条边。

接下来的m 行每行2 个整数x 、y ,之间用一个空格隔开,表示有一条边从点x 指向点y 。

最后一行有两个用一个空格隔开的整数s 、t ,表示起点为s ,终点为t 。

输出格式:

输出文件名为road .out 。

输出只有一行,包含一个整数,表示满足题目᧿述的最短路径的长度。如果这样的路径不存在,输出- 1 。

输入输出样例

输入样例#1:
3 2  
1 2  
2 1  
1 3  
输出样例#1:
-1
输入样例#2:
6 6  
1 2  
1 3  
2 6  
2 5  
4 5  
3 4  
1 5  
输出样例#2:
3

说明

解释1:

如上图所示,箭头表示有向道路,圆点表示城市。起点1 与终点3 不连通,所以满足题

目叙述的路径不存在,故输出- 1 。

解释2:

如上图所示,满足条件的路径为1 - >3- >4- >5。注意点2 不能在答案路径中,因为点2连了一条边到点6 ,而点6 不与终点5 连通。

对于30%的数据,0<n≤10,0<m≤20;

对于60%的数据,0<n≤100,0<m≤2000;

对于100%的数据,0<n≤10,000,0<m≤200,000,0<x,y,s,t≤n,x≠t。



图论第二题

课件里面给出的标算是bfs求最短路

【不知道怎么写bfs的心痛

于是开始自己折腾

考虑再建立一个反向的图

从终点做一遍dijkstra

找出终点走不到的点

将出边指向这些点的点标记

从起点做dijkstra

不选被标记的点

写起来很流畅

可调bug用了很长时间

每次调bug都会有很深的体会

1、赋初值真的要细致

2、dijkstra一开始千万不要写vis[sx]=1;【wa了很多次

3、memset函数将dis数组初始化为无穷大0x7f 事实0x7f直接输出为127 在判断点与终点是否相通时 不能把dis[i]与0x7f直接比较

4、出边指向某一特殊点的点不止一个 所以pre数组不能简单的定义为一维数组 可考虑结构体一维数组里面分装vector

嗯…codevs上ac了 luogu上t了一个点

可能是dijkstra没写堆优化吧【事实上已经不会写堆了…

明天继续加油吧


#include<bits/stdc++.h>
#define maxn 10001
#define maxm 200001
using namespace std;
template <typename T> void read(T &x){
	x=0;int f=1;char ch=getchar();
	for(;!isdigit(ch);ch=getchar())if(ch=='-')f=-1;
	for(;isdigit(ch);ch=getchar())x=(x<<1)+(x<<3)+ch-'0';
	x*=f;
}

int to1[maxm],nxt1[maxm],head1[maxn],total1;//正搜 
int to2[maxm],nxt2[maxm],head2[maxn],total2;//反搜 
int dis[maxn];
//int pre[maxn];
bool vis[maxn];
bool flag[maxn];
int n,m;
int sx,ex;

struct node{
	vector<int> vec;
}pre[maxn];

void add1(int u,int v){
	++total1;
	to1[total1]=v;
	nxt1[total1]=head1[u];
	head1[u]=total1;
}

void add2(int u,int v){
	++total2;
	to2[total2]=v;
	nxt2[total2]=head2[u];
	head2[u]=total2;
}

void dijkstra1(){
	memset(dis,0x7f,sizeof(dis));
	memset(vis,0,sizeof(vis));
	memset(flag,0,sizeof(flag));
	dis[ex]=0;
	for(int ti=1;ti<=n;++ti){
		int jd=-1;
		for(int i=1;i<=n;++i)
			if(!vis[i]&&(jd==-1||dis[i]<dis[jd])) jd=i;
		vis[jd]=1;
		for(int e=head2[jd];e;e=nxt2[e])
			dis[to2[e]]=min(dis[to2[e]],dis[jd]+1);
	}
	//for(int i=1;i<=n;++i) cout<<dis[i]<<" "; cout<<endl;
	//cout<<0x7f<<endl;
	for(int i=1;i<=n;++i)
		if(dis[i]>maxm){
			//cout<<i<<" "<<pre[i]<<endl;
			//flag[pre[i]]=1;
			for(int j=0;j<pre[i].vec.size();++j) flag[pre[i].vec[j]]=1;
		}
}

void dijkstra2(){
	memset(dis,0x7f,sizeof(dis));
	memset(vis,0,sizeof(vis));
	//vis[sx]=1;
	dis[sx]=0;
	for(int ti=1;ti<=n;++ti){
		int jd=-1;
		for(int i=1;i<=n;++i)
			if(!vis[i]&&(jd==-1||dis[i]<dis[jd])&&!flag[i]) jd=i;
		vis[jd]=1;
		for(int e=head1[jd];e;e=nxt1[e]){
			if(!flag[to1[e]]) 
			dis[to1[e]]=min(dis[to1[e]],dis[jd]+1);
		}
	}
}

int main(){
	read(n),read(m);
	memset(pre,0,sizeof(pre));
	for(int i=1;i<=m;++i){
		int u,v;
		read(u),read(v);
		pre[v].vec.push_back(u);
		add1(u,v);
		add2(v,u);
	}
	//for(int i=1;i<=n;++i) cout<<pre[i]<<" "; cout<<endl;
	read(sx),read(ex);
	dijkstra1();
	dijkstra2();
	//for(int i=1;i<=n;++i) cout<<flag[i]<<" "; cout<<endl;
	if(dis[ex]>maxm*100)  cout<<-1<<endl;
	else cout<<dis[ex]<<endl; 
 	return 0;
}



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

请填写红包祝福语或标题

红包个数最小为10个

红包金额最低5元

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

抵扣说明:

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

余额充值