Floyd 弗洛伊德算法 多源最短路

原理：动态规划

设中转点是ki，那么i->j的最短路是min(dis[i][j],dis[i][ki]+dis[ki][j])，就是他的最短路要么是两个点直接相连，要么是两个点经过一系列的中转点相连，一系列中转点是一个集合，可以是空集，也可以是k1,也可以是k1、k2，他可以是K的任意子集（K={k1,k2,......,kn}）由于最小值满足区间加法，因此我们可以用dp优化 时间复杂度O（N^3）

用途 

• 求多源最短路
• 判断负环

求任意两点的最短距离，如果有负环可以判断出来 

核心

bool floyd() {
	for(int k=1; k<=n; k++) //k必须放前面，否则比如说k=3、4,i=1,j=2，那就只有1->3->2,1->4->2,而1->3->4->2丢了
		//f[i][k](k-1)+f[k][j](k-1)<=f[i][j](k)，括号表示中转点你要丢掉k-1这一维度并保证无后效性，必须把k-1下所有情况都考虑到,所以k必须放到最外层
		for(int i=1; i<=n; i++) //先丢掉哪一维度，哪一维度套外面
			for(int j=1; j<=n; j++)
				if(dis[i][k]!=INF && dis[k][j]!=INF && dis[i][k]+dis[k][j]<dis[i][j]) {
					dis[i][j]=dis[i][k]+dis[k][j];//dp找到更短路径
				}
	for(int i=1;i<=n;i++)	if(dis[i][i]<0)	return 1;//有负环 
	return 0;//无负环
}

易错

 k必须在最外层，否则比如说k=3、4,i=1,j=2，那就只有1->3->2,1->4->2,而1->3->4->2丢了

f[i][k](k-1)+f[k][j](k-1)<=f[i][j](k)，括号表示中转。你要丢掉k-1这一维度并保证无后效性，必须把k-1下所有情况都考虑到，所以k必须放到最外层，先丢掉哪一维度，哪一维度套外面

模板

#include<cstdio>
const int N=210,INF=0x3f3f3f3f;
int dis[N][N];//表示最短距离
int n,m;//n个点，m条边
bool floyd() {
	for(int k=1; k<=n; k++) //k必须放前面，否则比如说k=3、4,i=1,j=2，那就只有1->3->2,1->4->2,而1->3->4->2丢了
		//f[i][k](k-1)+f[k][j](k-1)<=f[i][j](k)，括号表示中转点你要丢掉k-1这一维度并保证无后效性，必须把k-1下所有情况都考虑到,所以k必须放到最外层
		for(int i=1; i<=n; i++) //先丢掉哪一维度，哪一维度套外面
			for(int j=1; j<=n; j++)
				if(dis[i][k]!=INF && dis[k][j]!=INF && dis[i][k]+dis[k][j]<dis[i][j]) {
					dis[i][j]=dis[i][k]+dis[k][j];//dp找到更短路径
				}
	for(int i=1;i<=n;i++)	if(dis[i][i]<0)	return 1;//有负环 
	return 0;//无负环
}
void init() {
	for(int i=1; i<=n; i++)
		for(int j=1; j<=n; j++) {
			if(i==j)	dis[i][j]=0;
			else	dis[i][j]=INF;
		}
}
void print() {
	for(int i=1; i<=n; i++)
		for(int j=1; j<=n; j++) {
			if(dis[i][j]!=INF)	printf("the distance from %d to %d is %d\n",i,j,dis[i][j]);
		}
}
int main() {
	scanf(" %d%d",&n,&m);
	init();
	for(int i=0; i<m; i++) {
		int x,y,w;
		scanf(" %d%d%d",&x,&y,&w);
		dis[x][y]=w;
	}
	bool t=floyd();
	if(t)	printf("有负环\n");
	else print();
	return 0;
}
/*
4 4
1 2 1
2 3 2
3 4 3
4 1 -8
*/