最短路计数
题目描述
给出一个 N N N 个顶点 M M M 条边的无向无权图,顶点编号为 1 ∼ N 1\sim N 1∼N。问从顶点 1 1 1 开始,到其他每个点的最短路有几条。
输入格式
第一行包含 2 2 2 个正整数 N , M N,M N,M,为图的顶点数与边数。
接下来 M M M 行,每行 2 2 2 个正整数 x , y x,y x,y,表示有一条连接顶点 x x x 和顶点 y y y 的边,请注意可能有自环与重边。
输出格式
共 N N N 行,每行一个非负整数,第 i i i 行输出从顶点 1 1 1 到顶点 i i i 有多少条不同的最短路,由于答案有可能会很大,你只需要输出 $ ans \bmod 100003$ 后的结果即可。如果无法到达顶点 i i i 则输出 0 0 0。
样例 #1
样例输入 #1
5 7
1 2
1 3
2 4
3 4
2 3
4 5
4 5
样例输出 #1
1
1
1
2
4
提示
1 1 1 到 5 5 5 的最短路有 4 4 4 条,分别为 2 2 2 条 1 → 2 → 4 → 5 1\to 2\to 4\to 5 1→2→4→5 和 2 2 2 条 1 → 3 → 4 → 5 1\to 3\to 4\to 5 1→3→4→5(由于 4 → 5 4\to 5 4→5 的边有 2 2 2 条)。
对于
20
%
20\%
20% 的数据,
1
≤
N
≤
100
1\le N \le 100
1≤N≤100;
对于
60
%
60\%
60% 的数据,
1
≤
N
≤
1
0
3
1\le N \le 10^3
1≤N≤103;
对于
100
%
100\%
100% 的数据,
1
≤
N
≤
1
0
6
1\le N\le10^6
1≤N≤106,
1
≤
M
≤
2
×
1
0
6
1\le M\le 2\times 10^6
1≤M≤2×106。
思路
因为我们知道最短路至少可能有一条,而且我们又知道bfs的点入队仅一次(这其中出队就已经是最优解了),所以我们只需要在dist的被更新的时候当dist能被更新的话,说明什么?说明新的一条路不确定是最短路,因此就继承上一节点的cnt,反之,就加上(为什么呢?因为相等,证明新的一条路可能是最短路)。
- 细节1:记得用scanf和printf,要不然会超时。
- 细节2:最短路树一般都是拓扑图。
- 细节3:
代码
#include<iostream>
#include<algorithm>
#include<cstring>
#include<queue>
#define x first
#define y second
using namespace std;
const int N = 5e6+10,M = 2*N,mod=100003;
typedef pair<int,int>PII;
int e[M],ne[M],h[N],idx;
bool st[N];
int cnt[N];
int dist[N];
queue<int>q;
// priority_queue<PII,vector<PII>,greater<PII>>q;
int n,m;
void add(int a,int b){
e[idx]=b,ne[idx]=h[a],h[a]=idx++;
}
void bfs(){
memset(dist,0x3f,sizeof dist);
q.push(1);
dist[1]=0;
st[1]=true;
cnt[1]=1;//至少存在一条
while(q.size()){
int t=q.front();
q.pop();
for(int i=h[t];~i;i=ne[i]){
int j=e[i];
if(st[j])continue;
if(dist[j]>dist[t]+1){
dist[j]=dist[t]+1;
cnt[j]=cnt[t];//继承路径
q.push(j);
}else if(dist[j]==dist[t]+1){
cnt[j]=(cnt[j]+cnt[t])%mod;//加上之前的
//距离相等,说明存在两条距离相等的路径,所以条数就要累加
}
}
}
}
// void dijista(){
// memset(dist,0x3f,sizeof dist);
// dist[1]=0;
// cnt[1]=1;
// q.push({dist[1],1});
// while(q.size()){
// PII t=q.top();
// q.pop();
// int d=t.x,ver=t.y;
// if(st[ver])continue;
// st[ver]=true;
// for(int i=h[ver];~i;i=ne[i]){
// int j=e[i];
// if(dist[j]>d+1){
// dist[j]=d+1;
// cnt[j]=cnt[ver];
// q.push({dist[j],j});
// }else if(dist[j]==d+1){
// cnt[j]=(cnt[j]+cnt[ver])%mod;
// }
// }
// }
// }
int main(){
cin>>n>>m;
memset(h,-1,sizeof h);
for(int i=1;i<=m;i++){
int a,b;
scanf("%d%d",&a,&b);
add(a,b);
add(b,a);
}
bfs();
// dijista();
for(int i=1;i<=n;i++){
printf("%d\n",cnt[i]);
}
return 0;
}