bzoj1038 [ZJOI2008]瞭望塔

Description

　　致力于建设全国示范和谐小村庄的H村村长dadzhi，决定在村中建立一个瞭望塔，以此加强村中的治安。我们

将H村抽象为一维的轮廓。如下图所示 我们可以用一条山的上方轮廓折线(x1, y1), (x2, y2), …. (xn, yn)来描

述H村的形状，这里x1 < x2 < …< xn。瞭望塔可以建造在[x1, xn]间的任意位置, 但必须满足从瞭望塔的顶端可

以看到H村的任意位置。可见在不同的位置建造瞭望塔，所需要建造的高度是不同的。为了节省开支，dadzhi村长

希望建造的塔高度尽可能小。请你写一个程序，帮助dadzhi村长计算塔的最小高度。

Input

　　第一行包含一个整数n，表示轮廓折线的节点数目。接下来第一行n个整数, 为x1 ~ xn. 第三行n个整数，为y1

 ~ yn。

Output

　　仅包含一个实数，为塔的最小高度，精确到小数点后三位。

Sample Input

【输入样例一】

6

1 2 4 5 6 7

1 2 2 4 2 1

【输入样例二】

4

10 20 49 59

0 10 10 0

Sample Output

【输出样例一】

1.000

【输出样例二】

14.500

HINT

 N ≤ 300，输入坐标绝对值不超过106，注意考虑实数误差带来的问题。

正解：半平面交。

题目给了你一个半平面交，还要求一个半平面交。求出两个半平面交以后，直接抽象到坐标系上，对于两个分段函数，它们的距离最小值所在的位置只可能出现在两个函数的任意一个端点上。所以枚举每个端点，计算垂直距离，求出最小值就行了。

//It is made by wfj_2048~
#include <algorithm>
#include <iostream>
#include <cstring>
#include <cstdlib>
#include <cstdio>
#include <vector>
#include <cmath>
#include <queue>
#include <stack>
#include <map>
#include <set>
#define ls (x<<1)
#define rs (x<<1|1)
#define inf (1<<30)
#define il inline
#define RG register
#define ll long long

using namespace std;

struct edge{ int nt,to; }g[100010];

int sum[500010],mx[500010],head[100010],w[100010],fa[100010],son[100010],top[100010],size[100010],dep[100010],tid[100010],dfn[100010],n,q,num,dex;

il int gi(){
RG int x=0,q=0; RG char ch=getchar(); while ((ch<'0' || ch>'9') && ch!='-') ch=getchar();
if (ch=='-') q=1,ch=getchar(); while (ch>='0' && ch<='9') x=x*10+ch-48,ch=getchar(); return q ? -x : x;
}

il char gc(){ RG char ch=getchar(); while (ch<'A' || ch>'Z') ch=getchar(); return ch; }

il void Swap(RG int &a,RG int &b){ RG int t=a; a=b,b=t; return; }

il void insert(RG int from,RG int to){ g[++num]=(edge){head[from],to},head[from]=num; return; }

il void dfs1(RG int x,RG int p){
fa[x]=p,dep[x]=dep[p]+1,size[x]=1;
for (RG int i=head[x];i;i=g[i].nt){
RG int v=g[i].to; if (v==p) continue;
dfs1(v,x); size[x]+=size[v]; if (size[son[x]]<size[v]) son[x]=v;
}
return;
}

il void dfs2(RG int x,RG int p,RG int a){
dfn[++dex]=x,tid[x]=dex,top[x]=a; if (son[x]) dfs2(son[x],x,a);
for (RG int i=head[x];i;i=g[i].nt){ RG int v=g[i].to; if (v==p || v==son[x]) continue; dfs2(v,x,v); }
return;
}

il void build(RG int x,RG int l,RG int r){
if (l==r){ sum[x]=mx[x]=w[dfn[l]]; return; }
RG int mid=(l+r)>>1; build(ls,l,mid),build(rs,mid+1,r);
sum[x]=sum[ls]+sum[rs],mx[x]=max(mx[ls],mx[rs]); return;
}

il void update(RG int x,RG int l,RG int r,RG int p,RG int v){
if (l==r){ sum[x]=mx[x]=v; return; } RG int mid=(l+r)>>1;
if (p<=mid) update(ls,l,mid,p,v); else update(rs,mid+1,r,p,v);
sum[x]=sum[ls]+sum[rs],mx[x]=max(mx[ls],mx[rs]); return;
}

il int query1(RG int x,RG int l,RG int r,RG int xl,RG int xr){
if (xl<=l && r<=xr) return mx[x]; RG int mid=(l+r)>>1,res=0;
if (xr<=mid) res=query1(ls,l,mid,xl,xr);
else if (xl>mid) res=query1(rs,mid+1,r,xl,xr);
else res=max(query1(ls,l,mid,xl,mid),query1(rs,mid+1,r,mid+1,xr));
return res;
}

il int query2(RG int x,RG int l,RG int r,RG int xl,RG int xr){
if (xl<=l && r<=xr) return sum[x]; RG int mid=(l+r)>>1,res=0;
if (xr<=mid) res=query2(ls,l,mid,xl,xr);
else if (xl>mid) res=query2(rs,mid+1,r,xl,xr);
else res=query2(ls,l,mid,xl,mid)+query2(rs,mid+1,r,mid+1,xr);
return res;
}

il int Query1(RG int u,RG int v){
RG int res=-inf;
while (top[u]!=top[v]){
if (dep[top[u]]<dep[top[v]]) Swap(u,v);
res=max(res,query1(1,1,n,tid[top[u]],tid[u]));
u=fa[top[u]];
}
if (dep[u]>dep[v]) Swap(u,v); res=max(res,query1(1,1,n,tid[u],tid[v])); return res;
}

il int Query2(RG int u,RG int v){
RG int res=0;
while (top[u]!=top[v]){
if (dep[top[u]]<dep[top[v]]) Swap(u,v);
res+=query2(1,1,n,tid[top[u]],tid[u]);
u=fa[top[u]];
}
if (dep[u]>dep[v]) Swap(u,v); res+=query2(1,1,n,tid[u],tid[v]); return res;
}

il void work(){
n=gi(); for (RG int i=1;i<n;++i){ RG int u=gi(),v=gi(); insert(u,v),insert(v,u); }
for (RG int i=1;i<=n;++i) w[i]=gi(); dfs1(1,0),dfs2(1,0,1),build(1,1,n); q=gi();
for (RG int i=1;i<=q;++i){
RG char ch=gc();
if (ch=='C'){ RG int u=gi(),v=gi(); update(1,1,n,tid[u],v); }
if (ch=='Q'){ ch=gc(); RG int u=gi(),v=gi(); printf("%d\n",ch=='M' ? Query1(u,v) : Query2(u,v)); }
}
return;
}

int main(){
work();
return 0;
}