UOJ276 [清华集训2016] 汽水 【二分答案】【点分治】【树状数组】
题目分析:
这种乱七八糟的题目一看就是点分治,答案有单调性,所以还可以二分答案。
我们每次二分的时候考虑答案会不会大于等于某个值,注意到系数$k$是无意义的,因为我们可以通过转化使得$k=0$。
合并的过程相当于很多个向量,加起来后看斜率。
注意单个向量也要判定。
由于有了二分的答案$Ans$。判定变得简单多了,推一下。
令$(A,C)$是从一个点到重心的一个向量,$(B,D)$是从另一个点到重心的向量。其中$A$和$B$是重心到该点的路径权值和,$C$和$D$是经过的边数。
$-k \leq \frac{A+C}{B+D} \leq k \Rightarrow -k(B+D) \leq A+C \leq k(B+D)$.
进一步的$A+kB \geq -C-kD$且$A-kB \leq kD-C$。虽然有四元,但是顺序相互关联,所以实际只有两元,排序后树状数组就可以解决啦。
#include<bits/stdc++.h>
using namespace std; typedef long long ll; const int maxn = ; ll k,md; int flag = ,num,n,rnum;
vector <pair<int,ll> > g[maxn];
int arr[maxn],sz[maxn],imp[maxn],cnt[maxn];
struct node{ll A;int B,pla;}op[maxn]; int cmp(node X,node Y){
return -X.A-md*X.B < -Y.A-md*Y.B;
} struct Fenwick{
int C[maxn];
void Add(int now){
while(now <= rnum){C[now] ++; now += (now&-now);}
}
int query(int now){
int ans = ;
while(now){ans += C[now]; now -= (now&-now);}
return ans;
}
}T1; void read(){
scanf("%d%lld",&n,&k);
for(int i=;i<n;i++){
int x,y;long long v; scanf("%d%d%lld",&x,&y,&v); v -= k;
g[x].push_back(make_pair(y,v)); g[y].push_back(make_pair(x,v));
}
} void dfs1(int now,int fa,int dp){
sz[now] = ;imp[now] = ;
for(auto it : g[now]){
if((arr[it.first] && arr[it.first] < dp) || fa == it.first) continue;
dfs1(it.first,now,dp); sz[now] += sz[it.first];
}
} int dfs2(int now,int fa,int dp,int ssz){
int ans = ;
for(auto it : g[now]){
if((arr[it.first] && arr[it.first] < dp) || fa == it.first) continue;
int data = dfs2(it.first,now,dp,ssz);
if(ans== || imp[ans] > imp[data])ans = data;
imp[now] = max(sz[it.first],imp[now]);
}
imp[now] = max(imp[now],ssz-sz[now]);
if(ans== || imp[ans] > imp[now]) ans = now;
return ans;
} void dfs3(int now,int fa,int dp,int A,int B){
for(auto it : g[now]){
if(it.first == fa || (arr[it.first] && arr[it.first] < dp)) continue;
op[++num] = (node){A+it.second,B+,it.first};
dfs3(it.first,now,dp,A+it.second,B+);
}
} long long lisan[maxn];
void solve(int dr){
rnum = num;
for(int i=;i<=num;i++){lisan[i] = -op[i].A+md*op[i].B;}
sort(lisan+,lisan+num+);rnum = unique(lisan+,lisan+num+)-lisan-;
for(int i=;i<=rnum;i++) T1.C[i]=;
for(int i=num,j=;i>=;i--){
while(j <= num && (-op[j].A-md*op[j].B < op[i].A+md*op[i].B)){
T1.Add(lower_bound(lisan+,lisan+rnum+,-op[j].A+md*op[j].B)-lisan);
j++;
}
int ans=j--T1.query(upper_bound(lisan+,lisan+rnum+,op[i].A-md*op[i].B)-lisan-);
if(op[i].A-md*op[i].B < -op[i].A+md*op[i].B && j > i)ans--;
if(dr == ){
if(ans - cnt[op[i].pla]){flag = ;return;}
}else{cnt[op[i].pla] = ans;}
}
} void divide(int now,int dp,int lst,long long AA){
dfs1(now,,dp); int heavy = dfs2(now,,dp,sz[now]);arr[heavy] = dp;
for(auto it : g[heavy]){
if(arr[it.first]&&arr[it.first]<dp) continue;
divide(it.first,dp+,heavy,it.second);
if(flag == ) return;
}
num = ; dfs3(heavy,,dp,,); sort(op+,op+num+,cmp);
solve(); num = ;
if(lst) {
num = ;dfs3(now,,dp,AA,);
sort(op+,op+num+,cmp);
for(int i=;i<=num;i++) cnt[op[i].pla] = ;
solve();
}
} void work(){
long long l = ,r = 1e13;
while(l < r){
md = (l+r)/; flag = ;
memset(arr,,sizeof(arr));
divide(,,,);
if(flag) r = md; else l = md+;
}
printf("%lld",l-);
} int main(){
read();
work();
return ;
}
UOJ276 [清华集训2016] 汽水 【二分答案】【点分治】【树状数组】的更多相关文章
- 并不对劲的uoj276. [清华集训2016]汽水
想要很对劲的讲解,请点击这里 题目大意 有一棵\(n\)(\(n\leq 50000\))个节点的树,有边权 求一条路径使该路径的边权平均值最接近给出的一个数\(k\) 输出边权平均值下取整的整数部分 ...
- [UOJ#276][清华集训2016]汽水[分数规划+点分治]
题意 给定一棵 \(n\) 个点的树,给定 \(k\) ,求 \(|\frac{\sum w(路径长度)}{t(路径边数)}-k|\)的最小值. \(n\leq 5\times 10^5,k\leq ...
- BZOJ5321 JXOI2017加法(二分答案+贪心+堆+树状数组)
二分答案后得到每个位置需要被加的次数.考虑贪心.从左到右考虑每个位置,将以该位置为左端点的区间按右端点从大到小加进堆.看该位置还需要被加多少次,如果不需要加了就不管,否则取堆顶区间将其选择,BIT实现 ...
- [bzoj4009] [HNOI2015]接水果 整体二分+扫描线+dfs序+树状数组
Description 风见幽香非常喜欢玩一个叫做 osu!的游戏,其中她最喜欢玩的模式就是接水果. 由于她已经DT FC 了The big black, 她觉得这个游戏太简单了,于是发明了一个更 加 ...
- 洛谷P1527 [国家集训队] 矩阵乘法 [整体二分,二维树状数组]
题目传送门 矩阵乘法 题目描述 给你一个N*N的矩阵,不用算矩阵乘法,但是每次询问一个子矩形的第K小数. 输入输出格式 输入格式: 第一行两个数N,Q,表示矩阵大小和询问组数: 接下来N行N列一共N* ...
- csp-s模拟测试56(10.2)Merchant「二分」·Equation「树状数组」
又死了......T1 Merchant 因为每个集合都可以写成一次函数的形式,所以假设是单调升的函数,那么随着t越大就越佳 而单调减的函数,随着t的增大结果越小,所以不是单调的??? 但是我们的单调 ...
- BZOJ.4738.[清华集训2016]汽水(点分治 分数规划)
BZOJ UOJ 记\(val_i\)是每条边的边权,\(s\)是边权和,\(t\)是经过边数,\(k\)是给定的\(k\). 在点分治的时候二分答案\(x\),设\(|\frac st-k|=x\) ...
- UOJ#276. 【清华集训2016】汽水 二分答案 点分治
原文链接https://www.cnblogs.com/zhouzhendong/p/UOJ276.html 题解 首先,读入的时候就将所有的 $w_i$ 减掉 $k$ . 于是我们要求的就是平均值最 ...
- HDU 5877 2016大连网络赛 Weak Pair(树状数组,线段树,动态开点,启发式合并,可持久化线段树)
Weak Pair Time Limit: 4000/2000 MS (Java/Others) Memory Limit: 262144/262144 K (Java/Others) Tota ...
随机推荐
- 在vue2.0中引用element-ui组件库
element-ui是由饿了么团队开发的一套基于 Vue 2.0 的桌面端组件库. 官网:http://element.eleme.io/ 安装 npm i element-ui -S 引用完整的el ...
- ASP.NET MVC5+EF6+EasyUI 仓库管理系统
简介 此系统是在框架基础上叠加的仓库系统功能,同样是开源的,可以用于简单的仓库管理,您也可以在此基础上继续开发与完善! 仓库的框架功能会随着框架的升级而同样升级,仓库管理售价1666,已经购买过框架的 ...
- CSS 伪类 (Pseudo-classes)实例
CSS 伪类 (Pseudo-classes)实例CSS 伪类用于向某些选择器添加特殊的效果在支持 CSS 的浏览器中,链接的不同状态都可以不同的方式显示,这些状态包括:活动状态,已被访问状态,未被访 ...
- Ubuntu: 软件库(software repositories)
Linux 生态下,几乎每个发行版都有自己的软件库(software repositories),Ubuntu 当然也不例外.Ubuntu 提供了四个不同的软件库,分别是 main.restricte ...
- 《React Native 精解与实战》书籍连载「Node.js 简介与 React Native 开发环境配置」
此文是我的出版书籍<React Native 精解与实战>连载分享,此书由机械工业出版社出版,书中详解了 React Native 框架底层原理.React Native 组件布局.组件与 ...
- logstash grok内置规则
logstash grok 内置正则 https://github.com/elastic/logstash/blob/v1.4.2/patterns/grok-patterns USERNAME [ ...
- logstash采集与清洗数据到elasticsearch案例实战
原文地址:https://www.2cto.com/kf/201610/560348.html Logstash的使用 logstash支持把配置写入文件 xxx.conf,然后通过读取配置文件来采集 ...
- Python监控服务器利器--psutil
Python监控服务器利器--psutil 服务器的监控通过安装一些常用的监控软件之外,有时也需要运行一些shell或Python脚本:shell下可以使用系统自带的ps/free/top/df等sh ...
- PAT L2-016 愿天下有情人都是失散多年的兄妹
https://pintia.cn/problem-sets/994805046380707840/problems/994805061769609216 呵呵.大家都知道五服以内不得通婚,即两个人最 ...
- Windows10下安装VMware虚拟机并搭建CentOS系统环境
转载: http://blog.51cto.com/10085711/2069270 操作系统 Windows 10专业版(64位) VMware虚拟机 产品:VMware® Workstation ...