codevs3731 luogu2296 寻找道路

题目描述

在有向图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;
}