【NOIP模拟】从我背后出现
Description
给定n个点m条边的无向连通图,对于每条边求出强制选这条边后的最⼩⽣成树⼤⼩。
\(n\leq 10^5,m\leq 2*10^5\)
Input Format
第 1 行包含两个整数 n,m,表示点数和路径数。第 2~m+1 行每行三个整数 ui,vi,wi,表示有一条可以修的路连接 ui 和 vi ,费用为 wi ,wi ≤ 1e9。
Output Format
输出 m 行,每行一个整数,表示选择第i条路径的前提下的最小花费。
路径按输入的顺序编号为 1~m 。
Solution
原图构造最小生成树,对于一条边,如果他是最小生成树上的边,那么答案不变,
那如果不在上面,例如一条边u to v,
那么肯定是在原先树上u到v的路径上删去一条边,然后把这条边添加上去,可以发现删去的是最大的边
然后倍增预处理一下树上最大边即可
Code
#include <cstdio>
#include <algorithm>
#include <cmath>
#define N 100010
#define ll long long
using namespace std;
struct info{
int fr,to,w,nex,id;
friend bool operator < (info a,info b){
return a.w<b.w;
}
}ke[N*2],e[N*4];
int n,m,_log,tot,head[N*4],dep[N],f[N][20];
ll mx[N][20],Ans,Que[N*2];
bool used[N*2];
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*10+ch-'0';ch = getchar();}
return x*f;
}
inline void Link(int u,int v,int w){
e[++tot].to=v;
e[tot].w=w;
e[tot].nex=head[u];
head[u]=tot;
}
void dfs(int u,int fa){
for(int i=1;i<=_log;++i){
f[u][i]=f[f[u][i-1]][i-1];
mx[u][i]=max(mx[u][i-1],mx[f[u][i-1]][i-1]);
}
for(int i=head[u];i;i=e[i].nex){
int v=e[i].to;
if(v==fa) continue;
dep[v]=dep[u]+1;
mx[v][0]=e[i].w;
f[v][0]=u;
dfs(v,u);
}
}
int fa[N];
int Find(int u){
return (fa[u]==u)?u:fa[u]=Find(fa[u]);
}
inline void kruskal(){
for(int i=1;i<=n;++i) fa[i]=i;
sort(ke+1,ke+m+1);
int cnt=0;
for(int i=1;i<=m;++i){
int px=Find(ke[i].fr),py=Find(ke[i].to);
if(px!=py){
fa[px]=py;
used[i]=1;
Link(ke[i].fr,ke[i].to,ke[i].w);
Link(ke[i].to,ke[i].fr,ke[i].w);
cnt++;
Ans+=ke[i].w;
}
if(cnt==n-1) break;
}
}
inline ll LCA(int u,int v){
ll r=0;
if(dep[u]>dep[v]) swap(u,v);
int d=dep[v]-dep[u];
for(int i=0;i<=_log;++i)if(d&(1<<i)) r=max(r,mx[v][i]),v=f[v][i];
if(u==v) return r;
for(int i=_log;i>=0;--i)
if(f[u][i]!=f[v][i]) r=max(r,max(mx[u][i],mx[v][i])),u=f[u][i],v=f[v][i];
return max(r,max(mx[v][0],mx[u][0]));
}
int main(){
n=read(),m=read();
_log=log(n)/log(2);
for(int i=1;i<=m;++i)ke[i].fr=read(),ke[i].to=read(),ke[i].w=read(),ke[i].id=i;
kruskal();
dfs(1,0);
for(int i=1;i<=m;++i)
if(used[i]) Que[ke[i].id]=Ans;
else Que[ke[i].id]=Ans-LCA(ke[i].fr,ke[i].to)+ke[i].w;
for(int i=1;i<=m;++i)
printf("%lld\n",Que[i]);
return 0;
}
【NOIP模拟】从我背后出现的更多相关文章
- NOIP模拟赛20161022
NOIP模拟赛2016-10-22 题目名 东风谷早苗 西行寺幽幽子 琪露诺 上白泽慧音 源文件 robot.cpp/c/pas spring.cpp/c/pas iceroad.cpp/c/pas ...
- contesthunter暑假NOIP模拟赛第一场题解
contesthunter暑假NOIP模拟赛#1题解: 第一题:杯具大派送 水题.枚举A,B的公约数即可. #include <algorithm> #include <cmath& ...
- NOIP模拟赛 by hzwer
2015年10月04日NOIP模拟赛 by hzwer (这是小奇=> 小奇挖矿2(mining) [题目背景] 小奇飞船的钻头开启了无限耐久+精准采集模式!这次它要将原矿运到泛光之源的矿 ...
- 大家AK杯 灰天飞雁NOIP模拟赛题解/数据/标程
数据 http://files.cnblogs.com/htfy/data.zip 简要题解 桌球碰撞 纯模拟,注意一开始就在袋口和v=0的情况.v和坐标可以是小数.为保险起见最好用extended/ ...
- 队爷的讲学计划 CH Round #59 - OrzCC杯NOIP模拟赛day1
题目:http://ch.ezoj.tk/contest/CH%20Round%20%2359%20-%20OrzCC杯NOIP模拟赛day1/队爷的讲学计划 题解:刚开始理解题意理解了好半天,然后发 ...
- 队爷的Au Plan CH Round #59 - OrzCC杯NOIP模拟赛day1
题目:http://ch.ezoj.tk/contest/CH%20Round%20%2359%20-%20OrzCC杯NOIP模拟赛day1/队爷的Au%20Plan 题解:看了题之后觉得肯定是DP ...
- 队爷的新书 CH Round #59 - OrzCC杯NOIP模拟赛day1
题目:http://ch.ezoj.tk/contest/CH%20Round%20%2359%20-%20OrzCC杯NOIP模拟赛day1/队爷的新书 题解:看到这题就想到了 poetize 的封 ...
- CH Round #58 - OrzCC杯noip模拟赛day2
A:颜色问题 题目:http://ch.ezoj.tk/contest/CH%20Round%20%2358%20-%20OrzCC杯noip模拟赛day2/颜色问题 题解:算一下每个仆人到它的目的地 ...
- CH Round #52 - Thinking Bear #1 (NOIP模拟赛)
A.拆地毯 题目:http://www.contesthunter.org/contest/CH%20Round%20%2352%20-%20Thinking%20Bear%20%231%20(NOI ...
- CH Round #49 - Streaming #4 (NOIP模拟赛Day2)
A.二叉树的的根 题目:http://www.contesthunter.org/contest/CH%20Round%20%2349%20-%20Streaming%20%234%20(NOIP 模 ...
随机推荐
- 1)C语言简介(C自考学习)
C语言历史由来 世界上第一个高级语言是"ALFOL",而C的前身是ALGOL语言.1970年美国贝尔实验室的肯·汤普逊对BCPL(基本复合程序设计语言)进行了进一步的简化,突出了硬 ...
- 《Metasploit魔鬼训练营》第一章习题
书本p391.The time line is as follows: * May 7, 2007: Initial defect disclosure to the security@samba.o ...
- Swarm 如何实现 Failover?- 每天5分钟玩转 Docker 容器技术(98)
故障是在所难免的,容器可能崩溃,Docker Host 可能宕机,不过幸运的是,Swarm 已经内置了 failover 策略. 创建 service 的时候,我们没有告诉 swarm 发生故障时该如 ...
- 延迟执行之 Invoke 函数
Invoke 函数需要继承 MonoBehaviour 类后才能使用. Invoke(string str,float a):a 秒后执行名为 str 函数(只会调用一次). Invoke(strin ...
- 万能动态库调用工具IDMA(InvokeDllMethodsAdvance)
万能动态库调用工具IDMA 开发者:马兆瑞 QQ/微信:624762543 百度云下载链接:https://pan.baidu.com/s/1skW5W4H CSDN下载链接:http://d ...
- c语言贪吃蛇详解3.让蛇动起来
c语言贪吃蛇详解3.让蛇动起来 前几天的实验室培训课后作业我布置了贪吃蛇,今天有时间就来写一下题解.我将分几步来教大家写一个贪吃蛇小游戏.由于大家c语言未学完,这个教程只涉及数组和函数等知识点. 上次 ...
- java agent
cmd使用java -help可以看到关于agent参数: -agentlib:<libname>[=<选项>] 加载本机代理库 <libname>, 例如 -ag ...
- Mybatis分页插件PageHelper简单使用
一个好的讲解mybatis的博客地址http://www.jianshu.com/nb/5226994 引言 对于使用Mybatis时,最头痛的就是写分页,需要先写一个查询count的select语句 ...
- memcache调整value大小限制
> *事件背景: 当Redis有问题时按预案就会切换到本机memcache,但是我们首页 key:value现 在是1.5M同时memcache item限制是1M,导致首页写入memcache ...
- 实战-Mysql5.6.36脚本编译安装及初始化
概述 本文为centos7.3自动化编译安装mysql5.3.6的脚本及后续初始化操作,话不多少,直接上脚本. 安装脚本install.py如下: #coding=utf-8 #!/usr/bin/p ...