计蒜之道复赛A题 百度地图的实时路况（分治+floyd）

百度地图的实时路况功能相当强大，能方便出行的人们避开拥堵路段。一个地区的交通便捷程度就决定了该地区的拥堵情况。假设一个地区有 n 个观测点，编号从 1 到 n。定义 d(u,v,w) 为从 u号点出发，严格不经过 v 号点，最终到达 ww 号点的最短路径长度，如果不存在这样的路径，d(u,v,w) 的值为 -1。

那么这个地区的交通便捷程度 P 为：

$P = \sum_{1 \leq x,y,z \leq n , x \neq y , y \neq z}{d(x,y,z)}$

现在我们知道了该地区的 n 个点，以及若干条有向边，求该地区的交通便捷程度 P。

输入格式

第一行输入一个正整数 $n(4 \leq n \leq 300)$，表示该地区的点数。

接下来输入 n 行，每行输入 n 个整数。第 i 行第 j 个数$G_{i,j}(-1 \leq G_{i,j} \leq 10000;G_{i,i} = 0)$表示从 i 号点到 j 号的有向路径长度。如果这个数为 -1，则表示不存在从 i 号点出发到 j号点的路径。

输出格式

输出一个整数，表示这个地区的交通便捷程度。

样例输入

4
0 1 -1 -1
-1 0 1 -1
-1 -1 0 1
1 -1 -1 0
样例输出

4

---

题意如上，这题就是类似floyd，但是就是求不经过这个点的时候，所有点之间两两的最短路，正常方法应该是枚举不经过的点然后floyd，然后就是$O(n^4)$的复杂度，然后就GG了，但是考虑到每个不经过的点的时候，其实有很多的两点之间的距离并不会改变，就是有重复计算的，所以就是要想办法优化复杂度。
这里用到的是巧妙的分治，solve(l,r)表示当前要计算这个区间内的点不被经过时的最短路，然后考虑$[l,mid]$里这些点如果经过的话，那么可以用来计算$[mid+1,r]$这些区间里的答案，如果$[mid+1,r]$这些区间里的点经过的话，那么可以用来计算前半个区间的答案，考虑$[l,r]$的时候，其实前面的区间的点和后面区间的点都已经计算进去了，然后你进去前半段，就把后半段区间的点也算进去，如果考虑后半段，就把前面的点算进去，就是大概这么一个思路，其实实现起来还是比较麻烦的，逻辑上不太容易想清楚。
我仍需努力啊这种题就GG了。

---

代码：

#include <map>
#include <set>
#include <stack>
#include <queue>
#include <cmath>
#include <string>
#include <vector>
#include <cstdio>
#include <cctype>
#include <cstring>
#include <sstream>
#include <cstdlib>
#include <iostream>
#include <algorithm>
#pragma comment(linker, "/STACK:102400000,102400000")

using namespace std;
#define   MAX           200005
#define   MAXN          6005
#define   maxnode       15
#define   sigma_size    30
#define   lson          l,m,rt<<1
#define   rson          m+1,r,rt<<1|1
#define   lrt           rt<<1
#define   rrt           rt<<1|1
#define   middle        int m=(r+l)>>1
#define   LL            long long
#define   ull           unsigned long long
#define   mem(x,v)      memset(x,v,sizeof(x))
#define   lowbit(x)     (x&-x)
#define   pii           pair<int,int>
#define   bits(a)       __builtin_popcount(a)
#define   mk            make_pair
#define   limit         10000

//const int    prime = 999983;
const int    INF   = 0x3f3f3f3f;
const LL     INFF  = 0x3f3f;
const double pi    = acos(-1.0);
//const double inf   = 1e18;
const double eps   = 1e-8;
const LL    mod    = 1e9+7;
const ull    mx    = 133333331;

/*****************************************************/
inline void RI(int &x) {
    char c;
    while((c=getchar())<'0' || c>'9');
    x=c-'0';
    while((c=getchar())>='0' && c<='9') x=(x<<3)+(x<<1)+c-'0';
 }
/*****************************************************/

int dp[20][305][305];
LL ans;
int n;
void cdq(int num,int l,int r){
    if(l==r){
        for(int i=1;i<=n;i++){
            for(int j=1;j<=n;j++){
                if(i!=l&&j!=l&&i!=j){
                    if(dp[num][i][j]>=INF) ans-=1;
                    else ans+=dp[num][i][j];
                }
            }
        }
        return;
    }
    for(int i=1;i<=n;i++){
        for(int j=1;j<=n;j++) dp[num+1][i][j]=dp[num][i][j];
    }
    int mid=(l+r)/2;
    for(int k=l;k<=mid;k++){
        for(int i=1;i<=n;i++){
            if(dp[num+1][i][k]==-1) continue;
            for(int j=1;j<=n;j++){
                if(dp[num+1][k][j]==-1) continue;
                dp[num+1][i][j]=min(dp[num+1][i][j],dp[num+1][i][k]+dp[num+1][k][j]);
            }
        }
    }
    cdq(num+1,mid+1,r);
    for(int i=1;i<=n;i++){
        for(int j=1;j<=n;j++) dp[num+1][i][j]=dp[num][i][j];
    }
    for(int k=mid+1;k<=r;k++){
        for(int i=1;i<=n;i++){
            if(dp[num+1][i][k]==-1) continue;
            for(int j=1;j<=n;j++){
                if(dp[num+1][k][j]==-1) continue;
                dp[num+1][i][j]=min(dp[num+1][i][j],dp[num+1][i][k]+dp[num+1][k][j]);
            }
        }
    }
    cdq(num+1,l,mid);
}
int main(){
    cin>>n;
    ans=0;
    mem(dp,INF);
    for(int i=1;i<=n;i++){
        for(int j=1;j<=n;j++){
            scanf("%d",&dp[0][i][j]);
            if(dp[0][i][j]==-1) dp[0][i][j]=INF;
        }
    }
    cdq(0,1,n);
    cout<<ans<<endl;
    return 0;
}