BZOJ1568: [JSOI2008]Blue Mary开公司（洛谷P4254）

线段树

BZOJ题目传送门
洛谷题目传送门

这是一个叫李超线段树的东西。

这道题可以抽象成一个平面直角坐标系，有两种操作：加入一条直线或询问当$x=t$时所有直线的最大值（$t$为整数）。

对于每个区间记录一条直线$k_x$和$b_x$，在这个区间中这条直线取到的长度最长。

当插入一条直线时，对于一个区间分三种情况讨论：

1.没有最长直线或当前直线完全覆盖最长直线。此时直接把当前直线当成最长直线即可。

2.最长直线完全覆盖当前直线。此时不做修改直接返回。

3.这两条直线有交点$x_0$。那么根据$x_0$与区间中点的位置关系和两条直线在$x_0$左右的大小关系，最长直线取被覆盖少的那条，另一条递归做它比较大的那半个区间。

查询的话把当前区间的直线在$t$取到的值和递归后返回的值比较，返回大的那个即可。

代码：

#include<cctype>
#include<cstdio>
#include<cstring>
#include<algorithm>
#define N 50005
#define F inline
using namespace std;
typedef double DB;
typedef long long LL;
struct tree{ int l,r,f; DB k,b; }t[N<<2];
int n=50000,m; char s[10];
F char readc(){
	static char buf[100000],*l=buf,*r=buf;
	if (l==r) r=(l=buf)+fread(buf,1,100000,stdin);
	return l==r?EOF:*l++;
}
F DB _read(){
	DB x=0,f=1,p=1; char ch=readc();
	while (!isdigit(ch)&&!isupper(ch)) ch=='-'?f=-1:1,ch=readc();
	if (isupper(ch)) return ch;
	while (isdigit(ch)&&ch!='.') x=x*10+ch-48,ch=readc();
	if (ch!='.') return x*f; ch=readc();
	while (isdigit(ch)) x=x*10+ch-48,p*=10,ch=readc();
	return x*f/p;
}
void build(int x,int l,int r){
	t[x].l=l,t[x].r=r; if (l==r) return; int mid=l+r>>1;
	build(x<<1,l,mid),build(x<<1|1,mid+1,r);
}
void mdfy(int x,DB k,DB b){
	if (!t[x].f) return t[x].f=1,t[x].k=k,void(t[x].b=b);
	int mid=t[x].l+t[x].r>>1;
	DB l1=k*t[x].l+b,l2=t[x].k*t[x].l+t[x].b;
	DB r1=k*t[x].r+b,r2=t[x].k*t[x].r+t[x].b;
	if (l1<=l2&&r1<=r2) return;
	if (l1>l2&&r1>r2) return t[x].k=k,void(t[x].b=b);
	DB xn=(b-t[x].b)/(t[x].k-k),kk=k,bb=b; bool f1=l1>l2,f2=xn>mid;
	if (!(f1^f2)) kk=t[x].k,bb=t[x].b,t[x].k=k,t[x].b=b; mdfy(x<<1|f2,kk,bb);
}
DB srch(int x,int p){
	DB ans=t[x].k*p+t[x].b; int mid=t[x].l+t[x].r>>1;
	if (t[x].l==t[x].r) return ans;
	return max(ans,srch(x<<1|(p>mid),p));
}
int main(){
	for (build(1,1,n),m=_read();m;m--){		
		if (_read()=='Q') printf("%lld\n",LL(srch(1,_read())/100.0));
		else{ DB b=_read(),k=_read(); mdfy(1,k,b-k); }
	}
	return 0;
}