Ten Score Pretty‘s Bellman-Ford

好久没更新了，今天正好学图论，讲一讲（其实是因为没怎么写过图论）。

上周学的Dijsktra学的我是十分不爽，因为老师一上来就说他的弱点：不能有负权边。以至于什么上限高啊什么的，没兴趣，基础不打牢，绝对不是个好算法（我自己说的）

这一周，我就见到了一个人畜无害的新算法：Bellman-Ford。虽然上限一般般，可是地基打的牢固啊！

Bellman-Ford的思路大概是这样的：

1. 初始化:将除源点外的所有顶点的最短距离估计值 d[v] -->+∞, d[s]-->0;
2. 把从权值非无穷大的点能到的其他点（此处指一步）从这个点算花费，如果此花费小于原权值，则更新此点的权值。（专业一点叫松弛），迭代求解:反复对边集E中的每条边进行松弛操作，使得顶点集V中的每个顶点v的最短距离估计值逐步逼近其最短距离;(运行|v|-1次)，因为图的任意一条最短路径既不能包含负权回路，也不会包含正权回路，因此它最多包含|v|-1条边
3. 检验负权回路:判断边集E中的每一条边的两个端点是否收敛。如果存在未收敛的顶点，则算法返回false，表明问题无解;否则算法返回true，并且从源点可达的顶点v的最短距离保存在 d[v]中。

没什么特别的，思路简单，直接上代码。

#include<bits/stdc++.h>
using namespace std;
struct node{
    int from,to,len;
}a[15555];
int n,m,w,x,y,z,s,t;
int d[2505],cnt;
void add(int from,int to,int len){
    a[cnt].from=from,a[cnt].to=to,a[cnt].len=len;
    cnt++;
}
int main(){
    scanf("%d %d %d %d",&n,&m,&s,&t);
    for(int i=1;i<=n;i++) d[i]=1e9;
    d[s]=0;
    for(int i=0;i<m;i++){
        scanf("%d %d %d",&x,&y,&z);
        add(x,y,z);
    }
    int from,to,len;
    for(int k=0;k<n;k++){
        for(int i=0;i<cnt;i++){
            from=a[i].from,to=a[i].to,len=a[i].len;
            if(d[to]>d[from]+len) d[to]=d[from]+len;
        }
    }
    printf("%d\n",d[t]);
}

就很nice。（我个究极脑残都会写了！）

可是，一旦遇上了负权值，就离不开一个小东西：负权回路！

负权回路指的是：一个图里的一个小圈圈，这个圈圈整圈上的权值加起来的值是负的。之所以要考虑这个就是因为只要你一直不停的走这个圈圈，你的权值就这辈子都不会有最低。

所以，在程序里，我们还要排查负权回路，只要出现，就没有最小到达值，直接return掉。（所以负权回路做错了什么，要通缉他）

给大家参考一个排查法

#include<bits/stdc++.h>
using namespace std;
struct node{
    int from,to,len;
};
vector<node>v[505];
int n,m,w,x,y,z,d[505];
int main(){
    scanf("%d %d %d",&n,&m,&w);
    for(int i=2;i<=n;i++) d[i]=25000000;
    for(int i=0;i<m;i++){
        scanf("%d %d %d",&x,&y,&z);
        v[x].push_back((node){x,y,z});
        v[y].push_back((node){y,x,z});
    }
    for(int i=0;i<w;i++){
        scanf("%d %d %d",&x,&y,&z);
        v[x].push_back((node){x,y,-z});
    }
    int from,to,len;
    for(int k=0;k<n;k++){
        for(int i=1;i<=n;i++){
            for(int j=0;j<v[i].size();j++){
                from=v[i][j].from,to=v[i][j].to,len=v[i][j].len;
                if(d[to]>d[from]+len){
                	d[to]=d[from]+len;
                    if(k==n-1){
                        printf("YES\n");//有负权回路 
                        return 0;
                    }	
                }
            }
        }
    }
    printf("NO\n");//没有负权回路 
}

最后，再把这两个的思路合体：哇，金色传说！！！！！！！

#include<bits/stdc++.h>
using namespace std;
struct node{
    int from,to,len;
}a[5501];
int n,m,w,x,y,z,s,e,d[5001],cnt;
void add(int from,int to,int len){
    a[cnt].from=from,a[cnt].to=to,a[cnt].len=len;
    cnt++;
}
int main(){
    scanf("%d %d %d",&n,&m,&s);//有终点"前加一个 %d,并在s后面加,&e
    for(int i=1;i<=n;i++) d[i]=1e9;
    d[s]=0;
    for(int i=0;i<m;i++){
        scanf("%d %d %d",&x,&y,&z);
        add(x,y,z);//无向图后加add(y,x,z); 
    }
    int from,to,len;
    for(int k=0;k<n;k++){
        for(int i=0;i<cnt;i++){
            from=a[i].from,to=a[i].to,len=a[i].len;
            if(d[to]>d[from]+len){//边的松弛
            	d[to]=d[from]+len;
            	if(k==n-1){
            		printf("no solution\n");//有负权回路
            		return 0;
				}
			}
        }
    }
    for(int i=1;i<=n;i++){//每一个位置权值一一打印
    	if(d[i]==1e9) d[i]=-1;
    	printf("%d ",d[i]);
	}
	//有终点去上写：if(d[e]==1e9) d[e]=-1; \n printf("%d",d[e]);
	return 0;
}

好啦，点个关注，明天发一篇关于queue简化的Bellman-Ford。