洛谷 P1144 最短路计数 题解

原题地址

题目大意:给你一个无权无向图,求从1开始到每一个点的最短路有几条。
传说中的图论入门题啊…
bfs反着搜,珂以求出每个点到1的距离,第 i i i个点到1的距离用 d i s i dis_i disi表示。
f [ i ] f[i] f[i]为第 i i i个点的答案。
i i i个点的最短路一定是由 x ∣ x ∈ s o n [ i ] , d i s [ x ] + 1 = d i s [ i ] x|x\in son[i],dis[x]+1=dis[i] xxson[i],dis[x]+1=dis[i],这就好求了啊…
式子: f [ i ] = ∑ i = 1 s i z e i f [ s o n [ i ] ] ( 符 合 d i s [ s o n [ i ] ] + 1 = d i s [ i ] ) f[i]=\sum_{i=1}^{size_i}f[son[i]] (符合dis[son[i]]+1=dis[i]) f[i]=i=1sizeif[son[i]](dis[son[i]]+1=dis[i])
于是您的第一反应大多是从 1 − n 1-n 1n求值,可惜 W A \color{red}WA WA了。。。
为什么会 W A \color{red}{WA} WA呢?原因很简单,假如 d i s [ 6 ] = 2 , d i s [ 5 ] = 3 , 有 边 ( 5 , 6 ) dis[6]=2,dis[5]=3,有边(5,6) dis[6]=2,dis[5]=3,(5,6),在计算 f [ 5 ] f[5] f[5]时,其实是需要加上 f [ 6 ] f[6] f[6]的,但是由于 f [ 6 ] f[6] f[6]还没有求到,所以直接 W A \color{red}WA WA哭。。。
所以要根据深度来求。
深度从低到高求,就不会有漏求的情况了。
式子和前面一样, f [ 1 ] = 1 + ∑ i = 1 m e d g e ( i ) = = ( 1 , 1 ) f[1]=1+\sum_{i=1}^m edge(i)==(1,1) f[1]=1+i=1medge(i)==(1,1)
做完了, n i c e nice nice
C o d e \color{blue}Code Code

//unbfs
# include <bits/stdc++.h>
using namespace std;
struct edge
{
	int from,to;
	edge(){}
	edge(int u,int v) {from=u,to=v;}
};
vector <edge> g[1000000+10];
int n,m;
int dis[1000000+5];
bool vis[1000000+5];
int f[1000000+5];
vector <int> dep[1000005];
int maxl=INT_MIN;
void unbfs()
{
	queue <int> q;
	q.push(1) ;
	vis[1]=1;
	dis[1]=0;
	while(!q.empty())
	{
		int head=q.front() ;
		for(int i=0;i<g[head].size();i++) 
		{
			int nice=g[head][i].to;
			if(!vis[nice]) 
			{
				dis[nice]=dis[head]+1;
				vis[nice]=1;
				dep[dis[nice]].push_back(nice) ;
				maxl=max(maxl,dis[nice]) ;
				q.push(nice) ;
			}
		}
		q.pop() ;
	}
	return ;
}
void solve(int x) 
{
	for(int i=0;i<g[x].size();i++) 
	{
		int nice=g[x][i].to;
		if(dis[x]-dis[nice]!=1) continue ;
		f[x]+=f[nice];
		f[x]%=100003;
	}
	return ;
}
int main(void) 
{
	scanf("%d%d",&n,&m) ;
	f[1]=1;
	for(int i=1;i<=m;i++) 
	{
		int u,v;
		scanf("%d%d",&u,&v) ;
		g[u].push_back(edge(u,v)) ;
		g[v].push_back(edge(v,u)) ;
		if(u==v&&u==1) f[1]++;
	}
	unbfs() ;
	
	for(int i=1;i<=maxl;i++) 
	{
		for(int j=0;j<dep[i].size();j++) 
		{
			solve(dep[i][j]) ;
		}
	}
	for(int i=1;i<=n;i++) 
	{
		printf("%d\n",f[i]) ;
	}
	return 0;
}
  • 1
    点赞
  • 0
    收藏
    觉得还不错? 一键收藏
  • 1
    评论

“相关推荐”对你有帮助么?

  • 非常没帮助
  • 没帮助
  • 一般
  • 有帮助
  • 非常有帮助
提交
评论 1
添加红包

请填写红包祝福语或标题

红包个数最小为10个

红包金额最低5元

当前余额3.43前往充值 >
需支付:10.00
成就一亿技术人!
领取后你会自动成为博主和红包主的粉丝 规则
hope_wisdom
发出的红包
实付
使用余额支付
点击重新获取
扫码支付
钱包余额 0

抵扣说明:

1.余额是钱包充值的虚拟货币,按照1:1的比例进行支付金额的抵扣。
2.余额无法直接购买下载,可以购买VIP、付费专栏及课程。

余额充值