4246. 【五校联考6day2】san

Description

小明经常去N 个地点，其中有些地点之间有直接的无向道路（共M 条这样的道路），可以直接互相到达，这些道路的长短不一。由于小明对这些道路都很熟悉，无论起点和终点在哪里，总能走最短路。小明有严重的强迫症，认为奇数很不和谐，如果他某一天从一个地点去另一个地点走过的路程是奇数，就会很不爽，但他又不想白白多走路，所以遇到最短路长度是奇数的情况就只能忍了。
如果从某个地点A 到另一个地点B 的最短路径长度为奇数，则称这条最短路径为“不和谐最短路”。如果一条不和谐最短路上包含地点C，则称它为“经过C 的不和谐最短路”。现在请你编程求出对于每个地点，经过它的不同的不和谐最短路数量。两条最短路不同，当且仅当它们途径的地点的序列不同。

Input

第一行两个正整数N;M，含义见题面。
接下来M 行，每行三个正整数Ai;Bi;Li，表示一条无向道路的两端和长度。

Output

　N 行每行一个整数，第i 行表示经过第i 个点的不同的不和谐最短路条数。

Sample Input

4 4
1 4 1
1 2 1
3 4 100
2 3 2

Sample Output

6
4
2
2

样例说明

长度为奇数的最短路有:1 → 2; 1 → 2 → 3; 1 → 4; 2 → 1; 3 → 2 → 1; 4 → 1。
这些路径中四个点的经过次数分别为6, 4, 2, 2。
其它一些路，如1 → 4 → 3 不是最短路，2 → 3 是最短路但长度为2，是偶数。这些路都不计入答案。

Data Constraint

对于50% 的数据，N ≤ 100；
对于全部数据，N ≤ 1000;M ≤ 3000，每条路的长度不超过1000。
保证图连通，无自环重边。

Solution

50%：N = 100，可以枚举起点和终点，用各种算法枚举两点间所有最短路并累加 答案，复杂度一般为 O(N^3 )。

100%：题目中的图是一般图，结构复杂没有规律。考虑枚举起点并计算单源最短路，保留所有最短路中的边（有向），原图就变成了一个DAG，可以很方便地在上面进行拓扑排序，DP 等。再考虑长度为奇数的最短路，从一个点断开，一定是 一边长度为奇数，另一边长度为偶数。因此可以在DAG上DP计算从起点到 达每个点的奇偶最短路径分别有几条，再逆序 DP计算从每个点出发的奇偶最1短路径条数，在后面的DP过程中顺便统计答案。这样对每个起点算一遍。复杂度O(NM log N)。

然而我的做法就是枚举起点i，求出i到每个点的dis[]，然后就dfs走那些是最短路的边，不是的可以忽略，这样原无向图就变成了一个DAG(有向无环图)。然后我们就计算出每个点的size，答案就加上size[i]，表示从起点S走到当前这个点以及所有这个点的儿子路程和为奇数的点都要经过当前点。那么如何计算呢？若果从起点开始走到当前的点走过的路程是奇数，就将size[x]=1，否则=0。然后size[x]+=size[son_x]。注意size[S(起点)]=0。

Code

#include<cstdio>
#include<cstring>
#include<algorithm>
#define N 1010
#define M 3010
using namespace std;
int n,m,x,y,z,f[N],bz[N],d[N*10],t[M*2],nx[M*2],c[M*2],ls[N],ans[N],sz[N];
void add(int x,int y,int z){
	t[++t[0]]=y;
	c[t[0]]=z;
	nx[t[0]]=ls[x];
	ls[x]=t[0];
}
void dg(int x,int st){
	if(st%2==1) sz[x]=1;else sz[x]=0;
	for(int k=ls[x],son;k;k=nx[k]){
		if(f[son=t[k]]==f[x]+c[k]){
			dg(son,st+c[k]);
			sz[x]+=sz[son];
		}
	}
	ans[x]+=sz[x];
}
int main(){
	scanf("%d%d",&n,&m);
	for(int i=1;i<=m;i++){
		scanf("%d%d%d",&x,&y,&z);
		add(x,y,z);add(y,x,z);
	}
	for(int i=1;i<=n;i++){
		memset(f,127,sizeof(f));
		memset(sz,0,sizeof(sz));
		f[d[1]=i]=0;
		int l=0,r=1;
		while(l++<=r){
			int x=d[l];bz[x]=0;
			for(int k=ls[x];k;k=nx[k]){
				int y=t[k];
				if(f[y]>f[x]+c[k]){
					f[y]=f[x]+c[k];
					if(!bz[y]){
						bz[y]=1;
						d[++r]=y;
					}
				}
			}
		}
		dg(i,0);
	} 
	for(int i=1;i<=n;i++) printf("%d\n",ans[i]);
	return 0;
}

---

---

作者：zsjzliziyang 
QQ:1634151125 
转载及修改请注明 
本文地址:https://blog.csdn.net/zsjzliziyang/article/details/86696013