李超线段树
模板題:[JSOI2008]Blue Mary開公司https://www.luogu.org/problemnew/show/P4254
李超線段樹:動態維護一個平面直角坐標系,支持在中間插入一條線段(或直線),詢問與x=x0這條直線相交的所有線段中,交點的y坐標的最大(小)值。
線段樹上每個區間維護在mid處y坐標最大的線段。(要標記永久化)
考慮怎么插入一條直線,假設它當前處理到了某個區間:(摘自yyb的博客https://www.cnblogs.com/cjyyb/p/10560973.html)
- 如果這個區間沒有記錄最長的線段,那么直接把這個區間記錄的線段修改為這條線段,然后返回。
- 如果當前線段在這個區間內已經被這個區間記錄的線段覆蓋,那么直接gg,返回就好了。
- 反過來,如果完全覆蓋了之前記錄的線段,那么直接賦值、返回。
- 否則和已經記錄的直線有交,判斷哪根線段覆蓋的區域較長,把這個區間記錄的值給修改一下,然后把短的那一半丟下去遞歸。 #include<cstdio>
#include<cstring>
#include<algorithm>
using namespace std;
const int N=100005;
int n,i;
bool v[N*5];
double x,y,K[N*5],B[N*5];
char ch[20];
void update(int o,int l,int r,double k,double b)
{if(!v[o]){v[o]=true,K[o]=k,B[o]=b;return;}if(l==r){if(k*l+b>K[o]*l+B[o])K[o]=k,B[o]=b;return;}int mid=l+r>>1;double l1=k*l+b,r1=k*r+b,l2=K[o]*l+B[o],r2=K[o]*r+B[o];if(l1>=l2&&r1>=r2){K[o]=k,B[o]=b;return;}if(l1<=l2&&r1<=r2)return;double x=(b-B[o])/(K[o]-k);if(l1>=l2){if(x>mid){update(o<<1|1,mid+1,r,K[o],B[o]);K[o]=k,B[o]=b;}elseupdate(o<<1,l,mid,k,b);}else{if(x>mid)update(o<<1|1,mid+1,r,k,b);else{update(o<<1,l,mid,K[o],B[o]);K[o]=k,B[o]=b;}}
}
double query(int o,int l,int r,int x)
{if(l==r)return K[o]*l+B[o];int mid=l+r>>1;double rtn=(v[o]?K[o]*x+B[o]:0.0);if(x<=mid)return max(rtn,query(o<<1,l,mid,x));elsereturn max(rtn,query(o<<1|1,mid+1,r,x));
}
int main()
{freopen("company.in","r",stdin);freopen("company.out","w",stdout);scanf("%d",&n);while(n--){scanf("%s",ch);if(ch[0]=='P'){scanf("%lf%lf",&x,&y);update(1,1,50000,y,x-y);}else{scanf("%d",&i);printf("%d\n",max(0,(int)query(1,1,50000,i)/100));}}return 0;
}
?這題因為插入的是直線,所以插入、查詢都是log(n)的。
轉載于:https://www.cnblogs.com/pthws/p/11079795.html
總結
- 上一篇: d3.js 简易柱形图,入门demo
- 下一篇: 学习笔记84—[深度学习]神经网络反向传