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$

输出格式

输出文件包含两个数，分别是最少花费和花费最少的路径的总数。保证花费最少的路径的总数不超过 $2^{30}$。
两个不同的最短路方案要求：路径长度相同（均为最短路长度）且最短路经过的点的编号序列不同。若城市 $N$ 无法到达则只输出一个 No answer。

思路

首先去重，对于重边，可保留边权最小的那条：

• 如果边权大的那条边，显而易见无论对最短路还是计数都没有影响，可以省略。
• 如果存在边权相等且同为最小的边，则后续统计会出现重复，必须省略。

考虑Dijkstra算法，对于每次转移的时候，可以同时记录每个点的方案数，具体如下：

令 $wa{y}_i$ 表示到第 $i$ 个点的方案数，则：

• 从 $i$ 到 $j$ 存在新路径比原有到达 $j$ 路径要短： $wa{y}_j\leftarrow wa{y}_i$。
• 从 $i$ 到 $j$ 存在新路径同原来一样： $wa{y}_j\leftarrow wa{y}_j+wa{y}_i$。
  

代码

AC 记录

#include<bits/stdc++.h>
#define int long long
using namespace std;
int n,m;
int head[200005],nex[4000005],cnt = 0,a[2000005],way[2000005];
bool vis[200005];
struct edge{
	int to,w;
}node[4000005];
void add(int x,int y,int z) {
	nex[++cnt] = head[x];
	head[x] = cnt;
	node[cnt].to = y;
	node[cnt].w = z;
} 
struct point{
	int id,w;
	friend bool operator <(point x,point y) {
		return x.w > y.w;
	}
};
priority_queue<point>p; 
int ed[2005][2005];
signed main() {
	scanf("%lld %lld",&n,&m);
	for(int i = 1;i <= n;i++) a[i] = 2147483647;
	a[1] = 0,way[1] = 1;
	point b;
	b.id = 1,b.w = 0;
	p.push(b);
	for(int i = 1;i <= n;i++) for(int j = 1;j <= n;j++) ed[i][j] = 2147483647;
	for(int i = 1;i <= m;i++) {
		int x,y,z;
		scanf("%lld %lld %lld",&x,&y,&z);
		ed[x][y] = min(ed[x][y],z);
	}
	for(int i = 1;i <= n;i++) {
		for(int j = 1;j <= n;j++) {
			if(ed[i][j] != 2147483647) {
				add(i,j,ed[i][j]);
			}
		}
	}
	while(!p.empty()) {
		point now = p.top();
		p.pop();
		if(vis[now.id]) continue;
		vis[now.id] = 1;
		for(int i = head[now.id];i;i = nex[i]) {
			if(a[node[i].to] > a[now.id] + node[i].w) {
				a[node[i].to] = a[now.id] + node[i].w;
				point c;
				c.id = node[i].to;
				c.w = a[now.id] + node[i].w;
				p.push(c);
				way[node[i].to] = way[now.id];
			}
			else if(a[node[i].to] == a[now.id] + node[i].w) {
				way[node[i].to] += way[now.id];
			}
		}
	}
	//for(int i = 1;i <= n;i++) printf("%lld ",way[i]);
	//*
	if(a[n] != 2147483647) {
		printf("%lld %lld\n",a[n],way[n]);
	}
	else printf("No answer\n");
	//*/
	
    return 0;
}