货车运输(最大生成树+倍增LCA)
看到第一篇题解的神奇码风……我决定发一篇码风正常的题解造福人类
这题的做法也非常经典,最大生成树\(+LCA\),相当于先贪心一下,在LCA的时候记录一下当前最小的边权
顺便吐槽一下最后一个测试点:
testdata.in
7 8
1 2 2
1 3 5
3 4 4
4 4 2
3 5 3
6 7 4
1 3 3
4 5 8
8
1 2
1 4
1 3
1 5
1 6
2 5
3 5
6 7
testdata.out
2
4
5
4
-1
2
4
4
回到题面:注意: \(x\)不等于\(y\),两座城市之间可能有多条道路 。
虽然\(Kruskal\)不会挂掉……但是出题人造数据真不严谨
\(Code\ Below:\)
#include <bits/stdc++.h>
#define INF 99999999
using namespace std;
const int maxn=10000+10;
const int maxm=50000+10;
int n,m,q,head[maxn],f[maxn],tot;
//如题目所示,存储链式前向星和并查集
int fa[maxn][16],w[maxn][16],dep[maxn],rt,t;
//fa[i][j]表示第i个结点向上跳2^i个结点,dep[i]存储第i个结点的深度
struct Edge{
int x,y,w;
}e[maxm];
struct node{
int to,next,val;
}tree[maxm<<1];
inline bool cmp(Edge a,Edge b){//cmp
return a.w>b.w;
}
inline int read(){//读入优化
int x=0,f=1;char ch=getchar();
while(ch<'0'||ch>'9'){if(ch=='-')f=-1;ch=getchar();}
while(ch>='0'&&ch<='9'){x=(x<<3)+(x<<1)+ch-'0';ch=getchar();}
return x*f;
}
inline void swap(int &a,int &b){int t=a;a=b;b=t;}//交换
inline void add(int x,int y,int w){//存边
tree[++tot].to=y;
tree[tot].val=w;
tree[tot].next=head[x];
head[x]=tot;
}
int find(int x){//并查集
if(x!=f[x])
f[x]=find(f[x]);
return f[x];
}
inline void Kruskal(){//最大生成树
for(int i=1;i<=n;i++)
f[i]=i;
sort(e+1,e+m+1,cmp);
int ans=0;
for(int i=1;i<=m;i++){
int a=find(e[i].x),b=find(e[i].y);
if(a!=b){
f[a]=b;ans++;
add(e[i].x,e[i].y,e[i].w);
add(e[i].y,e[i].x,e[i].w);
}
}
}
void dfs(int x,int pre,int val){//LCA预处理
dep[x]=dep[pre]+1;
w[x][0]=val;
fa[x][0]=pre;
for(int i=1;i<=t;i++){
fa[x][i]=fa[fa[x][i-1]][i-1];
w[x][i]=min(w[x][i-1],w[fa[x][i-1]][i-1]);
}
for(int i=head[x];i;i=tree[i].next){
if(tree[i].to!=pre){
dfs(tree[i].to,x,tree[i].val);
}
}
}
inline int LCA(int x,int y){//倍增LCA
if(dep[x]<dep[y]) swap(x,y);
int ans=INF;
for(int i=t;i>=0;i--)
if(dep[fa[x][i]]>=dep[y]){
ans=min(ans,w[x][i]);
x=fa[x][i];
}
if(x==y) return ans;
for(int i=t;i>=0;i--)
if(fa[x][i]!=fa[y][i]){
ans=min(ans,min(w[x][i],w[y][i]));
x=fa[x][i],y=fa[y][i];
}
int val=w[x][0];
if(fa[x][0]!=y) val=min(val,w[y][0]);
return min(ans,val);
}
int main()
{
memset(fa,INF,sizeof(fa));
memset(w,INF,sizeof(w));
n=read(),m=read();t=log2(n)+1;
for(int i=1;i<=m;i++){
e[i].x=read(),e[i].y=read(),e[i].w=read();
}
Kruskal();
//记住!数据不一定联通!图是一个森林!
for(int i=1;i<=n;i++)
if(i==f[i])
dfs(i,i,INF);
scanf("%d",&q);
while(q--){
int x=read(),y=read();
//printf("%d %d\n",find(x),find(y));
if(find(x)!=find(y)) printf("-1\n");
else printf("%d\n",LCA(x,y));
}
return 0;
}
货车运输(最大生成树+倍增LCA)的更多相关文章
- TZOJ 4848 货车运输(最大生成树+倍增lca)
描述 A 国有 n 座城市,编号从 1 到 n,城市之间有 m 条双向道路.每一条道路对车辆都有重量限制,简称限重.现在有 q 辆货车在运输货物,司机们想知道每辆车在不超过车辆限重的情况下,最多能运多 ...
- $Noip2013/Luogu1967$ 货车运输 最大生成树+倍增$lca$
$Luogu$ $Sol$ 首先当然是构建一棵最大生成树,然后对于一辆货车的起点和终点倍增跑$lca$更新答案就好.记得预处理倍增的时候不仅要处理走了$2^i$步后是那个点,还有这中间经过的路径权值的 ...
- 【NOIP2013】货车运输 最大生成树+倍增
题目大意:给你一张n个点m条边的图,有q次询问,每次让你找出一条从x至y的路径,使得路径上经过的边的最小值最大,输出这个最大的最小值. 显然,经过的路径必然在这张图的最大生成树上. 我们求出这个图的最 ...
- 「NOIP2013」「LuoguP1967」货车运输(最大生成树 倍增 LCA
题目描述 AA国有nn座城市,编号从 11到nn,城市之间有 mm 条双向道路.每一条道路对车辆都有重量限制,简称限重.现在有 qq 辆货车在运输货物, 司机们想知道每辆车在不超过车辆限重的情况下,最 ...
- 【洛谷1967】货车运输(最大生成树+倍增LCA)
点此看题面 大致题意: 有\(n\)个城市和\(m\)条道路,每条道路有一个限重.多组询问,每次询问从\(x\)到\(y\)的最大载重为多少. 一个贪心的想法 首先,让我们来贪心一波. 由于要求最大载 ...
- 货车运输 noip2013 luogu P1967 (最大生成树+倍增LCA)
luogu题目传送门! 首先,题目让我们求每个货车的最大运输量,翻译一下就是求路径上边权最小的边. 利用一下贪心思想可知,所有货车肯定都会尽量往大的边走. 进一步翻译,即为有一些小边货车根本不会走,或 ...
- NOIP2013D1T3货车运输(最大生成树+倍增lca)
传送门 这道题,先用kruskal求一遍图中的最大生成树. 然后,倍增求lca,求lca的同时求出边权的最小值. #include <cstring> #include <cstdi ...
- 【NOIP2013】货车运输 最大生成树+LCA
题目描述 AA国有nn座城市,编号从 1到n,城市之间有m 条双向道路.每一条道路对车辆都有重量限制,简称限重.现在有 q 辆货车在运输货物, 司机们想知道每辆车在不超过车辆限重的情况下,最多能运多重 ...
- Luogu1967 NOIP2013 货车运输 最大生成树、倍增
传送门 题意:给出一个$N$个节点.$M$条边的图,$Q$次询问,每一次询问两个点之间的所有可行路径中经过的边的边权的最小值中的最大值.$N \leq 10000 , M \leq 50000 , Q ...
随机推荐
- 【转】C#异步的世界【上】
[转]C#异步的世界[上] 新进阶的程序员可能对async.await用得比较多,却对之前的异步了解甚少.本人就是此类,因此打算回顾学习下异步的进化史. 本文主要是回顾async异步模式之前的异步,下 ...
- 【Linux】zlib安装
zlib简介 zlib是提供数据压缩用的函式库,由Jean-loup Gailly与Mark Adler所开发,初版0.9版在1995年5月1日发表.zlib使用DEFLATE算法,最初是为libpn ...
- keras框架的MLP手写数字识别MNIST,梳理?
keras框架的MLP手写数字识别MNIST 代码: # coding: utf-8 # In[1]: import numpy as np import pandas as pd from kera ...
- Codeforces 1086 简要题解
文章目录 A题 B题 C题 D题 E题 传送门 这场比赛原地爆炸了啊!!! 只做了两道. A题 传送门 手贱没关freopenfreopenfreopen于是wawawa了一次,死活调不出错. 题意: ...
- Educational Codeforces Round 51 D. Bicolorings(dp)
https://codeforces.com/contest/1051/problem/D 题意 一个2*n的矩阵,你可以用黑白格子去填充他,求联通块数目等于k的方案数,答案%998244353. 思 ...
- C++STL deque
deque双端数组 deque<int> dq; deque<int>::iterator it; dq.push_back();//尾部插入元素 dq.push_front( ...
- mysql order by 中文 排序
mysql order by 中文 排序 1. 在MySQL中,我们经常会对一个字段进行排序查询,但进行中文排序和查找的时候,对汉字的排序和查找结果往往都是错误的. 这种情况在MySQL的很多版本中都 ...
- 学以致用十二-----YouCompeteMe巨坑
接上一篇,通过这几天的不断尝试,发现一个无法解决的问题.至于我安装成功的那台,我至今不知道是安装了哪一步导致成功的. 首先,我在.vimrc里加上了 Plugin 'Valloric/YouComp ...
- 1.2.2实现Runnable接口
使用Runnable创建线程 package com.cky.runner; /** * Created by chenkaiyang on 2017/12/2. */ public class My ...
- python在使用redis时zadd错误
最近在看<redis实战>,在写zadd时报错 Traceback (most recent call last): File "<stdin>", lin ...