bzoj1568: [JSOI2008]Blue Mary开公司 （超哥线段树）

bzoj1568: [JSOI2008]Blue Mary开公司

Description

Input

第一行 ：一个整数N ，表示方案和询问的总数。
接下来N行，每行开头一个单词“Query”或“Project”。
若单词为Query，则后接一个整数T，表示Blue Mary询问第T天的最大收益。
若单词为Project，则后接两个实数S，P，表示该种设计方案第一天的收益S，以及以后每天比上一天多出的收益P。
1 <= N <= 100000 1 <= T <=50000 0 < P < 100，| S | <= 10^6
提示：本题读写数据量可能相当巨大，请选手注意选择高效的文件读写方式。

Output

对于每一个Query，输出一个整数，表示询问的答案，并精确到整百元（以百元为单位，
例如：该天最大收益为210或290时，均应该输出2）。没有方案时回答询问要输出0

Sample Input

10
Project 5.10200 0.65000
Project 2.76200 1.43000
Query 4
Query 2
Project 3.80200 1.17000、
Query 2
Query 3
Query 1
Project 4.58200 0.91000
Project 5.36200 0.39000

Sample Output

0
0
0
0
0

知识点：超哥线段树

很中二的名字有木有。
超哥线段树其实就是用线段树来维护半平面交的上/下凸壳。但是函数值只能是整数。
每个线段树节点中保存的是使这个区间中点的函数值达到最大/最小的直线，这里采用的是标记永久化的思想（下面为了方便统一最大值）
具体的方法是

1. 插入直线，并判断这条直线在区间中点处的函数值与原本线段树节点的函数值大小，如果更优就替换。
2. 比较当前直线和最优直线的斜率大小。如果斜率大，就往线段树右区间扔（因为此时它对左区间没有贡献），否则往左区间扔。
3. 到叶子节点比较值的大小并判断是否替换。

可以注意到，一条直线不一定是这个区间最大的那条直线，但是对于某个点或区间，其最值一定在到根路径的某条直线中。这就是标记永久化。所以说我们查询的方法也出来了。

1. 更新当前区间直线对于当前询问的答案。
2. 向下询问左右子区间。
3. 叶子节点直接返回当前直线答案。

代码

这道题是模板题，单点修改单点询问，不多说了。

#include<iostream>
#include<cstdlib>
#include<cstdio>
#include<cstring>
#include<algorithm>
#include<cmath>
using namespace std;
const int N = 55000;
const int M = 110000;
int read() {
    char ch = getchar(); int x = 0, f = 1;
    while(ch < '0' || ch > '9') {if(ch == '-') f = -1; ch = getchar();}
    while(ch >= '0' && ch <= '9') {x = (x << 1) + (x << 3) + ch - '0'; ch = getchar();}
    return x * f;
}
double a[M], b[M];
int tr[N << 2], rs[N << 2], ls[N << 2], sz, root, top;
double Ans(int pos, double x) {return a[pos] + b[pos] * (x - 1);}
bool Jud(int p1, int p2, int x) {return Ans(p1, x) > Ans(p2, x);}
void Change(int &p, int L, int R, int pos) {
    if(!p) p = ++sz;
    if(L == R) {
        if(Jud(pos, tr[p], L)) tr[p] = pos;
        return ;
    }
    int mid = L + R >> 1, tp = tr[p];
    if(b[pos] >= b[tp]) swap(tp, pos);
    if(Jud(pos, tp, mid)) tr[p] = pos, Change(rs[p], mid + 1, R, tp);
    else tr[p] = tp, Change(ls[p], L, mid, pos);
}
double Query(int p, int L, int R, int pos) {
    if(!p) return 0;
    if(L == R) return Ans(tr[p], pos);
    int mid = L + R >> 1;
    double ans = Ans(tr[p], pos);
    if(pos <= mid) ans = max(ans, Query(ls[p], L, mid, pos));
    else ans = max(ans, Query(rs[p], mid + 1, R, pos));
    return ans;
} 
  
int main() {
    int n = read();
    for(int i = 1;i <= n; ++i) {
        char ch = getchar();
        while(ch != 'Q' && ch != 'P') ch = getchar();
        char opt = getchar(); while(opt >= 'a' && opt <= 'z') opt = getchar();
        if(ch == 'Q') printf("%d\n", (int)Query(root, 1, N, read()) / 100);
        else {
            ++top; scanf("%lf%lf", &a[top], &b[top]);
            Change(root, 1, N, top);
        }
    }
    return 0;
}