青青草原综合久久大伊人导航_色综合久久天天综合_日日噜噜夜夜狠狠久久丁香五月_热久久这里只有精品

算法學社
記錄難忘的征途
posts - 141,comments - 220,trackbacks - 0

題目描述:

   對一個序列進行維護,要求支持四種操作:
      1. 在結尾加入一個數。
      2. 詢問區間第K大的數
      3. 詢問大小為X的數在序列中的排名
      4. 詢問第K大的數

吐槽:

    1. 真是捉了N久啊。。。。。。
    2. 手動寫內存池,與MLE只差400K。。。。。

算法分析:

    操作3和4的話,離散化一下用樹狀數組/線段樹搞就可以了。問題如何維護1和2。
    好在更改的操作只是在結尾加入一個數,所以離線處理一下,用一個支持查詢區間第K大的數據結構就行了。
    于是這里選擇了主席樹。。。。
    二分答案來找區間第K小。。。。
  1 #include<iostream>
  2 #include<cassert>
  3 #include<cstring>
  4 #include<algorithm>
  5 #include<cstdio>
  6 using namespace std;
  7 const int N = 100005;
  8 int num[N],now;
  9 pair<int,int> hash[N];
 10 int seg[N<<2],M;
 11 // chairman tree
 12 const int MAXN = 2000000;
 13 int sum[MAXN],pl[MAXN],pr[MAXN],size,root[N];
 14 int make_tree(int l, int r, int val){
 15     int t = size++;
 16 //    cout<<l<<" "<<r<<" "<<val<<endl;
 17     if(l == r) { sum[t] = val; return t;}
 18     int m = l+r>>1;
 19     pl[t] = make_tree(l,m,val);
 20     pr[t] = make_tree(m+1,r,val);
 21     sum[t] = sum[pl[t]] + sum[pr[t]];
 22     return t;
 23 }
 24 int make_tree(int l,int r, int _pl,int _pr){
 25     int t = size ++;
 26     pl[t] = _pl; pr[t] =_pr;
 27     sum[t] = sum[_pl] + sum[_pr];    
 28     return t;
 29 }
 30 int change(int L, int R,int pos, int at, int v){
 31     if(L == R)
 32         return make_tree(L,R,v);
 33     int m = L+R >> 1;
 34     if(at <=m ) return make_tree(L,R, change(L,m,pl[pos] , at, v),pr[pos]);
 35     return make_tree(L, R , pl[pos] , change(m+1,R,pr[pos] ,at, v));
 36 }
 37 int ask(int l,int r, int pos, int L, int R){
 38     if(l <= L && r >= R) 
 39         return sum[pos];
 40     int mid = L+R >>1 , val = 0;
 41     if(l <= mid) val += ask(l,r,pl[pos], L, mid);
 42     if(r > mid) val += ask(l,r,pr[pos], mid+1,R);
 43     return val;
 44 }
 45 // init()
 46 int n;
 47 void init(){
 48     size = 0;
 49     for(int i=0;i<n;i++)
 50         hash[i] = make_pair(num[i] , i);
 51     sort(hash, hash + n);
 52     root[0] = make_tree(0,n-1,1);
 53     for(int i=1; i<n; i++){
 54         root[i] = change(0,n-1,root[i-1],hash[i-1].second , -1);
 55     }
 56 //    cout<<size<<endl;
 57     for(int i=30; i>=0; i--) if(n < (1<<i)) M = 1<<i;
 58     for(int i=0; i<M*2; i++) seg[i] = 0;
 59 }
 60 // segment tree
 61 void ins(int pos){
 62     for(pos += M; pos; pos >>= 1)
 63         seg[pos] ++;
 64 }
 65 int find_rank(int pos){
 66     int ans = 0;
 67     for(pos += M; pos; pos >>=1)
 68         if(pos &1) ans += seg[pos ^ 1];
 69     return ans +1;
 70     
 71 }
 72 int find_kth(int k){
 73     int ans=1;
 74     while(ans < M){
 75         int l = ans <<1;
 76         int r = ans <<1|1;
 77         assert(seg[ans] >= k);
 78         if(seg[l] >= k) ans = l;
 79         else {
 80             k -= seg[l];
 81             ans = r;
 82         }
 83     }
 84     return hash[ans-M].first;
 85 }
 86 // main
 87 int find(int val){
 88     int l = 0, r= n;
 89     while(l<r){
 90         int mid = l+r >>1;
 91         if(val <= hash[mid].first) r = mid;
 92         else l = mid+1;
 93     }
 94     return l;
 95 }
 96 int chk(int pos, int l, int r){
 97     int val = ask(l,r,root[pos],0,n-1);
 98     int ans = (r-l+1-val >> 1) + 1;
 99 //    cout<<(ans)<<endl;
100     return ans;
101 }
102 struct Query{
103     char cmd[20];
104     int l,r,v;
105 } query[N*5];
106 int main(){
107     int _test = 0,Q;
108     while(~scanf("%d",&Q)){
109         n = 0;
110         for(int i=0; i<Q;i++){
111             scanf("%s",query[i].cmd);
112             if(!strcmp(query[i].cmd, "Insert")){
113                 scanf("%d",&query[i].v);
114                 num[n++] = query[i].v;
115             }
116             else {
117                 char x = query[i].cmd[strlen(query[i].cmd) -1];
118                 if(x =='1')
119                     scanf("%d%d%d",&query[i].l , &query[i].r ,&query[i].v);
120                 else scanf("%d",&query[i].v);
121             }
122         }
123         init();
124         long long Q1 = 0, Q2= 0, Q3 = 0;
125         printf("Case %d:\n", ++ _test);
126         for(int i=0;i<Q;i++){
127             if(!strcmp(query[i].cmd,"Insert")){
128                 ins(find(query[i].v));
129             }
130             else {
131                 char x = query[i].cmd[strlen(query[i].cmd) -1];
132                 if(x == '1'){
133                     int L = query[i].l-1, R = query[i].r-1 , k = query[i].v;
134                     int l = 0, r = n;
135 //                    cout<<"qeury_1 "<<L<<" "<<R<<" "<<k<<endl;
136                     while(l < r){
137                         int m = l+r >>1;
138 //                        cout<<l<<" "<<r<<" "<<m<<" "<<hash[m].first<<endl;
139                         if(chk(m,L,R) > k) r = m ;
140                         else l = m+1;
141                     }
142                     l--;
143                     int v1 = hash[l].first;
144                     Q1 += v1;
145                     //cout<<v1<<endl;
146                 }
147                 else if(x=='2'){
148                     int u = find(query[i].v);
149                     int v2 = find_rank(u);
150                     //cout<<v2<<endl;
151                     Q2 += v2;
152                 }
153                 else {
154                     int v3 = find_kth(query[i].v);
155                     Q3 += v3;
156                     //cout<<v3<<endl;
157                 }
158             }
159         }
160         cout<<Q1<<endl<<Q2<<endl<<Q3<<endl;
161 //        cout<<size<<endl;
162     }
163 }
164 
posted on 2012-06-21 15:47 西月弦 閱讀(1145) 評論(4)  編輯 收藏 引用 所屬分類: 解題報告

