题意

4908 Race 0x49「数据结构进阶」练习

描述

给定一棵 N 个节点的树,每条边带有一个权值。

求一条简单路径,路径上各条边的权值和等于K,且路径包含的边的数量最少。

输入格式

第一行两个整数 N, K。

第2~N行每行三个整数x,y,z,表示一条无向边的两个端点x,y和权值z,点的编号从0开始。

输出格式

一个整数,表示最少边数量。如果不存在满足要求的路径,输出-1。

样例输入

4 3
0 1 1
1 2 2
1 3 4

样例输出

2

数据范围与约定

  • N <= 200000, K <= 1000000

来源

IOI2011

        </article>

分析

训练指南的配套代码有问题,双指针绝对是错的。学习了hzwer的做法。

开一个100W的数组t,t[i]表示权值为i的路径最少边数

找到重心分成若干子树后, 得出一棵子树的所有点到根的权值和x,到根c条边,用t[k-x]+c更新答案,全部查询完后

然后再用所有c更新t[x]

这样可以保证不出现点分治中的不合法情况

把一棵树的所有子树搞完后再遍历所有子树恢复T数组,如果用memset应该会比较慢

时间复杂度\(O(n \log^2 n)\)

代码

网上好多代码,点分治递归的时候没有重新统计以重心为根的整棵树每个节点的size,按道理找出来的重心是错的,复杂度没有保证。

但是出题人根本就不会卡这种东西,所以不用写了,反而常数小。

#include<bits/stdc++.h>
#define rg register
#define il inline
#define co const
template<class T>il T read(){
rg T data=0,w=1;rg char ch=getchar();
while(!isdigit(ch)) {if(ch=='-') w=-1;ch=getchar();}
while(isdigit(ch)) data=data*10+ch-'0',ch=getchar();
return data*w;
}
template<class T>il T read(rg T&x) {return x=read<T>();}
typedef long long ll;
using namespace std;
typedef pair<int,int> pii; co int N=2e5+1;
int n,k,ans,sum,root,max_size;
vector<pii> e[N];
int t[1000001],s[N],d[N],c[N],v[N];
void getroot(int x){
v[x]=1,s[x]=1;
int max_part=0;
for(int i=0,y;i<e[x].size();++i){
if(v[y=e[x][i].first]) continue;
getroot(y);
s[x]+=s[y],max_part=max(max_part,s[y]);
}
max_part=max(max_part,sum-s[x]);
if(max_part<max_size) root=x,max_size=max_part;
v[x]=0;
}
void cal(int x){
v[x]=1;
if(d[x]<=k) ans=min(ans,c[x]+t[k-d[x]]);
for(int i=0,y;i<e[x].size();++i){
if(v[y=e[x][i].first]) continue;
c[y]=c[x]+1,d[y]=d[x]+e[x][i].second;
cal(y);
}
v[x]=0;
}
void add(int x,int flag){
v[x]=1;
if(d[x]<=k){
if(flag) t[d[x]]=min(t[d[x]],c[x]);
else t[d[x]]=n;
}
for(int i=0,y;i<e[x].size();++i){
if(v[y=e[x][i].first]) continue;
add(y,flag);
}
v[x]=0;
}
void work(int x){
v[x]=1,t[0]=0;
for(int i=0,y;i<e[x].size();++i){
if(v[y=e[x][i].first]) continue;
c[y]=1,d[y]=e[x][i].second;
cal(y),add(y,1);
}
for(int i=0,y;i<e[x].size();++i){
if(v[y=e[x][i].first]) continue;
add(y,0);
}
for(int i=0,y;i<e[x].size();++i){
if(v[y=e[x][i].first]) continue;
sum=max_size=s[y];
getroot(y),work(root);
}
}
int main(){
// freopen(".in","r",stdin),freopen(".out","w",stdout);
read(n),read(k);
fill(t+1,t+k+1,n);
for(int i=1,u,v,w;i<n;++i){
read(u),read(v),read(w);
e[++u].push_back(pii(++v,w)),e[v].push_back(pii(u,w));
}
ans=sum=max_size=n;
getroot(1),work(root);
if(ans==n) puts("-1");
else printf("%d\n",ans);
return 0;
}

