扫描线|bzoj2161 布娃娃扫描线|线段树

Description
小时候的雨荨非常听话，是父母眼中的好孩子。在学校是老师的左右手，同学的好榜样。后来她成为艾利斯顿第二代考神，这和小时候培养的良好素质是分不开的。雨荨的妈妈也为有这么一个懂事的女儿感到高兴。一次期末考试，雨荨不知道第多少次，再次考了全年级第一名。雨荨的妈妈看到女儿100分的成绩单时，脸上又泛起了幸福的笑容，作为奖励，她给雨荨买了n个布娃娃。细心的雨荨发现，第i个布娃娃有一个耐心值P[i]以及一个魅力值C[i]，并且还有能够忍受的耐心值的上限R[i]以及下限L[i]。当一个布娃娃j满足L[j]<=P[i]并且P[i]<=R[j]，那么布娃娃j喜欢布娃娃i。雨荨还发现，一个布娃娃有可能喜欢它自己。每个布娃娃心中都有一个谜团，具体来说就是：第i个布娃娃想知道喜欢它的布娃娃中，魅力值第i大的布娃娃的魅力值是多少，并且称这个布娃娃的谜团答案为这个魅力值的大小，如果不存在，那么这个布娃娃的谜团答案为0。鉴于雨荨的上司栋栋不让题目的数据过大，下面给出数据的生成方法：给出16个参数：
Padd, Pfirst, Pmod, Pprod, Cadd, Cfirst, Cmod, Cprod, Ladd, Lfirst, Lmod, Lprod, Radd, Rfirst, Rmod, Rprod。
P[1] = Pfirst % Pmod, P[i] = (P[i-1] Pprod + Padd + i) % Pmod (i > 1)。
对于C、L、R数组也有类似的得到方式, %代表取余运算。注意：L和R数组生成完之后，如果某个布娃娃的忍耐度上限小于下限，那么交换它的上限和下限。当然，雨荨也不会让你告诉她每个布娃娃的谜团答案，因为那样会使输出数据很大。所以雨荨希望你告诉她，所有布娃娃谜团答案的和除以19921228的余数是多少。
Input
输入的第一行有一个整数n，代表布娃娃的个数。
输入的第二行有16个用空格隔开的整数
分别代表Padd,Pfirst,Pmod,Pprod,Cadd,Cfirst,Cmod,Cprod,Ladd,Lfirst,Lmod,Lprod,Radd,Rfirst,Rmod,Rprod。
16个参数均为1到100,000,000中的整数。
Output
输出一个整数，代表所有布娃娃谜团答案的和除以19921228的余数。
Sample Input
3
2 3 4 3 1 4 5 2 3 6 9 1 1 2 3 4
Sample Output
4
题面补充：n<=10^5
分析：
一眼主席树
但是需要注意的是“喜欢ta”的布娃娃，不是“ta喜欢”的布娃娃
简化了一下题意
给出一堆点，和一堆线段
求覆盖某一点的所有线段中权值第k大的线段
我们把每条线段分成加入和删除两个操作
按照x排序
建立权值线段树，记录C（魅力值）
每个线段树节点中，我们维护几个信息
l,r 左右端点
ml,mr 实际的C值
s在这个节点所管辖的区间内有几个数
P按照从小到大排序
从左向右扫
如果当前的线段端点 < 当前查找的P
则加入权值线段树
在查找的时候，如果当前节点的s根本达不到k，直接return 0
如果>=k我们就分成左右两个节点
如果右节点的s <= k，直接进右儿子
否则 k-=右儿子的s，进入左儿子
（有点像主席树）
我们在这里学到了一个新知识，建立权值线段树查找第k大的时候
需要记录一下区间内的元素个数，
进行类似主席树的查找
tip 注意：L和R数组生成完之后，
如果某个布娃娃的忍耐度上限小于下限，那么交换它的上限和下限
然而第一次写的时候，我是一边写一边换的，jj
将P排序时，要保存一下原先的排名
因为原先的排名就是查找是需要用的k
【扫描线|bzoj2161 布娃娃】我又犯了个zz的错误，
第i个布娃娃想知道喜欢它的布娃娃中，魅力值第i大的布娃娃的魅力值是多少
要找的是第k**大**
爆炸oj日常爆炸

