HDU-1823-Luck and Love-二维线段树入门

http://acm.hust.edu.cn/vjudge/problem/viewProblem.action?id=14799

题意：

Luck and Love

Time Limit: 1000MS Memory Limit: 32768KB 64bit IO Format: %I64d & %I64u

Submit Status

Description

世界上上最远的距离不是相隔天涯海角 
而是我在你面前 
可你却不知道我爱你 
                �� 张小娴 

前段日子，枫冰叶子给Wiskey做了个征婚启事，聘礼达到500万哦，天哪，可是天文数字了啊，不知多少MM蜂拥而至，顿时万人空巷，连扫地的大妈都来凑热闹来了。�_�||| 
由于人数太多，Wiskey实在忙不过来，就把统计的事情全交给了枫冰叶子，自己跑回家休息去了。这可够枫冰叶子忙的了，他要处理的有两类事情，一是得接受MM的报名，二是要帮Wiskey查找符合要求的MM中缘分最高值。 

 

Input

本题有多个测试数据，第一个数字M，表示接下来有连续的M个操作，当M＝0时处理中止。 
接下来是一个操作符C。 
当操作符为‘I’时，表示有一个MM报名，后面接着一个整数，H表示身高，两个浮点数，A表示活泼度，L表示缘分值。 （100<=H<=200， 0.0<=A，L<=100.0） 
当操作符为‘Q’时，后面接着四个浮点数，H1，H2表示身高区间，A1，A2表示活泼度区间，输出符合身高和活泼度要求的MM中的缘分最高值。 （100<=H1，H2<=200， 0.0<=A1，A2<=100.0） 
所有输入的浮点数，均只有一位小数。 

 

Output

对于每一次询问操作，在一行里面输出缘分最高值，保留一位小数。 
对查找不到的询问，输出-1。 

 

用二维线段树，母树上维护的是高度区间对应的子树

 特定高度区间对应一个子树，子树上维护的是 每个活跃度区间的最大缘份值

每次操作复杂度 logn*logm  n,m分别是子树母树的节点数

#include <cstdio>#include <cmath>#include <cstring>#include <string>#include <algorithm>#include <iostream>#include <queue>#include <map>#include <set>#include <vector>using namespace std;const int sub_N=1000*4;const int N=200*4;int max(int a,int b){return a<b?b:a;}struct  sub_node //子线段树节点{int sl,sr;//sub_l,sub_rint max;  //附加信息};struct node //母线段树节点{int l,r;sub_node T[sub_N]; };node TT[N];void sub_build(int rt,int sub_rt,int sl,int sr){TT[rt].T[sub_rt].sl=sl;TT[rt].T[sub_rt].sr=sr;TT[rt].T[sub_rt].max=-1;if (sr==sl) return;int mid=(sl+sr)>>1;sub_build(rt,sub_rt<<1,sl,mid);sub_build(rt,sub_rt<<1|1,mid+1,sr);}void build(int rt,int l,int r,int sl,int sr){TT[rt].l=l;TT[rt].r=r;sub_build(rt,1,sl,sr);if (l==r) return ;int mid=(l+r)>>1;build(rt<<1,l,mid,sl,sr);build(rt<<1|1,mid+1,r,sl,sr); }void sub_update(int rt,int sub_rt,int active,int love)  //在当前【高度区间母树枝】下 更新活跃度为active的子分支{if (TT[rt].T[sub_rt].sl==TT[rt].T[sub_rt].sr){TT[rt].T[sub_rt].max=max(TT[rt].T[sub_rt].max,love);return ;}int sl=TT[rt].T[sub_rt].sl;int sr=TT[rt].T[sub_rt].sr;int mid=(sl+sr)>>1;if (active<=mid)sub_update(  rt, sub_rt<<1, active,  love);elsesub_update(  rt, sub_rt<<1|1, active,  love);TT[rt].T[sub_rt].max=max(TT[rt].T[sub_rt].max,love);//当前rt管辖的活跃度区间包含了active,所以更新max值}void update(int rt,int hh,int active,int love)//寻找合适的高度区间的母树枝更新{sub_update(rt,1,active,love);if (TT[rt].l==TT[rt].r)  return ;int mid=(TT[rt].l+TT[rt].r)>>1;if (hh<=mid)update(rt<<1,hh,active,love);elseupdate(rt<<1|1,hh,active,love);} int sub_query(int rt,int sub_rt,int sl,int sr){if (TT[rt].T[sub_rt].sl>=sl &&TT[rt].T[sub_rt].sr<=sr)return TT[rt].T[sub_rt].max;int mid=(TT[rt].T[sub_rt].sl+TT[rt].T[sub_rt].sr)>>1;int ret1=-1,ret2=-1;if (sl<=mid)ret1=sub_query(rt,sub_rt<<1,sl,sr);if (sr>mid)ret2=sub_query(rt,sub_rt<<1|1,sl,sr); return max(ret1,ret2);}int  query(int rt,int h1,int h2,int sl,int sr){if (TT[rt].l>=h1 &&TT[rt].r<=h2)return sub_query(rt,1,sl,sr);int mid=(TT[rt].l+TT[rt].r)>>1;int ret1=-1,ret2=-1;if (h1<=mid)ret1=query(rt<<1,h1,h2,sl,sr);if (h2>mid)ret2=query(rt<<1|1,h1,h2,sl,sr);return max(ret1,ret2);} int main(){double active,love;int hh;int n,i;while(cin>>n){getchar();build(1,100,200,0,1000);char op;int h1,h2;double a1,a2;for (i=1;i<=n;i++){scanf("%c",&op);if (op=='I'){scanf("%d%lf%lf",&hh,&active,&love); int aa=(int)(active*10);int ll=(int)(love*10);update(1,hh,aa,ll); }else{scanf("%d%d%lf%lf",&h1,&h2,&a1,&a2);int aa1=(int)(a1*10);int aa2=(int)(a2*10);if (h1>h2)swap(h1,h2);if (aa1>aa2)swap(aa2,aa1);int ans=query(1,h1,h2,aa1,aa2);if (ans==-1)printf("-1\n");elseprintf("%.1lf\n",((double)ans)/10.0);}getchar();}}return 0;}