负环

(tips:这里有一道题我交了不下30次,刚刚ac。。。)
首先,我们来了解一下负环的定义。
负环,又叫负权回路,负权环,指的是一个图中存在一个环,里面包含的边的边权总和<0。在存在负环的图中,是求不出最短路径的,因为只要在这个环上不停的兜圈子,最短路径就会无限小。
看到洛谷里各种玄学算法(感觉洛谷里的题解都很迷 ),本蒟蒻还是给大家了解一下bfs版spfa判负环的这种方法(dfs?不存在的 )。
判断负权回路的方法有两种
1、如果一个点的入队的次数大于点的总数N,则存在负权回路。
2、如果每个顶点的最短路径超过N,则存在负权回路。(雾
所以,我先给大家讲解第一种方法。
首先,一个很基本的存图操作。

struct node
{
 int next,to,w;
}edge[maxn*2];
void addedge(int x,int y,int z)
{
    edge[++tot].to=y;
    edge[tot].next=head[x];
    edge[tot].w=z;
    head[x]=tot;
}

建完图后,跑一遍spfa
我们用circle数组保存每个点的入队次数。本人英语垃圾,环,circle…

bool spfa()
{
 queue<int>q;
 q.push(1);
 circle[1]=1;//1号点入队一次。
    while(!q.empty())
    {
        int u=q.front();
        q.pop();
  flag[u]=false;
        for(int i=head[u];i;i=edge[i].next)
        {//遍历与队首节点相连的点。
            int v=edge[i].to;
            int w=edge[i].w;
            if(dis[v]>dis[u]+w)
            {
                dis[v]=dis[u]+w;//松弛操作
                if(!flag[v])
                {
                 circle[v]=circle[u]+1;//关键,这个点的入队次数是他上一个相连的点的次数+1。
                 if(circle[v]>n)return true;//第一种方法的定义,存在负环,退出函数。
                  q.push(v);
                 flag[v]=true;
                }
            }
        }
    }
    return false;//不存在负环
}


推荐题目
1、【模板】负环
(tips:输出的是YE5,不是S,N0,不是O,不要被题解上那些玄学算法迷惑。。。)。
2、P2850 [USACO06DEC]虫洞Wormholes
3、洛谷搜标签:负权环。

Code:
1、【模板】负环

#include<bits/stdc++.h>
using namespace std;
#define maxn 6001
#define INF INT_MAX
int head[maxn];
int n,m,p,tot=0;
int dis[2*maxn];
bool flag[2*maxn];
int ans=0;
int circle[2*maxn];
queue<int>q;
struct node
{
 int next,to,w;
}edge[maxn*2];
void addedge(int x,int y,int z)
{
    edge[++tot].to=y;
    edge[tot].next=head[x];
    edge[tot].w=z;
    head[x]=tot;
}
bool spfa()
{
 queue<int>q;
 q.push(1);
 circle[1]=1;
    while(!q.empty())
    {
        int u=q.front();
        q.pop();
  flag[u]=false;
        for(int i=head[u];i;i=edge[i].next)
        {
            int v=edge[i].to;
            int w=edge[i].w;
            if(dis[v]>dis[u]+w)
            {
                dis[v]=dis[u]+w;
                if(!flag[v])
                {
                 circle[v]=circle[u]+1;
                 if(circle[u]>n)return true;
                  q.push(v);
                 flag[v]=true;
                }
            }
        }
    }
    return false;
}
int main()
{
    int T;
    scanf("%d",&T);
    while(T--)
    {
        scanf("%d%d",&n,&m);
        tot=0;
        for(int i=1;i<=n;i++)
        {
         flag[i]=false;
         dis[i]=0;
         circle[i]=0;
         head[i]=0;
  }
  dis[1]=0;
        for(int i=1,u,v,w;i<=m;i++) 
        {
            scanf("%d%d%d",&u,&v,&w);
            addedge(u,v,w);
            if(w>=0)addedge(v,u,w);
        }       
        if(spfa())puts("YE5");
        else puts("N0");
    }
    return 0;
}

2、虫洞

#include<bits/stdc++.h>
using namespace std;
#define maxn 500001
#define INF INT_MAX
int pre[maxn*2],now[maxn*2],son[maxn*2],len[maxn];
int n,m,p,tot=0;
int dis[maxn];
bool flag[maxn];
int ans=0;
int circle[maxn];
queue<int>q;
inline void edge(int x,int y,int z)
{
    pre[++tot]=now[x];
    now[x]=tot;
    son[tot]=y;
    len[tot]=z;
}
void spfa(int p)
{
 memset(dis,0x3f,sizeof(dis));
 q.push(p);
 dis[p]=0;
 flag[p]=true;
 while(!q.empty())
 {
  int u=q.front();
  q.pop();
  flag[u]=false;
  for(int i=now[u];i;i=pre[i])
  {
   int v=son[i];
   if(dis[v]>dis[u]+len[i])
   {
    dis[v]=dis[u]+len[i];
    circle[v]=circle[u]+1;
    if(circle[v]>n){puts("YES");return ;}
    if(!flag[v])
    {
     flag[v]=false;
     q.push(v);
    }
   }
  }
 }
 puts("NO");
} 
int main()
{
 int T,holes;
 scanf("%d",&T);
    for(register int i=1;i<=T;i++)
    {
        scanf("%d%d%d",&n,&m,&holes);
        for(register int i=1,u,v,w;i<=m;i++) 
        {
            scanf("%d%d%d",&u,&v,&w);
            edge(u,v,w);
            edge(v,u,w);
        }
        for(register int i=1,u,v,w;i<=holes;i++)
        {
            scanf("%d%d%d",&u,&v,&w);
            edge(u,v,-w);
        }
        spfa(1);
        memset(len,0,sizeof(len));
        memset(pre,0,sizeof(pre));
        memset(son,0,sizeof(son));
        memset(now,0,sizeof(now));
        memset(circle,0,sizeof(circle));
    }
 return 0;
}

不要在意两次存图方法不一样 )。

