【洛谷】P1608 路径统计

题目地址：

https://www.luogu.com.cn/problem/P1608

题目描述：
“RP餐厅”的员工素质就是不一般，在齐刷刷的算出同一个电话号码之后，就准备让HZH和TZY去送快餐了，他们将自己居住的城市画了一张地图，已知在他们的地图上，有$N$个地方，而且他们目前处在标注为$1$的小镇上，而送餐的地点在标注为$N$的小镇。（有点废话）除此之外还知道这些道路都是单向的，从小镇$I$到$J$需要花费$D[I,J]$的时间，为了更高效快捷的将快餐送到顾客手中，他们想走一条从小镇$1$到小镇$N$花费最少的一条路，但是他们临出发前，撞到因为在路上堵车而生气的FYY，深受启发，不能仅知道一条路线，万一。。。于是，他们邀请FYY一起来研究起了下一个问题：这个最少花费的路径有多少条？

输入格式：
输入文件第一行为两个空格隔开的数$N$，$E$，表示这张地图里有多少个小镇及有多少边的信息。下面$E$行，每行三个数$I$、$J$、$C$，表示从$I$小镇到$J$小镇有道路相连且花费为$C$（注意，数据提供的边信息可能会重复，不过保证$I\ne J,1\le I,J\le N$）。

输出格式：
输出文件包含两个数，分别是最少花费和花费最少的路径的总数。两个不同的最短路方案要求：路径长度相同（均为最短路长度）且至少有一条边不重合。若城市$N$无法到达则只输出一个No answer。

数据范围：
对于$30$的数据$N\le 20$；对于$100$的数据 $1\le N\le 2000,0\le E\le N\times (N-1),1\le C\le 10$。

注意这题的”边的信息可能会重复“这句话，它的意思是重复的边按最短的一条算，别的边都要当作不存在。所以用邻接矩阵存图更好。

思路是Dijkstra算法。注意到图是正权图，所以Dijkstra算最短路的时候，在最短路树上（当然不一定真的是棵树，有可能某个点的最短路可以由两个父亲走到，但是仍然可以将多的边加入成为一个拓扑图），点出堆的顺序实际上是一个拓扑序，所以整个Dijkstra算法可以在求出每个点最短路的同时，利用动态规划的递推思想，同时更新其能到达的点的最短路条数（如果找到了更短路，则直接赋值路径条数，如果找到了等长的最短路，则累加路径条数）。代码如下：

#include <iostream>
#include <cstring>
#include <queue>
using namespace std;

typedef pair<int, int> PII;

const int N = 2010, INF = 0x3f3f3f3f;
int n, m;
// 邻接矩阵存图
int g[N][N];
int dist[N], cnt[N];
bool st[N];

void dijkstra(int s) {
    memset(dist, 0x3f, sizeof dist);
    priority_queue<PII, vector<PII>, greater<PII> > heap;
    heap.push({0, s});

    dist[s] = 0;
    cnt[s] = 1;

    while (heap.size()) {
        auto t = heap.top();
        heap.pop();
        int v = t.second, d = t.first;

        if (st[v]) continue;
        st[v] = true;

        for (int i = 1; i <= n; i++) {
        	// 略过到不了的点，以及之前已经算过的点
            if (g[v][i] == INF || st[i]) continue;

            if (dist[i] > d + g[v][i]) {
                dist[i] = d + g[v][i];
                cnt[i] = cnt[v];
                heap.push({dist[i], i});
            } else if (dist[i] == d + g[v][i]) {
                cnt[i] += cnt[v];
            }
        }
    }
}

int main() {
    cin >> n >> m;
    memset(g, 0x3f, sizeof g);
    while (m--) {
        int a, b, c;
        cin >> a >> b >> c;
        // 有重边则取最短边
        g[a][b] = min(g[a][b], c);
    }

    dijkstra(1);

    if (dist[n] == INF) printf("No answer\n");
    else printf("%d %d\n", dist[n], cnt[n]);

    return 0;
}

时间复杂度$O(E\log N)$，空间$O(N)$。