[分治 Floyd] 计蒜之道2016复赛 .百度地图的实时路况

对 $y$ 分治

假设处理到 $[l,r]$ 区间，如果 $l=r$ ，$n^2$ 扫一遍 $dis$ 数组记录答案。

否则把 $[l,mid]$ 的边加入图中，递归做 $[mid+1,r]$，再把 $[mid+1,r]$ 的边加入图中，递归做 $[l,mid]$

#include <cstdio>
#include <iostream>
#include <algorithm>
#include <cstring>

using namespace std;

const int N=310,inf=1<<29;

typedef long long ll;

int n;
int dis[N][N];
ll ans;

void solve(int l,int r){
    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[i][j]==inf?-1:dis[i][j];
        return ;
    }
    int mid=l+r>>1,tmp[N][N];
    memcpy(tmp,dis,sizeof(tmp));
    for(int k=l;k<=mid;k++)
        for(int i=1;i<=n;i++)
            for(int j=1;j<=n;j++)
                dis[i][j]=min(dis[i][j],dis[i][k]+dis[k][j]);
    solve(mid+1,r);
    memcpy(dis,tmp,sizeof(dis));
    for(int k=mid+1;k<=r;k++)
        for(int i=1;i<=n;i++)
            for(int j=1;j<=n;j++)
                dis[i][j]=min(dis[i][j],dis[i][k]+dis[k][j]);
    solve(l,mid);
    memcpy(dis,tmp,sizeof(dis));
}

int main(){
    scanf("%d",&n);
    for(int i=1;i<=n;i++)
        for(int j=1;j<=n;j++){
            scanf("%d",&dis[i][j]);
            if(!~dis[i][j]) dis[i][j]=inf; 
        }
    solve(1,n);
    cout<<ans<<endl;
    return 0;
}