【noip2014】寻找道路

题目描述

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

15min想出来思路，但是一直调了两节课。正着不好想的话，可以反过来，建一个反向的图，从终点出发标记，这些点是可以走的。
还是有很多小错误：
1）首先是spfa，dis初始值应该是0x7ffffffff，起点的dis值是0，然后 if (b[edge[i].to]&&dis[edge[i].to]>dis[now]+1)
又分不清大于号小于号！！
2）在标记的时候有一个坑，还是审题的问题“1 ．路径上的所有点的出边所指向的点都直接或间接与终点连通。”是指路径上的点指向的点只要有一条能通到终点就可以。所以单独开一个can数组来标记，然后在复制到b数组上
3）！！！在函数内定义数组一定要memset清零！！或者不要在函数内定义

#include<cstdio>#include<cstring>#include<algorithm>using namespace std;const int MAXN=10001;const int MAXM=400001;int n,m,x,y,s,t;bool b[MAXN];struct Edge{    int next,to;}edge[MAXM];struct E{    int next,to;}e[MAXM];int read(){    int now=0,f=1; char ch=getchar();    while (ch<'0'||ch>'9') {if (ch=='-') f=-1; ch=getchar();}    while (ch>='0'&&ch<='9')    {        now=(now<<1)+(now<<3)+ch-'0';        ch=getchar();    }    return now*f;}int head[MAXM],num_edge;void add_edge(int from,int to){    edge[++num_edge].next=head[from];    edge[num_edge].to=to;    head[from]=num_edge;}int h[MAXM],num_e;void add_e(int from,int to){    e[++num_e].next=h[from];    e[num_edge].to=to;    h[from]=num_e;}void dfs(int x){    b[x]=true;    for (int i=h[x]; i!=0; i=e[i].next)        if (!b[e[i].to])        {            dfs(e[i].to);        }}int team[MAXM],dis[MAXM];bool vis[MAXM];void SPFA(int s){//  memset(dis,0x7f/3,sizeof(dis));    for (int i=1; i<=n; i++) if (i!=s) dis[i]=0x7fffffff/3;    int front=0,tail=1;    team[1]=s; vis[s]=true; //dis[s]=0;    while (front<=tail)    {        int now=team[++front];        vis[now]=false;        if (!b[now]) continue;        for (int i=head[now]; i!=0; i=edge[i].next)        {            if (dis[edge[i].to]>dis[now]+1)            {                dis[edge[i].to]=dis[now]+1;                if (!vis[edge[i].to])                {                    vis[edge[i].to]=true;                    team[++tail]=edge[i].to;                }            }        }    }}bool judge(int x){    for (int i=head[x]; i!=0; i=edge[i].next)        if (b[edge[i].to]==0) return false;//只是判断没有相邻的点而已     return true;}int main(){    n=read(); m=read();//  scanf("%d%d",&n,&m);    for (int i=1; i<=m; i++)    {//      scanf("%d%d",&x,&y);        x=read(); y=read();        add_edge(x,y); add_e(y,x);    }    s=read(); t=read();    dfs(t);    bool can[MAXN];//需要重新开一个数组标记      memset(can,0,sizeof(can));/*  for (int i=1; i<=n; i++)        if (!judge(i)) can[i]=1;*///还需要做进一步的判断 (不对，这样多删掉一些点)    for (int i=1; i<=n; i++)        if (!b[i])        for (int j=h[i]; j!=0; j=e[j].next)//i和j又写错了             can[e[j].to]=1;    for (int i=1; i<=n; i++) if (can[i]) b[i]=false;//  for (int i=1; i<=n; i++) if (b[i]) printf("%d ",i); return 0;    if (!b[s]) {printf("-1"); return 0;}    SPFA(s);    if (dis[t]==0x7fffffff/3) printf("-1");    else printf("%d",dis[t]);    return 0;}

%%%女神Cansult：https://cansult.github.io/2017/09/24/翻车笔记_%20noip2014D2T2%20[BFS]/