多源最短路（Floyd）

2.多源汇最短路 
2.1 Floyd              o(n^3) 

---

Floyd        o(n^3)

注意：可以有负权边，但不能有负权回路

int d[N][N];//邻接矩阵
void init()
{
	for(int i=1;i<=n;i++)
		for(int j=1;j<=n;j++)
			if(i==j) d[i][j]=0;
			else d[i][j]=INF;
}
void floyd()
{
	for(int k=1;k<=n;k++)
		for(int i=1;i<=n;i++)
			for(int j=1;j<=n;j++)
				d[i][j]=min(d[i][j],d[i][k]+d[k][j]);
}

---

牛的旅行

1125. 牛的旅行 - AcWing题库

题意：两个连通块之间加一条边，求新的连通块的直径的最小值

思路： 首先Floyd求出所有点之间的最短距离，计算每个点到所在连通块的点的最远距离maxd，然后两两枚举两个连通块之间的点i和j连线，直径可能为min(maxd[i]+maxd[j]+点i到点j的距离)，也可能为max(maxd[i])，两者较大的一个为直径

#include<bits/stdc++.h>
using namespace std;
const int N=155;
typedef pair<int,int> PII;
int n;
double g[N][N];
char w[N][N];
double maxd[N];
PII node[N]; 
double cal(int i,int j){
	int x1=node[i].first,y1=node[i].second;
	int x2=node[j].first,y2=node[j].second;
	return (double)sqrt((x1-x2)*(x1-x2)+(y1-y2)*(y1-y2));
}
int main(){
	cin>>n;
	for(int i=1;i<=n;i++){
		int x,y;cin>>x>>y;
		node[i]={x,y};
	}
	for(int i=1;i<=n;i++){
		cin>>w[i]+1;
	}
	for(int i=1;i<=n;i++){
		for(int j=1;j<=n;j++){
			if(i==j) g[i][j]=0;
			else if(w[i][j]=='1') g[i][j]=cal(i,j);
			else g[i][j]=1e20;
		}
	}
	for(int k=1;k<=n;k++){
		for(int i=1;i<=n;i++){
			for(int j=1;j<=n;j++){
				g[i][j]=min(g[i][j],g[i][k]+g[k][j]);
			}
		}
	}
	//点i在连通块内的最远距离 
	for(int i=1;i<=n;i++){
		for(int j=1;j<=n;j++){
			if(g[i][j]!=1e20){
				maxd[i]=max(maxd[i],g[i][j]);
			}
		}
	}
	//所有连通块的最大直径 
	double res1=0;
	for(int i=1;i<=n;i++) res1=max(res1,maxd[i]);
	//枚举在那两个不连通的点之间连边 
	//距离为 i到j的距离 + i所在连通块直径 + j所在连通块直径 
	double res2=1e20;
	for(int i=1;i<=n;i++){
		for(int j=1;j<=n;j++){
			if(g[i][j]==1e20){
				double d=cal(i,j)+maxd[i]+maxd[j];
				res2=min(res2,d);
			}
		}
	}
	cout<<fixed<<setprecision(6)<<max(res1,res2);
} 

---

 铁路与公路

 4074. 铁路与公路 - AcWing题库

        两个城市之间有长度为1的公路或铁路，汽车走公路，火车走铁路，同时从点1出发，且不能同时到达同一点，问最终都到达点n的最短时间

        发现性质：除了起点之外，汽车和火车到达某一点的最短时间不相等，故铁路和公路分别建图做Floyd即可