CH4908 Race的更多相关文章

  1. Promise.race

    [Promise.race] 返回最先完成的promise var p1 = new Promise(function(resolve, reject) { setTimeout(resolve, 5 ...

  2. golang中的race检测

    golang中的race检测 由于golang中的go是非常方便的,加上函数又非常容易隐藏go. 所以很多时候,当我们写出一个程序的时候,我们并不知道这个程序在并发情况下会不会出现什么问题. 所以在本 ...

  3. 【BZOJ-2599】Race 点分治

    2599: [IOI2011]Race Time Limit: 70 Sec  Memory Limit: 128 MBSubmit: 2590  Solved: 769[Submit][Status ...

  4. hdu 4123 Bob’s Race 树的直径+rmq+尺取

    Bob’s Race Time Limit: 5000/2000 MS (Java/Others)    Memory Limit: 32768/32768 K (Java/Others) Probl ...

  5. Codeforces Round #131 (Div. 2) E. Relay Race dp

    题目链接: http://codeforces.com/problemset/problem/214/E Relay Race time limit per test4 secondsmemory l ...

  6. 【多线程同步案例】Race Condition引起的性能问题

    Race Condition(也叫做资源竞争),是多线程编程中比较头疼的问题.特别是Java多线程模型当中,经常会因为多个线程同时访问相同的共享数据,而造成数据的不一致性.为了解决这个问题,通常来说需 ...

  7. Codeforces Round #328 (Div. 2) C. The Big Race 数学.lcm

    C. The Big Race Time Limit: 20 Sec Memory Limit: 256 MB 题目连接 http://codeforces.com/contest/592/probl ...

  8. HDU 4123 Bob’s Race 树的直径 RMQ

    Bob’s Race Time Limit: 20 Sec Memory Limit: 256 MB 题目连接 http://acm.hdu.edu.cn/showproblem.php?pid=41 ...

  9. [LOJ 1038] Race to 1 Again

    C - Race to 1 Again Time Limit:2000MS     Memory Limit:32768KB     64bit IO Format:%lld & %llu D ...

随机推荐

  1. TTL反相器的外部特性

    TTL反相器的外部特性 电压传输特性 输入端噪声容限特性 静态输入特性: 静态输出特性: 动态特性: 传输延迟时间:是由晶体管的延迟时间,电阻以及寄生电容元素引起的.包括俩部分:输入由低电平跳为高电平 ...

  2. 深入理解Connector

    上一篇主要是从各个容器的生命周期的角度讲了一下整个tomcat的运行流程,说明了各个容器之间的调用关系.但并没有太过详细的说明每一个组件并区分他们. 下面从功能的角度上详细的分析一下connector ...

  3. java的类class 和对象object

    java 语言的源代码是以类为单位存放在文件中,已public修饰的类名须和存放这个类的源文件名一样.而 一个源文件中只能有一个public的类,类名的首字母通常为大写. 使用public修饰的类可以 ...

  4. IDEA中自动生成serialVersionUID

    File  >>  Setting  >>  Inspections  >>  serializable 勾选上后,光标放在实现Serializable接口的类名上 ...

  5. JAVA 中的MessageDigest类和Mac类的使用

    MessageDigest 消息摘要 例子: MD5加密: try{ MessageDigest md5 = MessageDigest.getInstance("MD5"); m ...

  6. GD2模块-图像处理

    GD2模块-图像处理 1.图像处理模块的主要功能: a) 验证码 b) 加盖水印 c) 缩略图 d) 帖子图片签名 e) 在线LOGO制作 2确认PHP是否支持图像处理 检测PHPINFO文件中是否存 ...

  7. Unity3d代码及效率优化总结

    1.PC平台的话保持场景中显示的顶点数少于200K~3M,移动设备的话少于10W,一切取决于你的目标GPU与CPU. 2.如果你用U3D自带的SHADER,在表现不差的情况下选择Mobile或Unli ...

  8. Java语法基础学习DaySix

    一.JavaBean——可重用组件 1.JavaBean是指符合以下标准的Java类: (1)类是公共的 (2)有一个无参的公共的构造器 (3)有属性,且有对应的get.set方法 2.好处 用户可以 ...

  9. Python 数据结构--查找

    1 顺序查找O(n) def sequential_search(a_list, item): pos = 0 found = False while pos < len(a_list) and ...

  10. JetBrains Goland 2018.2.1最新版破解

    JetBrains GoLand 2018.2.1破解版是JetBrains的新商业IDE,旨在为Go开发提供符合人体工程学的环境.新的IDE通过特定于Go语言的编码辅助和工具集成扩展了Intelli ...