最短路计数
题目描述
给出一个 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。
关于 S P F A ,它死了
大致思路
- 显而易见,这个题要跑最短路, D i j k s t r a Dijkstra Dijkstra 与 S P F A S \space P\space F \space A S P F A都可以(本蒟蒻写的 S P F A S \space P\space F \space A S P F A)
- 在 S P F A S \space P\space F \space A S P F A的松弛过程中进行最短路条数的计数,若当前点成功被松弛,则说明当前点被加入到最短路路径中来,它当前的路径数量继承过去,即 a n s [ v ] = a n s [ u ] ans[v]=ans[u] ans[v]=ans[u]
- 若目标点的最短路长度等于当前点的最短路长度+1,说明最短路径重叠,那么目标最短路长度就等于当前最短路数量加上目标最短路数量,即 d i s [ v ] = = d i s [ u ] + 1 dis[v]== dis[u]+1 dis[v]==dis[u]+1时, a n s [ v ] + = a n s [ u ] ans[v]+=ans[u] ans[v]+=ans[u]
AC CODE
#include<bits/stdc++.h>
using namespace std;
const int N=1e6+123;
#define int long long int
int n,m;
int dis[N],ans[N];
vector<int> a[N];
bool vis[N];
void SPFA(){
memset(vis,0,sizeof(vis));
memset(dis,0x7f,sizeof(dis));
memset(ans,0,sizeof(ans));
queue<int> q;
dis[1]=0;ans[1]=1;vis[1]=1;
q.push(1);
while(!q.empty()){
int x=q.front();
q.pop();
vis[x]=0;
for(int i=0;i<a[x].size();i++){
int y=a[x][i];
if(dis[y]>dis[x]+1){
dis[y]=dis[x]+1;
ans[y]=ans[x];
if(!vis[y]){
q.push(y);
vis[y]=1;
}
}
else if(dis[y]==dis[x]+1){
ans[y]+=ans[x];
ans[y]%=100003;
}
}
}
}
signed main(){
cin>>n>>m;
for(int i=1;i<=m;i++){
int x,y;
cin>>x>>y;
a[x].push_back(y);
a[y].push_back(x);
}
SPFA();
for(int i=1;i<=n;i++){
cout<<ans[i]%100003<<endl;
}
return 0;
}