#include<bits/stdc++.h>
using namespace std;
const int N=405;
typedef pair<int,int> PII;
int n,m;
int g[N][N],f[N][N];
int floyd(int g[N][N]){
	for(int k=1;k<=n;k++){
		for(int i=1;i<=n;i++){
			for(int j=1;j<=n;j++){
				g[i][j]=min(g[i][j],g[i][k]+g[k][j]);
			}
		}
	}
	return g[1][n];
}
int main(){
	cin>>n>>m;
	memset(g,0x3f,sizeof g);
	memset(f,0x3f,sizeof f);	
	while(m--){
		int a,b;cin>>a>>b;
		g[a][b]=g[b][a]=1;
	}
	for(int i=1;i<=n;i++){
		for(int j=1;j<=n;j++){
			if(i==j) continue;
			if(g[i][j]==0x3f3f3f3f){
				f[i][j]=1;
			}
		}
	}
	int res;
	if(g[1][n]!=0x3f3f3f3f){
		res=floyd(f);
	}else{
		res=floyd(g);
	}
	cout<<(res==0x3f3f3f3f?-1:res);
} 

---

牛站

345. 牛站 - AcWing题库

思路：Floyd倍增

#include<bits/stdc++.h>
using namespace std;
typedef long long ll;
typedef pair<int,int> PII;
const int N=1e3+5;
int n,t,st,ed;
int g[N][N];
map<int,int> mp;
int cnt;
void floyd(int A[N][N],int B[N][N],int res[N][N]){
	int t[N][N];
	memset(t,0x3f,sizeof t);
	for(int k=1;k<=cnt;k++){
		for(int i=1;i<=cnt;i++){
			for(int j=1;j<=cnt;j++){
				t[i][j]=min(t[i][j],A[i][k]+B[k][j]);
			}
		}
	}
	memcpy(res,t,sizeof t);
}
void qmi(int res[N][N],int g[N][N],int k){
	while(k){
		if(k&1) floyd(res,g,res);
		floyd(g,g,g);
		k>>=1;
	}
}
int main(){
	cin>>n>>t>>st>>ed;
	//点离散化，否则会超时 
	mp[st]=++cnt;mp[ed]=++cnt;
	//g代表经过一条边后 各点之间的距离
	//所以g[i][i]不能初始化为0
	memset(g,0x3f,sizeof g);
	while(t--){
		int w,a,b;cin>>w>>a>>b;
		if(!mp[a]) mp[a]=++cnt;
		if(!mp[b]) mp[b]=++cnt;
		a=mp[a],b=mp[b];
		g[a][b]=min(g[a][b],w);
		g[b][a]=min(g[b][a],w);
	}
	int res[N][N];
	//初始res代表经过0条边 各点的距离 
	memset(res,0x3f,sizeof res);
    for(int i=1;i<=cnt;i++) res[i][i]=0;
	qmi(res,g,n);
	cout<<res[mp[st]][mp[ed]];
} 

---

最短路之和 

4872. 最短路之和 - AcWing题库

         n个点的图，n次操作，每次操作前统计剩下的点对之间最短路之和，然后将点x[i]删除，求每次操作前的最短路之和

        逆向思维，将正序删点改为倒序加点 

#include<bits/stdc++.h>
using namespace std;
typedef long long ll;
typedef pair<int,int> PII;
const int N=505,M=N<<1;
int g[N][N];
int n;
int x[N];
bool vis[N];
vector<ll> vc;
void floyd(){
	ll res=0;
	//正序删点 改为 倒序加点 
	for(int z=n;z>=1;z--){
		int k=x[z];
		vis[k]=true;//加入点k 
		res=0;
		for(int i=1;i<=n;i++){
			for(int j=1;j<=n;j++){
				g[i][j]=min(g[i][j],g[i][k]+g[k][j]);
			}
		}
		for(int i=1;i<=n;i++){
			if(!vis[i]) continue;
			for(int j=1;j<=n;j++){
				if(!vis[j]) continue;
				res+=g[i][j];
			}
		}
		vc.push_back(res);
	}
}
int main(){
	scanf("%d",&n);
	for(int i=1;i<=n;i++){
		for(int j=1;j<=n;j++){
			scanf("%d",&g[i][j]);
		}
	}
	for(int i=1;i<=n;i++){
		scanf("%d",&x[i]);
	}
	floyd();
	for(int i=vc.size()-1;i>=0;i--){
		printf("%lld ",vc[i]);
	}
}