[BZOJ2402]陶陶的难题II(树链剖分+线段树维护凸包+分数规划)
陶陶的难题II
时间限制:40s 空间限制:128MB
题目描述
输入格式
第一行包含一个正整数N,表示树中结点的个数。
第二行包含N个正实数,第i个数表示xi (1<=xi<=10^5)。
第三行包含N个正实数,第i个数表示yi (1<=yi<=10^5)。
第四行包含N个正实数,第i个数表示pi (1<=pi<=10^5)。
第五行包含N个正实数,第i个数表示qi (1<=qi<=10^5)。
下面有N-1行,每行包含两个正整数a,b(1<=a,b<=N),表示树中的边。
第N+5行包含一个正整数M,表示询问的个数。
最后M行,每行包含正整数a,b(1<=a,b<=N),表示一次询问。
输出格式
共M行,每行一个实数,第i行的数表示第i次询问的答案。
只要你的输出和我们的输出相差不超过0.001即为正确。
样例输入
5
3.0 1.0 2.0 5.0 4.0
5.0 2.0 4.0 3.0 1.0
1.0 3.0 2.0 4.0 5.0
3.0 4.0 2.0 1.0 4.0
1 2
1 3
2 4
2 5
4
2 3
4 5
2 4
3 5
样例输出
2.5000
1.5000
1.5000
2.5000
提示
100%的数据满足N,M≤ 30,000。
1<=Xi,Yi,Pi,Qi<=10^8
题目来源
没有写明来源
首先根据答案是个分数以及“你的输出和我们的输出相差不超过0.001即为正确”可以想到分数规划。然后化下式子可以发现i和j是可以分开考虑的,加上题目求的式子是一个分数可以自然想到凸包。但是我们怎么维护这个凸包呢?由于没有修改,我们可以通过线段树套std::vector维护,合并时由于两边横坐标都是单增的所以用类似归并排序的方法即可。最后,一般凸包问题都是用三分,但是仔细思考会发现其实二分也可以做。这样总的复杂度就是四个log的(分数规划,树链剖分,线段树,二分查找各一个)。由于复杂度并不满所以可以跑过(当然能用$O(n \log^4 n)$卡掉$O(n^2)$的题并不多见)
#include<cstdio>
#include<vector>
#include<algorithm>
#define ls (x<<1)
#define rs ((x<<1)|1)
#define lson ls,L,mid
#define rson rs,mid+1,R
#define rep(i,l,r) for (int i=l; i<=r; i++)
#define For(i,x) for (int i=h[x],k; i; i=nxt[i])
typedef long double ld;
using namespace std; const int N=;
int n,m,u,v,cnt,nd,dfn[N],d[N],top[N],sz[N],son[N],f[N],h[N],to[N<<],nxt[N<<];
ld vx[N],vy[N],vp[N],vq[N]; void add(int u,int v){ to[++cnt]=v; nxt[cnt]=h[u]; h[u]=cnt; } struct Seg{
ld px[N],py[N];
vector<int>V[N<<]; void insert(vector<int>&V,int p){
while (V.size()> && (px[V[V.size()-]]-px[p])*(py[V[V.size()-]]-py[p])-(py[V[V.size()-]]-py[p])*(px[V[V.size()-]]-px[p])<-1e-) V.pop_back();
V.push_back(p);
} void build(int x,int L,int R){
if (L==R){ V[x].push_back(L); return; }
int mid=(L+R)>>;
build(lson); build(rson);
int pl=,pr=,tl=V[ls].size(),tr=V[rs].size();
while (pl<tl || pr<tr)
if (pr==tr || (pl<tl && px[V[ls][pl]]<px[V[rs][pr]])) insert(V[x],V[ls][pl++]);
else insert(V[x],V[rs][pr++]);
} ld que(int x,int L,int R,int l,int r,ld k){
if (L==l && R==r){
int ll=,rr=V[x].size()-,mid,ans=;
while (ll<=rr){
mid=(ll+rr)>>;
if (py[V[x][mid]]-k*px[V[x][mid]]-1e->py[V[x][mid-]]-k*px[V[x][mid-]]) ans=mid,ll=mid+; else rr=mid-;
}
return py[V[x][ans]]-k*px[V[x][ans]];
}
int mid=(L+R)>>;
if (r<=mid) return que(lson,l,r,k);
else if (l>mid) return que(rson,l,r,k);
else return max(que(lson,l,mid,k),que(rson,mid+,r,k));
}
}A,B; ld solve(int x,int y,ld ans){
ld res1=-1e16,res2=-1e16;
while (top[x]!=top[y]){
if (d[top[x]]<d[top[y]]) swap(x,y);
res1=max(res1,A.que(,,n,dfn[top[x]],dfn[x],ans));
res2=max(res2,B.que(,,n,dfn[top[x]],dfn[x],ans));
x=f[top[x]];
}
if (d[x]>d[y]) swap(x,y);
res1=max(res1,A.que(,,n,dfn[x],dfn[y],ans));
res2=max(res2,B.que(,,n,dfn[x],dfn[y],ans));
return res1+res2;
} void dfs(int x,int fa){
sz[x]=; son[x]=;
For(i,x) if ((k=to[i])!=fa){
f[k]=x; d[k]=d[x]+; dfs(k,x); sz[x]+=sz[k];
if (sz[k]>sz[son[x]]) son[x]=k;
}
} void dfs1(int x,int tp){
dfn[x]=++nd; top[x]=tp; A.px[nd]=vx[x]; A.py[nd]=vy[x]; B.px[nd]=vp[x]; B.py[nd]=vq[x];
if (son[x]) dfs1(son[x],tp);
For(i,x) if ((k=to[i])!=f[x] && k!=son[x]) dfs1(k,k);
} int main(){
freopen("bzoj2402.in","r",stdin);
freopen("bzoj2402.out","w",stdout);
scanf("%d",&n);
rep(i,,n) scanf("%Lf",&vx[i]);
rep(i,,n) scanf("%Lf",&vy[i]);
rep(i,,n) scanf("%Lf",&vp[i]);
rep(i,,n) scanf("%Lf",&vq[i]);
rep(i,,n-) scanf("%d%d",&u,&v),add(u,v),add(v,u);
dfs(,); dfs1(,);
A.build(,,n); B.build(,,n);
for (scanf("%d",&m); m--; ){
scanf("%d%d",&u,&v);
ld L=,R=1e8; int c=;
while (c--){
ld mid=(L+R)/;
if (solve(u,v,mid)>1e-) L=mid; else R=mid;
}
printf("%.5Lf\n",L);
}
return ;
}
[BZOJ2402]陶陶的难题II(树链剖分+线段树维护凸包+分数规划)的更多相关文章
- 【bzoj2402】陶陶的难题II 分数规划+树链剖分+线段树+STL-vector+凸包+二分
题目描述 输入 第一行包含一个正整数N,表示树中结点的个数.第二行包含N个正实数,第i个数表示xi (1<=xi<=10^5).第三行包含N个正实数,第i个数表示yi (1<=yi& ...
- BZOJ 3672[NOI2014]购票(树链剖分+线段树维护凸包+斜率优化) + BZOJ 2402 陶陶的难题II (树链剖分+线段树维护凸包+分数规划+斜率优化)
前言 刚开始看着两道题感觉头皮发麻,后来看看题解,发现挺好理解,只是代码有点长. BZOJ 3672[NOI2014]购票 中文题面,题意略: BZOJ 3672[NOI2014]购票 设f(i)f( ...
- 【bzoj1959】[Ahoi2005]LANE 航线规划 树链剖分+线段树
题目描述 对Samuel星球的探险已经取得了非常巨大的成就,于是科学家们将目光投向了Samuel星球所在的星系——一个巨大的由千百万星球构成的Samuel星系. 星际空间站的Samuel II巨型计算 ...
- 【BZOJ-2325】道馆之战 树链剖分 + 线段树
2325: [ZJOI2011]道馆之战 Time Limit: 40 Sec Memory Limit: 256 MBSubmit: 1153 Solved: 421[Submit][Statu ...
- 【BZOJ2243】[SDOI2011]染色 树链剖分+线段树
[BZOJ2243][SDOI2011]染色 Description 给定一棵有n个节点的无根树和m个操作,操作有2类: 1.将节点a到节点b路径上所有点都染成颜色c: 2.询问节点a到节点b路径上的 ...
- BZOJ2243 (树链剖分+线段树)
Problem 染色(BZOJ2243) 题目大意 给定一颗树,每个节点上有一种颜色. 要求支持两种操作: 操作1:将a->b上所有点染成一种颜色. 操作2:询问a->b上的颜色段数量. ...
- POJ3237 (树链剖分+线段树)
Problem Tree (POJ3237) 题目大意 给定一颗树,有边权. 要求支持三种操作: 操作一:更改某条边的权值. 操作二:将某条路径上的边权取反. 操作三:询问某条路径上的最大权值. 解题 ...
- bzoj4034 (树链剖分+线段树)
Problem T2 (bzoj4034 HAOI2015) 题目大意 给定一颗树,1为根节点,要求支持三种操作. 操作 1 :把某个节点 x 的点权增加 a . 操作 2 :把某个节点 x 为根的子 ...
- HDU4897 (树链剖分+线段树)
Problem Little Devil I (HDU4897) 题目大意 给定一棵树,每条边的颜色为黑或白,起始时均为白. 支持3种操作: 操作1:将a->b的路径中的所有边的颜色翻转. 操作 ...
- Aizu 2450 Do use segment tree 树链剖分+线段树
Do use segment tree Time Limit: 1 Sec Memory Limit: 256 MB 题目连接 http://www.bnuoj.com/v3/problem_show ...
随机推荐
- 【CODEVS】1922 骑士共存问题
[算法]二分图最大匹配(最大流) [题解]按(i+j)奇偶性染色后,发现棋子跳到的地方刚好异色. 然后就是二分图了,对于每个奇点向可以跳到的地方连边,偶点不需连(可逆). 所以题目要求转换为求二分图上 ...
- Name That Number 命名那个数字
1.2.3 Name That Number 命名那个数字 Time Limit: 1 Sec Memory Limit: 64 MBSubmit: 183 Solved: 33[Submit][ ...
- css 格式中id与class共存
PHP文件中有一段:<div class="post-alt blog" id="post-alt"> CSS文件中有一段:.post-alt {X ...
- BBScan — 一个信息泄漏批量扫描脚本
github:https://github.com/lijiejie/BBScan 有些朋友手上有几十万甚至上百万个域名,比如,乌云所有厂商的子域名. 如果把这30万个域名全部扔给wvs,APPsca ...
- gmail注册时“此电话号码无法用于进行验证”
网上有几个方法,有说不要改默认地点,有说验证时直接写+86手机号,试了以后还是不行. 我的方法:换成IE浏览器,就可以验证了.
- mysql 配置数据库主从同步
参考:https://www.cnblogs.com/kevingrace/p/6256603.html http://www.51testing.com/html/00/130600-243651. ...
- 画弧线DrawArc的研究-我自己 -- 直线交接圆角
procedure TForm4.Button7Click(Sender: TObject); var pwith: Integer; //画笔的宽度 hx1, hy1: Integer; //横线第 ...
- JAVA封装消息中间件调用二(kafka消费者篇)
上一遍我简单介绍了kafka的生成者使用,调用方式比较简单,今天我给大家分享下封装kafka消费者,作为中间件,我们做的就是最大程度的解耦,使业务方接入我们依赖程度降到最低. 第一步,我们先配置一个消 ...
- golang之终端操作,文件操作
终端操作 操作终端相关的文件句柄常量os.Stdin:标准输入os.Stdout:标准输出os.Stderr:标准错误输出 关于终端操作的代码例子: package main import " ...
- C#窗口矩形区域着色
C#写的一个GUI窗口,有几百个矩形区域.每个矩形区域的颜色随时都可能改变,并且多次改变. 我放弃使用label绘制矩形,因为效果不好.拖控件的界面使用power packs中的rectanglesh ...