2022 Hubei Provincial Collegiate Programming Contest

bigstrength

已于 2022-07-12 11:39:40 修改

阅读量909

点赞数

分类专栏： ICPC训练赛图论文章标签：算法贪心算法图论

于 2022-07-12 11:37:50 首次发布

本文链接：https://blog.csdn.net/apple_52197802/article/details/125739947

版权

图论同时被 2 个专栏收录

9 篇文章 0 订阅

订阅专栏

ICPC训练赛

2 篇文章 0 订阅

订阅专栏

2022 Hubei Provincial Collegiate Programming Contest

A.Nucleic Acid Test

题意

有n个城市，m条双向道路，n个城市中k个城市为核酸点，每次核酸的有效时间为t，要访问所有城市并且从核酸点出发并且以核酸点结束，且中途一直保证上次做过的核酸在有效期之内，问达到上述要求的最小的行进速度是多少。

数据范围

$\leqslant n \leqslant 300， 0 \leqslant m \leqslant \frac{n * (n - 1)}{2}, 1 \leqslant k \leqslant n， 0 \leqslant t \leqslant 10^9，1 \leqslant a_i,b_i\leqslant n，1 \leqslant c_i \leqslant 10^9$

考点

贪心， $f l o y d$ ， $K r u s k a l$ ， (二分)

思路

按照贪心的思路分类讨论，为了找到中途不做核酸走的最短距离 $s$ ，我们会想到把核酸点和非核酸点分开讨论：

对于非核酸点，我们是找到距离它们最近的核酸点 $x$ ，然后从 $x$ 出发，走到非核酸点然后回来，设所有非核酸点的距离各自 $x$ 的两倍距离中最大距离设为 $s_1$ 。
对于核酸点，我们只要对所有核酸点求一次最小生成树，然后取其中最大的那条边 $s_2$ 。
最后， $s = max(s_1, s_2)$
我们可以二分 $s$ ，答案是 $\frac{s}{t}$ 向上取整
答案为 $- 1$ 的情况：当所有的点不连通或者 $t = 0$ 的时候（坑点）

代码

#include<bits/stdc++.h>
using namespace std;
typedef long long LL;
const int N = 310, M = 100000;
int n, m, k, t;
int fa[N];
LL d[N][N];
bool st[N], q[N];
int findset(int x){
	if(fa[x] == x) return x;
	fa[x] = findset(fa[x]);
	return fa[x];
}
bool check(LL u){
	for(int i = 1; i <= n; i ++ ) fa[i] = i;
	memset(q, 0, sizeof q);
	int cnt = 0;
	for(int i = 1; i <= n; i ++ ){
		if(!st[i]) continue;
		q[i] = true;
		for(int j = 1; j <= n; j ++ ){
			if(st[j]){
				if(d[i][j] <= u){
					int fi = findset(i), fj = findset(j);
					if(fi != fj){
						fa[fi] = fj;
						cnt ++;
					}
				}
			}else{
				if(d[i][j] * 2 <= u){
					q[i] = q[j] = true;
				}
			}
		}
	}
	for(int i = 1; i <= n; i ++ ) if(!q[i]) return false;
	if(cnt != k - 1) return false;
	return true;
}
int main(){
	ios::sync_with_stdio(false);cin.tie(0);
	cin >> n >> m >> k >> t;
	if(t == 0){
		puts("-1");
		return 0;
	}
	for(int i = 1; i <= n; i ++ ) fa[i] = i;
	for(int i = 1; i <= n; i ++)
		for(int j = 1; j <= n; j ++ )
			if(i != j) d[i][j] = 1e18;
	for(int i = 1; i <= m; i ++ ){
		int a, b;
		LL w;
		cin >> a >> b >> w;
		d[a][b] = d[b][a] = min(d[a][b], w);
	}
	for(int p = 1; p <= n; p ++ ) for(int i = 1; i <= n; i ++ ) for(int j = 1; j <= n; j ++ ) d[i][j] = min(d[i][j], d[i][p] + d[p][j]);
	for(int i = 1; i <= k; i ++ ){
		int x;
		cin >> x;
		st[x] = true;
	}
	LL l = 1, r = 1e17;
	while(l < r){//最小的距离
		LL mid = l + r >> 1;
		if(check(mid)) r = mid;
		else l = mid + 1;
	}
	if(l == 1e17){
		 puts("-1");
		 return 0;
	}
	if(l % t == 0) l = l / t;
	else l = l / t + 1;
	cout << l << endl;
	return 0;
}