http://acm.hust.edu.cn/vjudge/problem/viewProblem.action?id=14799
题意:
Time Limit: 1000MS | Memory Limit: 32768KB | 64bit IO Format: %I64d & %I64u |
Description
世界上上最远的距离不是相隔天涯海角
而是我在你面前
可你却不知道我爱你
�� 张小娴
前段日子,枫冰叶子给Wiskey做了个征婚启事,聘礼达到500万哦,天哪,可是天文数字了啊,不知多少MM蜂拥而至,顿时万人空巷,连扫地的大妈都来凑热闹来了。�_�|||
由于人数太多,Wiskey实在忙不过来,就把统计的事情全交给了枫冰叶子,自己跑回家休息去了。这可够枫冰叶子忙的了,他要处理的有两类事情,一是得接受MM的报名,二是要帮Wiskey查找符合要求的MM中缘分最高值。
而是我在你面前
可你却不知道我爱你
�� 张小娴
前段日子,枫冰叶子给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)
所有输入的浮点数,均只有一位小数。
接下来是一个操作符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。
对查找不到的询问,输出-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_r
int 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);
else
sub_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);
else
update(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");
else
printf("%.1lf\n",((double)ans)/10.0);
}
getchar();
}
}
return 0;
}