这里写代码片 #include #include #include #includeusing namespace std; const int mod=19921228; const int N=100010; struct node{ int x,y,ff; }; node li[N*3]; int Padd,Pfirst,Pmod,Pprod,Cadd,Cfirst,Cmod,Cprod,Ladd,Lfirst,Lmod,Lprod,Radd,Rfirst,Rmod,Rprod; int L[N],R[N],C[N],n,tot=0; struct nd{ int l,r,ml,mr,s; }; nd t[N<<3]; struct po{ int x,id; }; po P[N]; int cmp(const node &a,const node &b) { if (a.x!=b.x) return a.x>1; build(bh<<1,l,mid); build(bh<<1|1,mid+1,r); }void add(int bh,int mb,int val) { if (t[bh].l==t[bh].r) { t[bh].s+=val; return; } int mid=(t[bh].l+t[bh].r)>>1; if (mb<=mid) add(bh<<1,mb,val); else add(bh<<1|1,mb,val); update(bh); }int ask(int bh,int k) { if (t[bh].s=k) return ask(bh<<1|1,k); else return ask(bh<<1,k-t[bh<<1|1].s); }int main() { scanf("%d",&n); scanf("%d%d%d%d%d%d%d%d%d%d%d%d%d%d%d%d",&Padd,&Pfirst,&Pmod,&Pprod,&Cadd,&Cfirst,&Cmod,&Cprod,&Ladd,&Lfirst,&Lmod,&Lprod,&Radd,&Rfirst,&Rmod,&Rprod); P[1].x=Pfirst%Pmod; C[1]=Cfirst%Cmod; P[1].id=1; L[1]=Lfirst%Lmod; R[1]=Rfirst%Rmod; for (int i=2; i<=n; i++) { P[i].x=(P[i-1].x*Pprod+Padd+i)%Pmod; C[i]=(C[i-1]*Cprod+Cadd+i)%Cmod; P[i].id=i; L[i]=(L[i-1]*Lprod+Ladd+i)%Lmod; R[i]=(R[i-1]*Rprod+Radd+i)%Rmod; }for (int i=1; i<=n; i++) { if (L[i]>R[i]) swap(L[i],R[i]); tot++; li[tot].x=L[i]; li[tot].y=C[i]; li[tot].ff=1; tot++; li[tot].x=R[i]+1; li[tot].y=C[i]; li[tot].ff=-1; }sort(li+1,li+1+tot,cmp); sort(C+1,C+1+n); sort(P+1,P+1+n,cmp2); build(1,1,n); int j=1; int ans=0; for (int i=1; i<=n; i++) { while (li[j].x<=P[i].x) add(1,li[j].y,li[j].ff),j++; ans+=ask(1,P[i].id); ans%=mod; } printf("%d",ans); return 0; }

扫描线|bzoj2161 布娃娃

推荐阅读

用CDR矩形工具设计箭头的方法

江西新余有多少个县

根轨迹分析法,时域分析法

烟熏妆一般在什么场合会画呢如何画烟熏妆呢

光油和清漆一样吗、有何区别？光油作用是什么？

IOS系统中好玩的手机游戏和实用的软件都有哪些？

地震后会不会有更大的地震地震后是否会有更大的地震

全文 2023年武汉社保最低缴费标准基数及比例武汉个人社保缴费价格表

抖音直播流行什么饰品卖，抖音直播流行什么饰品卖的

微信二

mysql将数字转换为时间格式，mysql把时间转换为字符串

职务侵占罪请律师要花多少钱职务侵占罪请律师要花多少钱

盆栽的菊花有什么好处，养菊花的好处和坏处

新媒体运营数据分析

喝黄芪泡水一次放多少

大麦若叶青汁适合什么人群喝

鹿茸泡水喝的功效与作用鹿茸泡水喝的禁忌

赛睿gg驱动怎么用赛睿鼠标驱动，赛睿310驱动设置教程

月球与地球的平均距离约38 月亮表面温度约多少摄氏度

metleb对数据分析