学步园

题意：Wiskey招女友，每个女生看其身高、活泼度和缘分值。现在执行两种操作，1、I，加入一位女生的身高，活泼度和缘分值；2、Q，查询身高在H1, H2之间，活泼度在A1, A2之间的女生的最高缘分值。

题目链接：http://acm.hdu.edu.cn/showproblem.php?pid=1823

——>>查询某个区间的最值，若是一维，可用RMQ解法。。也可用线段树解法。。

现在要查身高限一个区间，活泼度限一个区间，是一个二维的情景。。将线段树扩至二维。。时间复杂度：O(N+M*log(N)^2)。。

——>>坑1：G++的BUG。。。。以G++提交数次皆WA。。改交C++即过。。有朋友指出是代码触发了未定义行为，也有朋友说是因为G++的O2 BUG。。

——>>坑2：题目中说是1位小数，那么，处理方法用scanf("%d.%d", &A1, &A2)，再进行A1 * 10 + A2，会比用scanf("%lf", &A)，再进行(int)(A * 10)更精确（自己测下0.0到0.9就可以发现）。。可是，用第一种处理方式却会WA，偏偏要用第二种不那么精确的处理方式才会AC。。

第一道二维线段树题目，做到想哭了。。

#include <cstdio>
#include <algorithm>

using namespace std;

#define lc (o<<1)
#define rc ((o<<1)|1)

const int maxh = 100 + 10;
const int maxa = 1000 + 10;
int mmax[maxh<<2][maxa<<2];

void build_A(int O, int o, int L, int R) {      //活泼度建树
    mmax[O][o] = -1;
    if(L == R) return;
    int M = (L + R) >> 1;
    build_A(O, lc, L, M);
    build_A(O, rc, M+1, R);
}

void build_H(int o, int L, int R) {     //身高建树
    build_A(o, 1, 0, 1000);
    if(L == R) return;
    int M = (L + R) >> 1;
    build_H(lc, L, M);
    build_H(rc, M+1, R);
}

void Insert_A(int O, int o, int L, int R, int A, int D) {       //根据活泼度建树
    if(L == R) {
        mmax[O][o] = max(mmax[O][o], D);
        return;
    }
    int M = (L + R) >> 1;
    if(A <= M) Insert_A(O, lc, L, M, A, D);
    else Insert_A(O, rc, M+1, R, A, D);
    mmax[O][o] = max(mmax[O][lc], mmax[O][rc]);
}

void Insert(int o, int L, int R, int H, int A, int D) {
    Insert_A(o, 1, 0, 1000, A, D);
    if(L == R) return;
    int M = (L + R) >> 1;
    if(H <= M) Insert(lc, L, M, H, A, D);
    else Insert(rc, M+1, R, H, A, D);
}

int query_A(int O, int o, int L, int R, int A1, int A2) {       //根据活泼度查询
    if(A1 <= L && R <= A2) return mmax[O][o];
    int M = (L + R) >> 1;
    int Max1 = -1, Max2 = -1;
    if(A1 <= M) Max1 = query_A(O, lc, L, M, A1, A2);
    if(A2 > M) Max2 = query_A(O, rc, M+1, R, A1, A2);
    return (Max1 > Max2) ? Max1 : Max2;
}

int query(int o, int L, int R, int H1, int H2, int A1, int A2) {
    if(H1 <= L && R <= H2) return query_A(o, 1, 0, 1000, A1, A2);
    int M = (L + R) >> 1;
    int Max1 = -1, Max2 = -1;
    if(H1 <= M) Max1 = query(lc, L, M, H1, H2, A1, A2);
    if(H2 > M) Max2 = query(rc, M+1, R, H1, H2, A1, A2);
    return (Max1 > Max2) ? Max1 : Max2;
}

int main()
{
    int M;
    char op;
    while(scanf("%d", &M) == 1 && M) {
        build_H(1, 100, 200);
        for(int i = 0; i < M; i++) {
            getchar();
            op = getchar();
            if(op == 'I') {
//                int H, A, A1, A2, D, D1, D2;
//                scanf("%d %d.%d %d.%d", &H, &A1, &A2, &D1, &D2);
//                A = A1 * 10 + A2;
//                D = D1 * 10 + D2;
                int H, A, D;
                double AA, DD;
                scanf("%d%lf%lf", &H, &AA, &DD);
                A = (int)(AA * 10);
                D = (int)(DD * 10);
                Insert(1, 100, 200, H, A, D);
            }
            else {
//                int H1, H2, A[6];
//                scanf("%d %d %d.%d %d.%d", &H1, &H2, A+2, A+3, A+4, A+5);
//                A[0] = A[2] * 10 + A[3];
//                A[1] = A[4] * 10 + A[5];
                int H1, H2, A1, A2;
                double AA, BB;
                scanf("%d%d%lf%lf", &H1, &H2, &AA, &BB);
                A1 = (int)(AA * 10);
                A2 = (int)(BB * 10);
                if(A1 > A2) swap(A1, A2);
                if(H1 > H2) swap(H1, H2);
//                if(A[0] > A[1]) swap(A[0], A[1]);
//                int ret = query(1, 100, 200, H1, H2, A[0], A[1]);
                int ret = query(1, 100, 200, H1, H2, A1, A2);
                ret != -1 ? printf("%.1f\n", ret/10.0) : puts("-1");
            }
        }
    }
    return 0;
}