CF416E President‘s Path 题解

最新推荐文章于 2024-09-25 13:54:13 发布

MingJunYi

最新推荐文章于 2024-09-25 13:54:13 发布

阅读量659

点赞数 16

分类专栏：题解文章标签： floyd 信息竞赛算法 c++ 最短路

本文链接：https://blog.csdn.net/qq_49785217/article/details/142282858

版权

题解专栏收录该内容

29 篇文章 0 订阅

订阅专栏

CF416E President’s Path

题面：

题面翻译

对于 $n$ 个点 $m$ 条边的简单无向图（无重边、自环），记 $f (s, t)$ 为 $R_{s,t}$ 类边的数量。 $R_{s,t}$ 类边是指这样一类边：存在一条从 $s$ 到 $t$ 的最短路，满足这条边在最短路上。

对于每一对 $s, t$ ，你都要计算出 $f (s, t)$ 的值。

输入格式

先是 $n,m(2\le n\le 500,0\le m\le \frac{n(n-1)}{2})$ ，含义见题面。

接下来 $m$ 行表述边 $x_i,y_i,l_i)$ ，分别为连接的两个点的编号，边的长度。

输出格式

仅输出一行，包含 $\frac{n(n-1)}{2}$ 个数。前 $n - 1$ 个数为 $f(1,2),f(1,3),\dots,f(1,n)$ ；接下来 $n - 2$ 个数为 $f(2,3),f(2,4),\dots,f(2,n)$ ；以此类推。

样例 #1

样例输入 #1
5 6
1 2 1
2 3 1
3 4 1
4 1 1
2 4 2
4 5 4
样例输出 #1
1 4 1 2 1 5 6 1 2 1

观察 $n$ 的数据范围，自然想到跑 $f l oy d$ .

我们如果枚举每一条边，判断是否在最短路上，

对于最短路边 $(u, v)$ ：

$dis_{i,j} = dis_{i,u} + dis_{v,j} + w(u,v)$

时间复杂度是 $\mathcal{O}(n^4)$ 不能接受.

考虑优化到 $\mathcal{O}(n^3)$ ,

枚举点是否在最短路上，

对于最短路点 $k$ ：

$dis_{i,j} = dis_{i,k} + dis_{k,j}$

令 $dp_{i,j}$ 为 $\rightsquigarrow j$ 上最短边的个数，

根据最短路（一条链）的性质，对于 $\rightsquigarrow k \rightsquigarrow j$ ，我们考虑以计算 $(x, k)$ 边的个数，累加到 $dp_{i,j}$ 。

我们如何计算 $(x, k)$ 呢？

我们将边备份给 $g_{i,j}$ ,然后对于 $i\rightsquigarrow k$ 计算 $\rightsquigarrow x \rightarrow k$ 的个数，用：

$dis_{i,k} = dis_{i,x} + g_{x,k}$

判断。

然后枚举起点 $i$ ，以上述方法更新 $dp_{i,j}$

AC-code:

#include<bits/stdc++.h>
using namespace std;
#define int long long
int rd() {
	int x = 0, w = 1;
	char ch = 0;
	while (ch < '0' || ch > '9') {
		if (ch == '-') w = -1;
		ch = getchar();
	}
	while (ch >= '0' && ch <= '9') {
		x = x * 10 + (ch - '0');
		ch = getchar();
	}
	return x * w;
}

void wt(int x) {
	static int sta[35];
	int f = 1;
	if(x < 0) f = -1,x *= f;
	int top = 0;
	do {
		sta[top++] = x % 10, x /= 10;
	} while (x);
	if(f == -1) putchar('-');
	while (top) putchar(sta[--top] + 48);
}
constexpr int N = 505,inf = 0x3f3f3f3f3f3f3f3fLL;
int n,m,f[N][N],g[N][N],dp[N][N],con[N];

signed main() {
	n = rd(),m = rd();
	for(int i = 1;i<=n;i++)
		for(int j = 1;j<=n;j++)
			f[i][j] = f[j][i] = g[i][j] = g[j][i] = inf;
	for(int i = 1;i<=m;i++) {
		int u = rd(),v = rd(),w = rd();
		g[u][v] = g[v][u] = f[u][v] = f[v][u] = min(w,f[u][v]);
	}
	for(int i = 1;i<=n;i++) f[i][i] = 0;
	for(int k = 1;k<=n;k++)
		for(int i = 1;i<=n;i++)
			for(int j = 1;j<=n;j++)
				f[j][i] = f[i][j] = min(f[i][j],f[i][k] + f[k][j]);
	for(int i = 1;i<=n;i++)  {
		for(int j = 1;j<=n;j++) con[j] = 0;
		for(int j = 1;j<=n;j++)
			for(int k = 1;k<=n;k++)
				if(f[i][k] + g[k][j] == f[i][j]) 
					con[j]++;
		for(int j = 1;j<=n;j++)
			for(int k = 1;k<=n;k++)
				if(f[i][k] + f[k][j] == f[i][j])
					dp[i][j] += con[k];
	}
	for(int i = 1;i<=n;i++)
		for(int j = i + 1;j<=n;j++) {
			if(f[i][j] >= inf) dp[i][j] = 0;
			wt(dp[i][j]),putchar(' ');
		}
	
					
	return 0;
}