NOIP 模拟 $17\; \rm weight$
题解 \(by\;zj\varphi\)
一道树剖的题
先对于原图求出一棵最小生成树,求出来的这棵树中的边定为树边,其它边叫非树边
那么对于一条非树边,它要成为最小生成树上的边,权值只能为连接它两个端点的树上路径最小值减 \(1\),因为,他和这些树边构成了一个环,最小生成树要在环上去掉最大的一条边,所以只要这条边比树上路径最大值小即可
而对于树边,发现它要比它参与构成的所有环中的非树边权值要小,因为非树边在其环中权职一定为最大的,若树边比它大,就不会在最小生成树中了
对于这两种操作,我们都可以用树剖来做,就是一个板子
注意:维护最小值的线段树维护的是边下放到它终点后的点权,记得在 \(lca\) 处不要多修改
Code
#include<bits/stdc++.h>
#define ri register signed
#define p(i) ++i
using namespace std;
namespace IO{
char buf[1<<21],*p1=buf,*p2=buf;
#define gc() p1==p2&&(p2=(p1=buf)+fread(buf,1,1<<21,stdin),p1==p2)?EOF:*p1++
template<typename T>inline void read(T &x) {
ri f=1;x=0;register char ch=gc();
while(ch<'0'||ch>'9') {if (ch=='-') f=0;ch=gc();}
while(ch>='0'&&ch<='9') {x=(x<<1)+(x<<3)+(ch^48);ch=gc();}
x=f?x:-x;
}
}
using IO::read;
namespace nanfeng{
#define FI FILE *IN
#define FO FILE *OUT
template<typename T>inline T cmax(T x,T y) {return x>y?x:y;}
template<typename T>inline T cmin(T x,T y) {return x>y?y:x;}
static const int N=7e4+7,M=1e5+7,INF=2e9+7;
int first[N],fa[N],W[N],hs[N],siz[N],dfn[N],top[N],bc[N],ans[M],dep[N],vis[M],tot,t=1,n,m,opt,num;
struct edge{int v,nxt,w;}e[N<<1];
struct edgei{int u,v,w,id;}G[M],AS[M];
inline int cmp(edgei e1,edgei e2) {return e1.w<e2.w;}
inline void add(int u,int v,int w) {
e[t].v=v,e[t].w=w,e[t].nxt=first[u],first[u]=t++;
e[t].v=u,e[t].w=w,e[t].nxt=first[v],first[v]=t++;
}
int find(int x) {return fa[x]==x?x:fa[x]=find(fa[x]);}
inline void kruskal() {
sort(G+1,G+m+1,cmp);
for (ri i(1);i<=m;p(i)) {
int u=G[i].u,v=G[i].v,w=G[i].w;
if (find(u)==find(v)) continue;
vis[G[i].id]=1;
fa[find(u)]=find(v);
add(u,v,w);
}
}
struct Segmx{
#define ls(x) (x<<1)
#define rs(x) (x<<1|1)
#define up(x) T[x].mx=cmax(T[ls(x)].mx,T[rs(x)].mx);
struct segmenttree{int mx;}T[N<<2];
void build(int x,int l,int r) {
if (l==r) return (void)(T[x].mx=W[bc[l]]);
int mid(l+r>>1);
build(ls(x),l,mid);
build(rs(x),mid+1,r);
up(x);
}
int querymx(int x,int l,int r,int lt,int rt) {
if (l<=lt&&rt<=r) return T[x].mx;
int mid(lt+rt>>1),res(0);
if (l<=mid) res=cmax(res,querymx(ls(x),l,r,lt,mid));
if (r>mid) res=cmax(res,querymx(rs(x),l,r,mid+1,rt));
return res;
}
}Tx;
struct Segmn{
#define ls(x) (x<<1)
#define rs(x) (x<<1|1)
struct segmenttree{int mn,lz;segmenttree(){mn=INF;}}T[N<<2];
inline void down(int x) {
if (!T[x].lz) return;
T[ls(x)].mn=cmin(T[x].mn,T[ls(x)].mn);
T[rs(x)].mn=cmin(T[x].mn,T[rs(x)].mn);
T[ls(x)].lz=T[rs(x)].lz=1;
T[x].lz=0;
}
void update(int x,int k,int l,int r,int lt,int rt) {
if (l<=lt&&rt<=r) {
T[x].mn=cmin(T[x].mn,k);
return (void)(T[x].lz=1);
}
down(x);
int mid(lt+rt>>1);
if (l<=mid) update(ls(x),k,l,r,lt,mid);
if (r>mid) update(rs(x),k,l,r,mid+1,rt);
}
int querymn(int x,int k,int l,int r) {
if (l==r) return T[x].mn;
int mid(l+r>>1);
down(x);
if (k<=mid) return querymn(ls(x),k,l,mid);
else return querymn(rs(x),k,mid+1,r);
}
}Tn;
void dfs1(int x,int f) {
fa[x]=f;
siz[x]=1;
for (ri i(first[x]),v;i;i=e[i].nxt) {
if ((v=e[i].v)==f) continue;
W[v]=e[i].w,dep[v]=dep[x]+1;
dfs1(v,x);
siz[x]+=siz[v];
if (siz[v]>siz[hs[x]]) hs[x]=v;
}
}
void dfs2(int x,int tp) {
bc[dfn[x]=p(tot)]=x;
top[x]=tp;
if (hs[x]) dfs2(hs[x],tp);
for (ri i(first[x]),v;i;i=e[i].nxt) {
if ((v=e[i].v)==fa[x]||v==hs[x]) continue;
dfs2(v,v);
}
}
inline int solve(int u,int v,int w) {
ri res(0);
while(top[u]!=top[v]) {
if (dep[top[u]]<dep[top[v]]) swap(u,v);
res=cmax(res,Tx.querymx(1,dfn[top[u]],dfn[u],1,n));
Tn.update(1,w,dfn[top[u]],dfn[u],1,n);
u=fa[top[u]];
}
if (dep[u]<dep[v]) swap(u,v);
if (dfn[u]>dfn[v]) {
res=cmax(res,Tx.querymx(1,dfn[v]+1,dfn[u],1,n));
Tn.update(1,w,dfn[v]+1,dfn[u],1,n);
}
return res-1;
}
inline int main() {
// FI=freopen("nanfeng.in","r",stdin);
// FO=freopen("nanfeng.out","w",stdout);
read(n),read(m),read(opt);
for (ri i(1);i<=n;p(i)) fa[i]=i;
for (ri i(1);i<=m;p(i))
read(G[i].u),read(G[i].v),read(G[i].w),G[i].id=i,AS[i]=G[i];
kruskal();
dfs1(1,0),dfs2(1,1);
Tx.build(1,1,n);
for (ri i(1);i<=m;p(i)) {
int u=AS[i].u,v=AS[i].v,w=AS[i].w;
if (!dfn[u]) continue;
if (vis[i]) continue;
ans[i]=solve(u,v,w);
}
for (ri i(1);i<=m;p(i)) {
int u=AS[i].u,v=AS[i].v,w=AS[i].w,lw=dep[u]<dep[v]?v:u;
if (!dfn[u]) continue;
if (vis[i]) ans[i]=Tn.querymn(1,dfn[lw],1,n)-1;
}
for (ri i(1);i<=m;p(i)) printf("%d ",ans[i]==INF-1?-1:ans[i]);
return 0;
}
}
int main() {return nanfeng::main();}
NOIP 模拟 $17\; \rm weight$的更多相关文章
- NOIP 模拟 $17\; \rm 时间机器$
题解 \(by\;zj\varphi\) 一道贪心的题目 我们先将节点和电阻按左边界排序,相同的按右边界排序 对于每一个节点,我们发现,选取左边界小于等于它的电阻中右边界大于它且最接近的它的一定是最优 ...
- NOIP 模拟 $17\; \rm 世界线$
题解 \(by\;zj\varphi\) 此题经简单观察可发现,一个点的贡献就是这个点所能到的点减去它的出度 那么我们就可以暴力搜索,但是显然会超时,所以我们可以使用一个黑科技 \(\rm bitse ...
- NOIP模拟17.9.22
NOIP模拟17.9.22 前进![问题描述]数轴的原点上有一只青蛙.青蛙要跳到数轴上≥
- NOIP模拟17.9.21
NOIP模拟17.9.21 3 58 145 201 161.5 样例输出21.6 数据规模及约定对于40% 的数据,N <= 20对于60% 的数据,N <= 1000对于100% 的数 ...
- NOIP模拟 17.8.20
NOIP模拟17.8.20 A.阶乘[题目描述]亲爱的xyx同学正在研究数学与阶乘的关系,但是他喜欢颓废,于是他就制作了一个和阶乘有关系的数学游戏:给出两个整数 n,m,令 t = !n,每轮游戏的流 ...
- NOIP模拟 17.8.18
NOIP模拟17.8.18 A.小菜一碟的背包[题目描述]Blice和阿强巴是好朋友但萌萌哒Blice不擅长数学,所以阿强巴给了她一些奶牛做练习阿强巴有 n头奶牛,每头奶牛每天可以产一定量的奶,同时也 ...
- NOIP模拟 17.8.15
NOIP模拟17.8.15 A 债务文件名 输入文件 输出文件 时间限制 空间限制debt.pas/c/cpp debt.in debt.out 1s 128MB[题目描述]小 G 有一群好朋友,他们 ...
- NOIP模拟 17.8.17
NOIP模拟17.8.17 A 小 G 的字符串文件名 输入文件 输出文件 时间限制 空间限制str.pas/c/cpp str.in str.out 1s 128MB[题目描述]有一天,小 L 给小 ...
- NOIP模拟 17.8.16
NOIP模拟17.8.16 A 债务文件名 输入文件 输出文件 时间限制 空间限制debt.pas/c/cpp debt.in debt.out 1s 128MB[题目描述]小 G 有一群好朋友,他们 ...
随机推荐
- 关于vector.size()的一些常见错误总结
1. 问题引入 通过查看[https://www.cplusplus.com/reference/vector/vector/] 的vector.size()说明,即 member type defi ...
- 『心善渊』Selenium3.0基础 — 27、unittest跳过测试的使用
目录 1.什么是跳过测试 2.常用的跳过测试方法和装饰器 3.跳过测试示例 4.TestCase.skipTest()方法 1.什么是跳过测试 当测试用例写完后,有些模块有改动时候,会影响到部分用例的 ...
- Python3.7 lxml引入etree
用xml代替lxml,Python3.7中已经没有etree这个模块了 import xml.etree.ElementTree as etree from lxml import etree 这种方 ...
- java集合(3)-Java8新增的Stream操作集合
Java8新增了Stream,IntStream,LongStream,DoubleStream等流式API,这些API代表多个支持串行和并行聚集操作的元素.上面的4个接口中,Stream是一个通用的 ...
- Jupyter使用快捷键
命令行模式(按 Esc 生效) F: 查找并且替换 Ctrl-Shift-F: 打开命令配置 Ctrl-Shift-P: 打开命令配置 Enter: 进入编辑模式 P: 打开命令配置 Shift-En ...
- 【动态规划】石子合并 luogu-1880
分析 简单的区间DP AC代码 #include <bits/stdc++.h> using namespace std; #define ms(a,b) memset(a,b,sizeo ...
- vue点击复制功能
复制功能,选中复制或者点击复制(不使用插件的情况下) 1.选中复制 这个比点击复制简单点 <template> <div> <el-button type=& ...
- MERCY靶机
仅供个人娱乐 靶机信息 下载地址:https://drive.google.com/uc?id=1YzsW1lCKjo_WEr6Pk511DXQBFyMMR14y&export=downloa ...
- 你的ES数据备份了吗?
前言: 无论使用哪种存储软件,定期的备份数据都是重中之重,在使用ElasticSearch的时候,随着数据日益积累,存放es数据的磁盘空间也捉襟见肘, 此时对于业务功能使用不到的索引数据,又不能直接删 ...
- 痞子衡嵌入式:在IAR开发环境下将关键函数重定向到RAM中执行的三种方法
大家好,我是痞子衡,是正经搞技术的痞子.今天痞子衡给大家分享的是在IAR开发环境下将关键函数重定向到RAM中执行的三种方法. 嵌入式项目里应用程序代码正常是放在 Flash 中执行的,但有时候也需要将 ...