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

链接：https://nanti.jisuanke.com/t/11217

题意：中文题。

分析：首先我们会想到一个最暴力的做法，枚举v然后做Floyd，这样是O(n^4)的，时间不够。那我们仔细观察一下Floyd这个过程，我们会发现k时增量，我们每次用dis[i][k]+dis[k][j]来更新dis[i][j]。这里我们用分治处理这个k，分治处理k的区间[l,r]，那么我们就能将复杂度降到O(n^3*logn)。

代码：

#include<map>
#include<set>
#include<cmath>
#include<queue>
#include<bitset>
#include<math.h>
#include<vector>
#include<string>
#include<stdio.h>
#include<cstring>
#include<iostream>
#include<algorithm>
#pragma comment(linker, "/STACK:102400000,102400000")
using namespace std;
const int N=305;
const int mod=100000000;
const int MOD1=1000000007;
const int MOD2=1000000009;
const double EPS=0.00000001;
typedef long long ll;
const ll MOD=998244353;
const int INF=1000000010;
const ll MAX=1000000000000;
const double pi=acos(-1.0);
typedef double db;
typedef unsigned long long ull;
ll n;
ll ans=0,a

,dis

,now[10]

;
void cdq(int l,int r,int de) {
int i,j,k,mid=(l+r)>>1;
if (l==r) {
for (i=1;i<=n;i++)
for (j=1;j<=n;j++)
if (i!=l&&j!=l) {
if (dis[i][j]==MAX) ans-=1;
else ans+=dis[i][j];
}
return ;
}
for (i=1;i<=n;i++)
for (j=1;j<=n;j++) now[de][i][j]=dis[i][j];
for (k=l;k<=mid;k++)
for (i=1;i<=n;i++)
for (j=1;j<=n;j++)
dis[i][j]=min(dis[i][j],dis[i][k]+dis[k][j]);
cdq(mid+1,r,de+1);
for (i=1;i<=n;i++)
for (j=1;j<=n;j++) dis[i][j]=now[de][i][j];
for (k=mid+1;k<=r;k++)
for (i=1;i<=n;i++)
for (j=1;j<=n;j++)
dis[i][j]=min(dis[i][j],dis[i][k]+dis[k][j]);
cdq(l,mid,de+1);
}
int main()
{
int i,j;
scanf("%d", &n);
for (i=1;i<=n;i++)
for (j=1;j<=n;j++) {
scanf("%lld", &a[i][j]);
if (a[i][j]==-1) a[i][j]=MAX;
dis[i][j]=a[i][j];
}
cdq(1,n,0);
printf("%lld\n", ans);
return 0;
}