2016计蒜之道复赛A 百度地图的实时路况 CDQ分治

题目在这里呀！

题意

定义d(u,v,w)为从u号点出发，严格不经过v号点，最终到达w号点的最短路径长度，如果不存在这样的路径，d(u,v,w)的值为-1。
计算每个d(u,v,w)的和。

题解

一开始没什么想法啊，枚举严格不经过的点然后做floyd复杂度O(n^4)，为什么我感觉用dijkstra可以？？
不管不管还是用正经做法，这题可以用CDQ分治解决，l,r表示最后严格不经过的点的区间，然后就当然分成[l,mid]和[mid+1,r]做呀。
1、对于$[l,mid]$的部分，那么$[mid+1,r]$的点都是可以经过的，对$[mid+1,r]$做一遍floyd，然后分治$[l,mid]$。
2、对于$[mid+1,r]$的部分，对$[l,mid]$处理，然后分治$[mid+1,r]$。
3、当$l=r$时，即这个点就是严格不能经过的点，统计答案即可。

//Suplex
#include <cstdio>
#include <iostream>
#include <cstring>
#include <algorithm>
#include <cmath>
#define N 500
#define ll long long
using namespace std;
int n,mat[N][N];

ll solve(int l,int r)
{
    ll ans=0;
    if(l==r){
        for(int i=1;i<=n;i++)
            for(int j=1;j<=n;j++)
                if(i!=r && j!=r) ans+=mat[i][j];
        return ans;
    }
    int tmp[N][N];
    memcpy(tmp,mat,sizeof(tmp));
    int mid=(l+r)>>1;
    for(int k=mid+1;k<=r;k++)
        for(int i=1;i<=n;i++)
            for(int j=1;j<=n;j++)
                if(i!=k && j!=k && i!=j && ~mat[i][k] && ~mat[k][j]){
                    if(mat[i][j]==-1) mat[i][j]=mat[i][k]+mat[k][j];
                    else mat[i][j]=min(mat[i][j],mat[i][k]+mat[k][j]);
                }
    ans+=solve(l,mid);
    memcpy(mat,tmp,sizeof(mat));
    for(int k=l;k<=mid;k++)
        for(int i=1;i<=n;i++)
            for(int j=1;j<=n;j++)
                if(i!=k && j!=k && i!=j && ~mat[i][k] && ~mat[k][j]){
                    if(mat[i][j]==-1) mat[i][j]=mat[i][k]+mat[k][j];
                    else mat[i][j]=min(mat[i][j],mat[i][k]+mat[k][j]);
                }
    ans+=solve(mid+1,r);
    memcpy(mat,tmp,sizeof(mat));
    return ans;
}

int main()
{
    scanf("%d",&n);
    for(int i=1;i<=n;i++)
        for(int j=1;j<=n;j++) scanf("%d",&mat[i][j]);
    printf("%lld\n",solve(1,n));
    return 0;
}