FeedBack:
# re: hdu 3727 主席樹+ 線段樹[未登錄]
2012-06-22 13:00 | David
主席樹跟線段樹有區別嗎?有沒有參考資料?
謝謝。  回復  更多評論
  
# re: hdu 3727 主席樹+ 線段樹
2012-06-22 15:01 | 西月弦
@David
主席樹就是可持久化線段樹,見wjmzbmr的pdf
http://hi.baidu.com/wjbzbmr/blog/item/039405ece58450092cf5341a.html  回復  更多評論
  
# re: hdu 3727 主席樹+ 線段樹[未登錄]
2012-06-24 10:16 | David
@西月弦
謝謝!  回復  更多評論
  
# re: hdu 3727 主席樹+ 線段樹
2012-06-25 02:07 | wuyiqi
果斷去學學  回復  更多評論
  
青青草原综合久久大伊人导航_色综合久久天天综合_日日噜噜夜夜狠狠久久丁香五月_热久久这里只有精品
  • <ins id="pjuwb"></ins>
    <blockquote id="pjuwb"><pre id="pjuwb"></pre></blockquote>
    <noscript id="pjuwb"></noscript>
          <sup id="pjuwb"><pre id="pjuwb"></pre></sup>
            <dd id="pjuwb"></dd>
            <abbr id="pjuwb"></abbr>
            欧美成人三级在线| 欧美专区亚洲专区| 亚洲高清av| 牛牛国产精品| 欧美激情视频一区二区三区免费| 欧美黄色成人网| 亚洲一区二区在| 欧美一区二区三区四区在线| 麻豆精品在线播放| 影音先锋成人资源站| 亚洲婷婷免费| 亚洲在线视频观看| 国内精品久久久久久久影视麻豆 | 亚洲激情在线激情| 欧美在线免费看| 亚洲一区精品视频| 一区二区视频欧美| 亚洲天堂av在线免费| 午夜精品久久久久久久男人的天堂 | 久久av一区二区三区| 欧美二区在线播放| 久久久一本精品99久久精品66| 国产精品视频成人| 国产欧美视频在线观看| 欧美一区精品| 国产一区99| 亚洲高清123| 麻豆成人综合网| 日韩系列欧美系列| 美女被久久久| 狠狠色丁香婷婷综合影院| 亚洲亚洲精品三区日韩精品在线视频| 午夜视频一区在线观看| 久久久久国产精品一区| 亚洲欧美日韩视频一区| 久久人体大胆视频| 久久综合给合久久狠狠色| 国产亚洲精品aa午夜观看| 亚洲大胆人体在线| 久久精品女人天堂| 黑人操亚洲美女惩罚| 亚洲综合国产| 蜜臀久久99精品久久久久久9| 久久国产福利| 亚洲国内高清视频| 日韩一区二区精品在线观看| 国产精品视频免费| 麻豆精品一区二区av白丝在线| 一区二区日韩欧美| 亚洲第一福利在线观看| 欧美一区永久视频免费观看| 国产精品国产三级国产aⅴ9色| 亚洲日本视频| 欧美第一黄网免费网站| 久久久www免费人成黑人精品 | 午夜精品久久久久久久99热浪潮| 欧美国产欧美亚州国产日韩mv天天看完整| 国内精品免费午夜毛片| 欧美一区亚洲| 午夜激情一区| 国产伦精品一区二区三区视频孕妇 | 久久精品首页| 香蕉久久夜色精品国产| 国产精品视频区| 欧美在线观看视频一区二区| 国产精品久久久久久久久搜平片 | 欧美成人一区二区三区在线观看 | 欧美成人午夜视频| 亚洲人成毛片在线播放| 欧美夫妇交换俱乐部在线观看| 久久精品国产96久久久香蕉| 黄色成人免费观看| 欧美国产综合| 欧美一区二区三区在线观看视频 | 亚洲电影在线| 午夜在线a亚洲v天堂网2018| 亚洲午夜女主播在线直播| 欧美成人精品在线| 欧美高清自拍一区| 狠狠色丁香久久婷婷综合丁香| 午夜日韩电影| 久久精品国产欧美激情| 国产精品一级二级三级| 制服丝袜激情欧洲亚洲| 一区二区三区日韩欧美| 欧美另类综合| 一本一本a久久| 亚洲一区二区高清| 国产精品爱啪在线线免费观看| 亚洲日韩欧美视频| 日韩午夜一区| 欧美日韩精品一区| 99在线视频精品| 亚洲欧美日韩国产成人| 国产精品欧美经典| 亚洲欧美日韩在线播放| 久久九九99视频| 伊人婷婷欧美激情| 你懂的国产精品| 亚洲美女视频网| 亚洲欧美日韩在线一区| 国产精品亚洲激情| 久久av二区| 亚洲高清不卡av| 亚洲性图久久| 国产欧美大片| 夜夜嗨av色综合久久久综合网| 亚洲一级黄色av| 日韩视频在线免费| 久久久成人精品| 欧美一区二区三区精品电影| 欧美国产日韩二区| 久久综合网络一区二区| 国产老肥熟一区二区三区| 亚洲欧洲精品一区二区三区波多野1战4| 国产精品视频你懂的| 亚洲午夜在线视频| 国产精品夜夜嗨| 久久理论片午夜琪琪电影网| 亚洲春色另类小说| 欧美一区二区三区视频免费| 久久久久久亚洲精品中文字幕| 极品少妇一区二区三区| 欧美xxx在线观看| 亚洲视频一起| 久久精品30| 亚洲精品网站在线播放gif| 欧美麻豆久久久久久中文| 亚洲尤物在线| 欧美肥婆在线| 欧美中文字幕视频| 亚洲精品韩国| 你懂的一区二区| 日韩一二三区视频| 欧美激情1区| 亚洲国产成人精品女人久久久 | 黄色一区二区在线观看| 性做久久久久久免费观看欧美 | 久久久av毛片精品| 国产字幕视频一区二区| 欧美在线观看一区| 狂野欧美一区| 亚洲精品久久久蜜桃| 欧美国产日本高清在线| 亚洲人成网站精品片在线观看| 免费不卡在线观看| 亚洲欧美自拍偷拍| 亚洲日韩中文字幕在线播放| 久久久久成人精品免费播放动漫| 一本大道av伊人久久综合| 欧美性做爰猛烈叫床潮| 欧美专区亚洲专区| aa成人免费视频| 国产精品视频你懂的| 欧美国产一区二区| 久久精品水蜜桃av综合天堂| 亚洲精品欧美在线| 久久综合中文字幕| 欧美在线一二三四区| 在线中文字幕日韩| 亚洲黄色成人| 在线看欧美日韩| 国产亚洲制服色| 国产欧美日韩免费看aⅴ视频| 欧美日韩日本视频| 欧美国产一区二区在线观看 | 久久久久久久成人| 一区二区在线观看视频| 久久九九国产精品| 欧美成人按摩| 亚洲免费视频网站| 久久精品中文| 亚洲欧美高清| 麻豆成人精品| 亚洲字幕在线观看| 午夜在线精品偷拍| 欧美一级午夜免费电影| 亚洲国产精品激情在线观看| 一区二区激情小说| 亚洲国产精品美女| 欧美在线视频a| 亚洲线精品一区二区三区八戒| 欧美日韩在线综合| 日韩视频一区| 欧美成年视频| 欧美在线黄色| 香蕉尹人综合在线观看| 欧美一区二区三区另类| 亚洲少妇最新在线视频| 久久中文字幕一区| 欧美一区二区三区在| 久久综合色影院| 中国成人亚色综合网站| 免费在线视频一区| 欧美成人亚洲| 亚洲黄色有码视频| 亚洲免费观看高清完整版在线观看| 99精品久久| 午夜一级久久| 免费日韩av电影|