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

https://nanti.jisuanke.com/t/A1108

way1:

应该很多同学的做法都是对于每次y，每次x，dijkstra+堆优化

n^2 * nlogn

其实log(300)很小。。。

way2:

题解方法真心优秀！！！

cdq分治

有助于理解floyd：每次加入一个可以使用的点，加n次，顺序可以调整

本题：对于每个y，仅仅是y点不能使用

每次计算1/2，cdq剩下的1/2(left,right分别来一次)，直到只剩下一个数，这个点不用于求最短路。所有点都有且仅有一次出现“只剩下该点”的情况。

#include <cstdio>
#include <cstdlib>
#include <cmath>
#include <cstring>
#include <string>
#include <algorithm>
#include <iostream>
using namespace std;
#define ll long long

const double eps=1e-8;
const int inf=1e9;
const ll mod=1e9+7;
const int maxn=3e2+10;

int dist[maxn][maxn];
ll sum;
int n;

void cdq(int l,int r)
{
    int i,j;
    if (l==r)
    {
        for (i=1;i<=n;i++)
            for (j=1;j<=n;j++)
                if (i!=l && j!=l)
                    sum+=dist[i][j]==inf?-1:dist[i][j];
        return;
    }
    int m=(l+r)>>1,k,temp[maxn][maxn];
    memcpy(temp,dist,sizeof(dist));
    for (k=l;k<=m;k++)
        for (i=1;i<=n;i++)
            for (j=1;j<=n;j++)
                if (dist[i][j]>dist[i][k]+dist[k][j])
                    dist[i][j]=dist[i][k]+dist[k][j];
    cdq(m+1,r);

    memcpy(dist,temp,sizeof(temp));
    for (k=m+1;k<=r;k++)
        for (i=1;i<=n;i++)
            for (j=1;j<=n;j++)
                if (dist[i][j]>dist[i][k]+dist[k][j])
                    dist[i][j]=dist[i][k]+dist[k][j];
    cdq(l,m);
}

int main()
{
    int i,j;
    scanf("%d",&n);
    for (i=1;i<=n;i++)
        for (j=1;j<=n;j++)
        {
            scanf("%d",&dist[i][j]);
            if (dist[i][j]==-1)
                dist[i][j]=inf;
        }
    cdq(1,n);
    printf("%lld",sum);
    return 0;
}