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