计蒜客 百度地图的实时路况

Description

Solution

设 (solve(l,r)) 为不经过 ([l,r]) 之间的点的最短路
考虑 (floyd) 需要枚举一个 (k) 作为中转站,我们只需要保证 (k) 不是 ([l,r]) 之间的就可以了

具体的我们用 ([l,mid]) 更新完之后去处理 (solve(mid+1,r)),用 ([mid+1,r]) 更新完之后去处理 (solve(l,mid)),就可以答案这个效果了

#include<bits/stdc++.h>
using namespace std;
const int N=305,inf=707406378;
int n,dis[10][N][N];long long ans=0;
inline void solve(int l,int r,int t){
	if(l==r){
		for(int i=1;i<=n;i++)
			for(int j=1;j<=n;j++)
				if(i!=j && i!=l && j!=l)
					ans+=(dis[t][i][j]>=inf?-1:dis[t][i][j]);
		return ;
	}
	memset(dis[t+1],127/3,sizeof(dis[t+1]));
	int mid=(l+r)>>1;
	for(int i=1;i<=n;i++)
		for(int j=1;j<=n;j++)
			dis[t+1][i][j]=dis[t][i][j];
	for(int k=l;k<=mid;k++)
		for(int i=1;i<=n;i++)
			for(int j=1;j<=n;j++)
				dis[t+1][i][j]=min(dis[t+1][i][j],dis[t+1][i][k]+dis[t+1][k][j]);
	solve(mid+1,r,t+1);
	
	memset(dis[t+1],127/3,sizeof(dis[t+1]));
	for(int i=1;i<=n;i++)
		for(int j=1;j<=n;j++)
			dis[t+1][i][j]=dis[t][i][j];
	for(int k=mid+1;k<=r;k++)
		for(int i=1;i<=n;i++)
			for(int j=1;j<=n;j++)
				dis[t+1][i][j]=min(dis[t+1][i][j],dis[t+1][i][k]+dis[t+1][k][j]);
	solve(l,mid,t+1);	
}
int main(){
	freopen("pp.in","r",stdin);
	freopen("pp.out","w",stdout);
	scanf("%d",&n);
	for(int i=1;i<=n;i++)
		for(int j=1;j<=n;j++){
			scanf("%d",&dis[0][i][j]);
			if(dis[0][i][j]==-1)dis[0][i][j]=inf;
		}
	solve(1,n,0);
	cout<<ans<<endl;
	return 0;
}