题目：
http://acm.hdu.edu.cn/showproblem.php?pid=1823
题意：
Problem 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。
Sample Input
8
I 160 50.5 60.0
I 165 30.0 80.5
I 166 10.0 50.0
I 170 80.5 77.5
Q 150 166 10.0 60.0
Q 166 177 10.0 50.0
I 166 40.0 99.9
Q 166 177 10.0 50.0
0
Sample Output
80.5
50.0
99.9
思路
二维线段树，第一维维护高度，第二维维护活泼度，注意精度问题。注意若有两个女孩身高活泼度一样的话， 此时缘分值应该取其中较大的，而不是直接赋值（我比较蠢。。。）。另外这个题的活泼度数据貌似有问题，最大值可能是100.8，超出了题目所给范围，WA了一上午。。。
#include <iostream>
#include <cstdio>
#include <cstring>
#include <algorithm>
#include <functional>
using namespace std;
const int N = 200 + 10, M = 1000 + 50, INF = 0x3f3f3f3f;
const double eps = 1e-6;
int maxv[N<<2][M<<2];
int n = 200 + 5, m = 1000 + 20;
void push_upy(int ky, int kx)
{
maxv[kx][ky] = max(maxv[kx][ky<<1], maxv[kx][ky<<1|1]);
}
void push_upx(int ky, int kx)
{
maxv[kx][ky] = max(maxv[kx<<1][ky], maxv[kx<<1|1][ky]);
}
void updatey(int y, int L, int R, int ky, int kx, int val, int f)
{
if(L == R)
{
if(f) maxv[kx][ky] = max(maxv[kx][ky], val);
else push_upx(ky, kx);
return;
}
int mid = (L + R) >> 1;
if(y <= mid) updatey(y, L, mid, ky << 1, kx, val, f);
else updatey(y, mid + 1, R, ky << 1|1, kx, val, f);
push_upy(ky, kx);
}
void updatex(int x, int y, int L, int R, int kx, int val)
{
if(L == R)
{
updatey(y, 1, m, 1, kx, val, 1); return;
}
int mid = (L + R) >> 1;
if(x <= mid) updatex(x, y, L, mid, kx << 1, val);
else updatex(x, y, mid + 1, R, kx << 1|1, val);
updatey(y, 1, m, 1, kx, val, 0);
}
int queryy(int ly, int ry, int L, int R, int ky, int kx)
{
if(ly <= L && R <= ry) return maxv[kx][ky];
int mid = (L + R) >> 1;
int ans = -1;
if(ly <= mid) ans = max(ans, queryy(ly, ry, L, mid, ky << 1, kx));
if(ry > mid) ans = max(ans, queryy(ly, ry, mid + 1, R, ky << 1|1, kx));
return ans;
}
int queryx(int lx, int rx, int ly, int ry, int L, int R, int kx)
{
if(lx <= L && R <= rx) return queryy(ly, ry, 1, m, 1, kx);
int mid = (L + R) >> 1;
int ans = -1;
if(lx <= mid) ans = max(ans, queryx(lx, rx, ly, ry, L, mid, kx << 1));
if(rx > mid) ans = max(ans, queryx(lx, rx, ly, ry, mid + 1, R, kx << 1|1));
return ans;
}
int main()
{
int t;
while(scanf("%d", &t), t)
{
memset(maxv, -1, sizeof maxv);
char op;
for(int i = 1; i <= t; i++)
{
scanf(" %c", &op);
if(op == 'I')
{
int a;
double b, c;
scanf("%d%lf%lf", &a, &b, &c);
b += 1;
updatex(a, int(b*10 + eps), 1, n, 1, int(c*10 + eps));
}
else
{
double a, b, c, d;
scanf("%lf%lf%lf%lf", &a, &b, &c, &d);
int ta = (int)ceil(a), tb = (int)floor(b);
c += 1, d += 1;
if(ta > tb) swap(ta, tb);
if(c > d) swap(c, d);
int ans = queryx(ta, tb, int(c*10 + eps), int(d*10 + eps), 1, n, 1);
if(ans == -1) printf("-1\n");
else printf("%.1f\n", 1.0 * ans / 10);
}
}
}
return 0;
}
