BZOJ1568: [JSOI2008]Blue Mary开公司

BZOJ1568: [JSOI2008]Blue Mary开公司

超哥线段树

题解：

超哥线段树裸题。

一张著名的图：

就是每个节点保存一个一次函数y，当插入新函数yy的时候，如果yy在整个区间都比y优，则替换成yy，反之保留y。当在区间中有交点的时候，哪个在终点位置的值好保留哪个，另一个往下推，推到交点所在的区间。

查询的时候一路取max即可。

Code：

#include <iostream>
#include <cstring>
#include <cstdio>
#include <cmath>
#define D(x) cout<<#x<<" = "<<x<<"  "
#define E cout<<endl
using namespace std;
const int N = 100005;
const int T = 50005;

int n;
double K[T*4],B[T*4]; 

bool better(double k1,double b1,double k2,double b2,double x){
    return k1*x+b1 > k2*x+b2;
}

void insert(int x,int l,int r,double k,double b){
    if(l==r){
        if(better(k,b,K[x],B[x],l)){ K[x]=k; B[x]=b; }
        return;
    }
    else{
        int mid=(l+r)>>1;
        if(k > K[x]){
            if(better(k,b,K[x],B[x],mid)){
                swap(K[x],k); swap(B[x],b);
                insert(x*2,l,mid,k,b);
            }
            else insert(x*2+1,mid+1,r,k,b);
        }
        else{
            if(better(k,b,K[x],B[x],mid)){
                swap(K[x],k); swap(B[x],b);
                insert(x*2+1,mid+1,r,k,b);
            }
            else insert(x*2,l,mid,k,b);
        }
    }
}

double query(int x,int l,int r,int p){
//  D(x); D(l); D(r); E;
    double d=K[x]*p+B[x];
//  D(d); E;
    if(l==r){ return d; }
    else{
        int mid=(l+r)>>1;
        if(p<=mid) return max(d,query(x*2,l,mid,p));
        else return max(d,query(x*2+1,mid+1,r,p));
    }
}

int main(){
    freopen("a.in","r",stdin);
    scanf("%d",&n);
    double b,k; int a; char op[15];
    for(int i=1;i<=n;i++){
        scanf("%s",op);
        if(op[0]=='P'){
            scanf("%lf%lf",&b,&k);
            b-=k;
//          D(k); D(b); E;
            insert(1,1,T,k,b);
        }
        else{
            scanf("%d",&a);
            printf("%d\n",(int)query(1,1,T,a)/100);
        }
    }
}