BZOJ 1568: [JSOI2008]Blue Mary开公司

李超线段树裸题，每个区间维护覆盖当前区间中线的线段，发现答案必然在这些线段中。

复杂度听说是两只log，但感觉实际上只有一只？

#include<bits/stdc++.h>
#define N 100100
using namespace std;
double s[N<<2],p[N<<2];
int done[N<<2];
int n;
char op[10];
double cross(double b1,double b2,double k1,double k2){
	return (b2-b1)/(k1-k2);
}
void insert(double a,double b,int l,int r,int k){
	if(!done[k]){
		done[k]=1;
		s[k]=a,p[k]=b;
		return;
	}
	double f1=s[k]+p[k]*l,f2=s[k]+p[k]*r;
	double f3=a+b*l,f4=a+b*r;
	int mid=l+r>>1;
	if(f1<=f3&&f2<=f4) s[k]=a,p[k]=b;
	else if(f3<=f1&&f4<=f2) return;
	else{
		double X=cross(s[k],a,p[k],b);
		if(f1>=f3){
			if(X<=mid) insert(s[k],p[k],l,mid,k<<1),s[k]=a,p[k]=b;
			else insert(a,b,mid+1,r,k<<1|1);
		}
		else{
			if(X<=mid) insert(a,b,l,mid,k<<1);
			else insert(s[k],p[k],mid+1,r,k<<1|1),s[k]=a,p[k]=b;
		}
	}
}
double query(int t,int l,int r,int k){
	double bes=0;
	if(done[k]) bes=max(bes,s[k]+t*p[k]);
	if(l==r) return bes;
	int mid=l+r>>1;
	if(t<=mid) bes=max(bes,query(t,l,mid,k<<1));
	else bes=max(bes,query(t,mid+1,r,k<<1|1));
	return bes;
}
int main(){
	scanf("%d",&n);
	for(int i=0;i<n;i++){
		int t;double a,b;
		scanf("%s",op);
		if(op[0]=='Q'){
			scanf("%d",&t);
			printf("%.0lf\n",floor(query(t,1,50000,1)/100));
		}
		else{
			scanf("%lf%lf",&a,&b);
			insert(a-b,b,1,50000,1);
		}
	}
	return 0;
}