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]$[l,mid]$ 的部分，那麼[mid+1,r]$[mid+1,r]$ 的點都是可以經過的，對[mid+1,r]$[mid+1,r]$ 做一遍floyd，然後分治[l,mid]$[l,mid]$ 。
2、對於[mid+1,r]$[mid+1,r]$ 的部分，對[l,mid]$[l,mid]$ 處理，然後分治[mid+1,r]$[mid+1,r]$ 。
3、當l=r$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;
}