【JZOJ 3808】道路值守

Description

Crossbell 自治州有着四通八达的现代化交通。时值独立庆典之际，随着来自周边国家旅客的日益增
多，犯罪行为也悄无声息开始滋长起来。
特别任务支援科的警察们从总部收到了关于调查伪装在游客中的犯罪分子的请求。通过调查，他们
得到了一张地图，记载了Crossbell 自治州内每一条道路的长度。
显然，为了减少犯罪行为被发现的可能性，犯罪分子总是会选择最短的路径来行动。为了方便安排
人手和推测犯罪分子采取的路线，他们希望得知任意两个地点之间，有多少条犯罪分子可能会选择的道
路。

Solution

首先肯定先跑一遍Floyd，

有一种很暴力的算法，枚举一个起点，选定一个终点，
再跑一遍全局，看看有多少条边在最短路上，
可惜会TLE（我就是这么打的QWQ），（复杂度：$O(n^2(n+m))$）

那就不这么暴力，加一个优化：
枚举一个点i，把所有点按到它的距离排序，
设$d_x$表示所有它连出去的边中，有多少条在最短路上，
再枚举一个点j和k，那么如果k在最短路上，它对答案的贡献就是$d_k$，
最后再加上自己的$d_j$，就是i,j的答案了。

复杂度：$O(n^2\log(n)+n^2(n+m))$（常数小(笑)）

Code

#include <iostream>
#include <cstdio>
#include <cstdlib>
#include <algorithm>
#include <cstring>
#define fo(i,a,b) for(int i=a;i<=b;i++)
#define min(q,w) ((q)<(w)?(q):(w))
using namespace std;
const int N=505;
int read(int &n)
{
    scanf("%d",&n);return n;
    char ch=' ';int q=0,w=1;
    for(;(ch!='-')&&((ch<'0')||(ch>'9'));ch=getchar());
    if(ch=='-')w=-1,ch=getchar();
    for(;ch>='0' && ch<='9';ch=getchar())q=q*10+ch-48;n=q*w;return n;
}
int n,m,Ti;
int a[N][N],b[N][N];
int f[N][N];
int ans[N];
struct qqww
{int x,v;} dis[N];
int z[N];
bool PX(qqww q,qqww w){return q.v<w.v;}
int main()
{
    int q,w,e;
    read(n),read(m);
    memset(f,60,sizeof(f));
    fo(i,1,m)read(q),read(w),a[q][++a[q][0]]=w,a[w][++a[w][0]]=q,
    b[q][a[q][0]]=b[w][a[w][0]]=f[q][w]=read(f[w][q]);
    fo(k,1,n)
        fo(i,1,n)if(i!=k)
        fo(j,1,n)if(i!=j&&j!=k)f[i][j]=min(f[i][j],f[i][k]+f[k][j]);
    fo(i,1,n)
        if(a[i][0]>2&&a[i][a[i][0]-1]==a[i][a[i][0]])a[i][0]--;
    fo(i,1,n-1)
    {
        f[i][i]=0;
        fo(j,1,n)dis[j].v=f[i][j],dis[j].x=j,z[j]=ans[j]=0;
        sort(dis+1,dis+1+n,PX);
        fo(j,2,n)if(dis[j].v<10000000)
        {
            q=dis[j].x;
            fo(k,1,j)if(dis[j].v==dis[k].v+f[dis[k].x][q])ans[q]+=z[k];
            fo(k,1,a[q][0])
                if(f[i][a[q][k]]+b[q][k]==dis[j].v)
                        z[j]++;
            ans[q]+=z[j];
        }
        fo(j,i+1,n)printf("%d ",ans[j]);
    }
    printf("\n");
    return 0;
}