HDU 1823 Luck and Love 二维线段树

题目：

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;}