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

我又犯了個zz的錯誤，
第i個布娃娃想知道喜歡它的布娃娃中，魅力值第i大的布娃娃的魅力值是多少
要找的是第k**大**

爆炸oj日常爆炸

這里寫代碼片 #include<cstdio> #include<cstring> #include<iostream> #include<algorithm>using 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<b.x;else return a.ff<b.ff; //先刪后加 }int cmp2(const po &a,const po &b) {return a.x<b.x; }void update(int bh) {t[bh].s=t[bh<<1].s+t[bh<<1|1].s; }void build(int bh,int l,int r) //權值線段樹（魅力值C） {t[bh].l=l; t[bh].r=r;t[bh].ml=C[l]; t[bh].mr=C[r];if (l==r) return;int mid=(l+r)>>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 0;if (t[bh].l==t[bh].r) return t[bh].ml;if (t[bh<<1|1].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; }

轉載于:https://www.cnblogs.com/wutongtong3117/p/7673317.html

總結

以上是生活随笔為你收集整理的bzoj2161 布娃娃的全部內容，希望文章能夠幫你解決所遇到的問題。

如果覺得生活随笔網站內容還不錯，歡迎將生活随笔推薦給好友。