niop 2014寻找道路

最新推荐文章于 2021-10-19 06:26:35 发布

转载最新推荐文章于 2021-10-19 06:26:35 发布 · 105 阅读

0 ·

CC 4.0 BY-SA版权

原文链接：http://www.cnblogs.com/yanlifneg/p/5418866.html

本文介绍了一种针对 SPFA 最短路径算法的优化方案，确保路径上所有节点的出边指向的节点都能直接或间接与终点相连。通过逆向 BFS 标记可达节点并筛选有效节点的方式进行优化。

/*
乍一看就是个最短路 SFPA 
但是要保证路径上的所有点的出边所指向的点都直接或间接与终点连通。
这一点就蛋疼了0.0
开始想的是正着跑一边 每一个点的所有边都能符合条件 那这个点就符合条件0.0  可惜WA +RE 了
后来xyd大神说 可以先到这从终点跑一边 能到的标记 然后在枚举每个点的所有边指向的点 全被标记那么这个点就ok 
最后SPFA 
 注意 正反建边 数组要大！！！！！！！！
 还有 从终点跑一边 能到的标记时Dfs会TLE 所以0.0 BFs了 
*/
#include<iostream>
#include<cstdio>
#include<cstring>
#include<queue>
#define maxn 10010
#define maxx 400010
using namespace std;
queue<int>q;
int n,m,head1[maxx],head2[maxx],num,s,t,f[maxn],ans[maxn],cn1[maxn],cn2[maxn];
struct node
{
    int v,dis,pre;
}e1[maxx],e2[maxx];
/*int Dfs(int p)
{
    f[p]=1;
    for(int i=head2[p];i;i=e2[i].pre)
      {
          cn1[e2[i].v]=1;
          if(f[e2[i].v]==0)
            Dfs(e2[i].v);
      }
}*/
void Bfs(int p)
{
    q.push(p);
    f[p]=1;
    int i;
    while(!q.empty())
      {
          int k=q.front();
          q.pop();
        for(i=head2[k];i;i=e2[i].pre)
        {
            cn1[e2[i].v]=1;
            if(f[e2[i].v]==0)
              {
                q.push(e2[i].v);
                f[e2[i].v]=1;
            }
        }
      }
}
void Biu()
{
    int i,j;
    for(i=1;i<=n;i++)
      {
          int g=0;
          for(j=head1[i];j;j=e1[j].pre)
          if(cn1[e1[j].v]==0)
            g=1;
        if(g==0)
          cn2[i]=1;
      }
}
int main()
{
    scanf("%d%d",&n,&m);
    int i,x,y;
    memset(ans,127,sizeof(ans));
    for(i=1;i<=m;i++)
      {
          scanf("%d%d",&x,&y);
          num++;
          e1[num].v=y;
          e1[num].dis=1;
          e1[num].pre=head1[x];
          head1[x]=num;
          num++;
          e2[num].v=x;
          e2[num].dis=1;
          e2[num].pre=head2[y];
          head2[y]=num;
      }
    scanf("%d%d",&s,&t);
    Bfs(t);
    cn1[t]=1;
    Biu();
    memset(f,0,sizeof(f));
    q.push(s);
    ans[s]=0;
    f[s]=1;
    while(!q.empty())
      {
          int k=q.front();
          q.pop();
        f[k]=0;
          for(i=head1[k];i;i=e1[i].pre)
            if(ans[e1[i].v]>ans[k]+e1[i].dis&&cn2[k]==1)
              {
                ans[e1[i].v]=ans[k]+e1[i].dis;
                if(f[e1[i].v]==0)
                  {
                    q.push(e1[i].v);
                    f[e1[i].v]=1;
                }
            }
      }
    if(ans[t]<0x7fffffff)printf("%d",ans[t]);
    else printf("-1");
    return 0;
}

转载于:https://www.cnblogs.com/yanlifneg/p/5418866.html