本篇题解到此结束,谢谢大家的查看。
提供hack数据请在评论区评论,谢谢。

### SPFA算法检测的实现及原理 SPFA(Shortest Path Faster Algorithm)是一种用于解决单源最短路径问题的有效算法,尤其适用于含有权边但不包含的情况。然而,在实际应用中,如果图中存在,SPFA可能陷入无限循,因此需要一种机制来检测这种特殊情况。 #### 的存在条件 指的是在一个图中存在一条回路,该回路上所有边的权重之和为数。由于每次经过这条回路都会使路径长度变得更小,理论上可以不断绕圈从而得到无穷小的距离值[^1]。 #### SPFA算法的工作方式 SPFA的核心思想是利用队列进行广度优先搜索(BFS),并通过松弛操作更新节点到起点之间的距离。具体来说,对于每一个从队列弹出的顶点u,遍历它所有的邻接点v,并尝试用当前已知的最佳路径加上(u,v)这条边的成本去改进到达v点的距离。一旦发现新的更优解,则将v重新放入队列等待进一步探索[^3]。 #### 检测的关键技术 尽管标准版SPFA不具备内置功能自动识别,但可以通过增加额外计数器变量`cnt[]`记录每个结点被访问(即入队)的最大次数来间接完成此任务: - **初始化阶段**: 设置数组 `dis[]` 表示各顶点至源点间的最小估计成本;设初值均为正无穷大(+∞),除了起始位置设置成零外(`dis[s]=0`)。另外定义辅助数组 `inQueue[]`,用来标记哪些顶点正处于队列之中以防重复插入相同元素浪费时间资源。 - **核心逻辑修改**: - 当某顶点再次进入队列前检查它的累计入队数目是否已经超过整个网络中的总定点数量V; 如果满足上述条件则立即终止运算并报告发现了不可接受状况—存在至少一个权闭合链表结构[^2]. 以下是具体的伪代码描述如何扩展基础版本以支持探测: ```python from collections import deque def spfa_with_negative_cycle_detection(graph, start_node, num_nodes): dis = [float('inf')] * (num_nodes + 1) cnt = [0] * (num_nodes + 1) in_queue = [False] * (num_nodes + 1) queue = deque() dis[start_node] = 0 queue.append(start_node) in_queue[start_node] = True while queue: u = queue.popleft() in_queue[u] = False for v, weight in graph[u]: if dis[v] > dis[u] + weight: dis[v] = dis[u] + weight if not in_queue[v]: cnt[v] += 1 # If a node is relaxed more than the number of nodes times, # then there must be a negative cycle. if cnt[v] >= num_nodes: return "Negative Cycle Detected" queue.append(v) in_queue[v] = True return "No Negative Cycles Found" ``` 在这个增强型函数里,我们引入了一个名为`cnt[]`的新列表跟踪各个节点经历过的放松迭代轮次。每当某个特定节点准备第二次加入工作序列之前都要先核查其对应统计数值是不是已经达到了预设界限——也就是等于总的节点个数N。如果是的话就意味着系统内部必然隐藏着某种形式上的面反馈闭现象发生,此时应该立刻停止后续计算动作并向外界发出警告信号表明遇到了非法情形[^4]。 ### 结论 综上所述,虽然原始形态下的SPFA并不自带针对形拓扑结构的有效甄别手段,不过借助简单的附加措施比如设立专门监控指标就可以轻松弥补这项缺陷进而构建更加健壮可靠的解决方案出来供人们日常开发实践当中选用。
评论 6
添加红包

请填写红包祝福语或标题

红包个数最小为10个

红包金额最低5元

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

抵扣说明:

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

余额充值