​https://www.lydsy.com/JudgeOnline/problem.php?id=1568​

李超线段树模板


#include <cstdio>
#include <cstring>
#include <algorithm>
using namespace std;
#define N 500010

double a[1000010],b[1000010];
int tr[2000010];
int n,m;

double gety(int k,int x)
{
return a[k]+(x-1)*b[k];
}

int judge(int x,int y,int pos)
{
return gety(x,pos)>gety(y,pos);
}

void update(int now,int l,int r,int x)
{
int m;
if(l==r)
{
if(judge(x,tr[now],l)) tr[now]=x;
return;
}
m=(l+r)/2;
if(b[x]>b[tr[now]])
{
if(judge(x,tr[now],m))
{
update(now<<1,l,m,tr[now]);
tr[now]=x;
}
else update(now<<1|1,m+1,r,x);
}
if(b[x]<b[tr[now]])
{
if(judge(x,tr[now],m))
{
update(now<<1|1,m+1,r,tr[now]);
tr[now]=x;
}
else update(now<<1,l,m,x);
}
}

double query(int now,int l,int r,int x)
{
double res;
int m;
if(l==r) return gety(tr[now],x);
m=(l+r)/2;
res=gety(tr[now],x);
if(x<=m) res=max(res,query(now<<1,l,m,x));
else res=max(res,query(now<<1|1,m+1,r,x));
return res;
}

int main()
{
double t;
int q,k,x;
char op[10];
scanf("%d",&q);
while(q--)
{
scanf("%s",op);
if(op[0]=='P')
{
m++;
scanf("%lf%lf",&a[m],&b[m]);
update(1,1,N,m);
}
else
{
scanf("%d",&x);
t=query(1,1,N,x);
k=t;
printf("%d\n",k/100);
}
}
return 0;
}