对 y 分治
假设处理到
否则把 [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;
}