POJ 1086 Currency Exchange(SPFA判断正环)

最新推荐文章于 2021-09-30 21:03:28 发布

KGW_Yuan

最新推荐文章于 2021-09-30 21:03:28 发布

阅读量114

点赞数

分类专栏：笔记文章标签：算法

本文链接：https://blog.csdn.net/m0_50647144/article/details/109721547

版权

笔记专栏收录该内容

5 篇文章 0 订阅

订阅专栏

本文深入探讨了SPFA（Shortest Path Faster Algorithm）算法，它是Bellman-Ford算法的一种优化。SPFA通过保留可能更新最短路径信息的节点来提高效率，避免了Bellman-Ford中无效的边扫描。在判断负环的基础上，通过将距离数组初始化为0并调整更新条件，SPFA也可用于判断正环。文章提供了算法的基本思路和代码实现。

摘要由CSDN通过智能技术生成

spfa的深入理解

https://blog.csdn.net/qq_43093454/article/details/88812117

spfa本身是对Bellman-ford算法的优化。在Bellman-Ford算法中，n轮扫描每一条边看一看各个点有没有希望被松弛。

spfa之所以能用来处理Bellman-ford并进行优化，其原因是在Bellman-ford中，每一轮扫描的边有很多都没有可能更新其它点。spfa指出，当一条边有希望更新它的一个端点的最短路值，当且仅当它的另一个端点的最短路值被更新过。spfa使用一个队列保留这些可能被更新的信息，通常，spfa在队列中保存的更新信息的意义是"这个点被更新过，因而有可能通过与它连接的边更新这些边的其它端点"。直到队列为空，此局面表示"没有一个结点可能更新其它的结点"。

这样保存更新信息的原因是因为，在通常情况下，一个点可以通过与它相连的边更新其它点，每一个结点的更新则是可能由多个通过边连向它的结点决定，即一对多的更新。即因而我们保存的信息是"可能更新其他多个点的点"而不是"可能被其他点更新的点"因为这样做更方便。(想一想，用前者则出队一个结点扫描更新可以将所有由此结点被更新造成的可能被更新的结点的改变情况全部解决，而对于后者，一个结点的更新将所有连向它的结点入队，则当这些结点出队的时候，需要扫描它的所有前驱来决定结果，后者的复杂度显然比前者高很多，相当于n^2和n的区别)。

基本思路
我们知道spfa判断负环的原理，同理在判断正环时我们仅需将dis数组全部初始化为0，更新条件变为if(dis[v]<dis[u]+val)dis[v]=dis[u]+val;即可。

代码实现

#include<iostream>
#include<cstdio>
#include<cstring>
#include<queue>
#include<cstdlib>
#include<algorithm>
using namespace std;
#define maxn 100010
#define ll long long
#define clear(a) memset(a,0,sizeof a)

int n,m,s,tot;
double v;
double val[maxn];
int first[maxn],vis[maxn];
struct edge{
    int nextx,to;
    double rate,com;
}e[maxn<<1];
queue<int>Q;

void add(int u,int v,double r,double c){
    tot++;
    e[tot].nextx=first[u];
    first[u]=tot;
    e[tot].to=v;
    e[tot].rate=r;
    e[tot].com=c;
}

int spfa(int x){
    val[x]=v;
    Q.push(x);
    while(!Q.empty()){
        int u=Q.front();
        Q.pop();
        vis[u]=0;
        if(u==s&&val[u]>v)return 1;
        for(int i=first[u];i;i=e[i].nextx){
            int v=e[i].to;
            
            if(val[v]<(val[u]-e[i].com)*e[i].rate){
                val[v]=(val[u]-e[i].com)*e[i].rate;
                if(!vis[v]){
                    vis[v]=1;
                    Q.push(v);
                }
            }
        }
    }
    return 0;
}

int main(){
    scanf("%d%d%d%lf",&n,&m,&s,&v);
    for(int i=1;i<=m;i++){
        int x,y;
        double a,b,c,d;
        scanf("%d%d%lf%lf%lf%lf",&x,&y,&a,&b,&c,&d);
        add(x,y,a,b);add(y,x,c,d);
    }
    if(spfa(s))printf("YES\n");
    else printf("NO\n");
    return 0;
}