权值分块,离散化非常蛋疼,只能离散化搞……

需要支持操作:删除<=某个值得所有权值==打标记 O(sqrt(n))

码长和我的平衡树差不多……速度快3倍左右。

  1 #include<cstdio>
  2 #include<cmath>
  3 #include<algorithm>
  4 #include<cstring>
  5 using namespace std;
  6 #define N 201001
  7 struct Point{int v,p;}tmp[N];
  8 bool operator < (const Point &a,const Point &b){return a.v<b.v;}
  9 int n,m,Infu,a[N],c[N],leave,
 10 en/*插入的权值数*/,en2/*离散化之后的权值种类数*/,ma[N];
 11 char op[N];
 12 int num[N],l[460],CH[12],r[460],Num,sumv[460],sum,sz,b[N],all;
 13 bool delta[N];
 14 inline void R(int &x){
 15     char c=0;int f=1;
 16     for(;c<'0'||c>'9';c=getchar())if(c=='-')f=-1;
 17     for(x=0;c>='0'&&c<='9';c=getchar())(x*=10)+=(c-'0');
 18     x*=f;
 19 }
 20 inline void P(int x)
 21 {
 22     if(!x){putchar('0');puts("");return;}
 23     if(x<0){putchar('-');x=-x;}Num=0;
 24     while(x>0)CH[++Num]=x%10,x/=10;
 25     while(Num)putchar(CH[Num--]+48);puts("");
 26 }
 27 void makeblock()
 28 {
 29     sz=sqrt(en2); if(!sz) sz=1;
 30     for(sum=1;sum*sz<en2;sum++)
 31       {
 32           l[sum]=r[sum-1]+1;
 33           r[sum]=sum*sz;
 34           for(int i=l[sum];i<=r[sum];i++) num[i]=sum;
 35       }
 36     l[sum]=r[sum-1]+1;
 37     r[sum]=en2;
 38     for(int i=l[sum];i<=r[sum];i++) num[i]=sum;
 39 }
 40 void pushdown(const int &p)
 41 {
 42     if(delta[p])
 43       {
 44           for(int i=l[p];i<=r[p];i++) b[i]=0;
 45           delta[p]=0;
 46       }
 47 }
 48 inline void Insert(const int &x){pushdown(num[x]); b[x]++; sumv[num[x]]++; all++;}
 49 inline void Delete(const int &v)//删除小于等于v的所有权值
 50 {
 51     int used=all; pushdown(num[v]);
 52     for(int i=v;i>=l[num[v]];i--)
 53       {
 54           sumv[num[i]]-=b[i];
 55           all-=b[i];
 56           b[i]=0;
 57       }
 58     for(int i=num[v]-1;i>=1;i--) if(sumv[i])
 59       {
 60           delta[i]=1;
 61           all-=sumv[i];
 62           sumv[i]=0;
 63       } leave+=(used-all);
 64 }
 65 inline int Kth(const int &x)
 66 {
 67     int cnt=0;
 68     for(int i=sum;;i--)
 69       {
 70         cnt+=sumv[i];
 71         if(cnt>=x)
 72           {
 73             cnt-=sumv[i];
 74             for(int j=r[i];;j--)
 75             {cnt+=b[j]; if(cnt>=x) return j;}
 76           }
 77       }
 78 }
 79 int main()
 80 {
 81     R(n); R(m);
 82     for(int i=1;i<=n;i++)
 83       {
 84           op[i]=getchar(); R(a[i]);
 85           if(op[i]=='I')
 86             {
 87                 tmp[++en].v=a[i]-Infu;//为新员工消除之前工资变化的影响
 88                 tmp[en].p=en;
 89             }
 90           else if(op[i]=='A') Infu+=a[i];
 91           else if(op[i]=='S')
 92           {
 93               Infu-=a[i];
 94               tmp[++en].v=m-Infu-1;//每次删除<=m-Infu-1的权值
 95               tmp[en].p=en;
 96           }
 97       }
 98     sort(tmp+1,tmp+en+1);
 99     ma[c[tmp[1].p]=++en2]=tmp[1].v;
100     for(int i=2;i<=en;i++)
101       {
102           if(tmp[i].v!=tmp[i-1].v) en2++;
103           ma[c[tmp[i].p]=en2]=tmp[i].v;
104       } Infu=en=0; makeblock();
105     for(int i=1;i<=n;i++)
106       {
107           if(op[i]=='I') {en++; if(a[i]>=m) Insert(c[en]);}
108           else if(op[i]=='A') Infu+=a[i];
109           else if(op[i]=='S')
110           {
111             Infu-=a[i];
112             Delete(c[++en]);
113           }
114           else
115             {
116                 if(a[i]>all) puts("-1");
117                 else P(ma[Kth(a[i])]+Infu);
118             }
119       } P(leave);
120     return 0;
121 }

转载于:https://www.cnblogs.com/autsky-jadek/p/4107496.html

【权值分块】bzoj1503 [NOI2004]郁闷的出纳员相关推荐

  1. [BZOJ1503][NOI2004]郁闷的出纳员 无旋Treap

    1503: [NOI2004]郁闷的出纳员 Time Limit: 5 Sec  Memory Limit: 64 MB Description OIER公司是一家大型专业化软件公司,有着数以万计的员 ...

  2. bzoj1503: [NOI2004]郁闷的出纳员

    地址:http://www.lydsy.com/JudgeOnline/problem.php?id=1503 题目: 1503: [NOI2004]郁闷的出纳员 Time Limit: 5 Sec  ...

  3. BZOJ1503[NOI2004]郁闷的出纳员——treap

    OIER公司是一家大型专业化软件公司,有着数以万计的员工.作为一名出纳员,我的任务之一便是统计每位员工的工资.这本来是一份不错的工作,但是令人郁闷的是,我们的老板反复无常,经常调整员工的工资.如果他心 ...

  4. 【题解】 bzoj1503: [NOI2004]郁闷的出纳员 (Splay)

    bzoj1503,懒得复制,戳我戳我 Solution: 我知不知道我是那根筋抽了突然来做splay,调了起码\(3h+\),到第二天才改出来(我好菜啊),当做训练调错吧 一个裸的splay,没啥好说 ...

  5. [题解]bzoj1503(NOI2004)郁闷的出纳员

    Description OIER公司是一家大型专业化软件公司,有着数以万计的员工.作为一名出纳员,我的任务之一便是统计每位员工的工资.这本来是一份不错的工作,但是令人郁闷的是,我们的老板反复无常,经常 ...

  6. BZOJ1503 [NOI2004]郁闷的出纳员

    AC通道:http://www.lydsy.com/JudgeOnline/problem.php?id=1503 [我的感受] 这题郁闷了很久,因为discuss里讨论的问题和我都不是一个问题... ...

  7. 【BZOJ1503】郁闷的出纳员

    1503: [NOI2004]郁闷的出纳员 Time Limit: 5 Sec  Memory Limit: 64 MB Submit: 11250  Solved: 3962 [Submit][St ...

  8. 洛谷P1486 [NOI2004] 郁闷的出纳员 题解

    P1486 [NOI2004] 郁闷的出纳员 题目链接:P1486 [NOI2004] 郁闷的出纳员 题意:维护一个数据结构,支持 插入一个大小为 kkk​ 的值,小于下界时不插入 所有元素加上 kk ...

  9. 【BZOJ4129】Haruna’s Breakfast,树上带修莫队+权值分块求mex

    Time:2016.09.08 Author:xiaoyimi 转载注明出处谢谢 思路: 这道题相当于把昨天学的树上莫队和带修莫队融合了一下,顺便加了一个mex(未出现的最小自然数) 那么主要问题就是 ...

最新文章

  1. exists的用法 python_10 个 Python 开发技巧
  2. H.264/AVC视频压缩编码标准的新进展
  3. linux常用网络命令ping和arping
  4. 魔方微店商城系统 v1
  5. 文件操作命令(TYPE)
  6. MAC ideaeclipse快捷键--笔记
  7. yum安装Apache2.4
  8. 我的面试标准:能干活、基础要好、有潜力!
  9. 【转】javascript Object使用Array的方法
  10. 龙之谷linux一键端1网段手机可玩,龙之谷手游一键端怎么用
  11. 外包,你是外包,麻烦你不要偷吃公司零食,注意素质...
  12. 计算机中URL是指什么 ?
  13. ExpandableListView 去掉默认箭头,设置成别的图片
  14. 深度学习训练经验之迭代次数与准确率之惑
  15. 鸿蒙系统 google应用
  16. 发展科技到底有什么用,转NASA专家给一位修女的一封信
  17. 谷歌二次验证器手机里不能使用怎么办?
  18. linux查看和修改文件与目录的权限
  19. docsify部署静态文件服务器,云开发 Docsify 文档部署
  20. 《美国新闻与世界报道》公布2021年全球最佳大学排行榜

热门文章

  1. DNA复制体结构和工作原理首次被揭示
  2. 百余位中外学者探讨神经科技挑战:伦理担忧与监管难题并存
  3. 当科学遇上众包:9个值得关注的前沿科技算力众包平台
  4. 全球机器换人排行榜!这八个国家遥遥领先
  5. 他研究了5000家AI公司,说人工智能应用该这么做!
  6. 为什么当今的企业都需要人工智能战略?
  7. 远看像乱序执行,近看是内存屏障的 BUG 是如何解决的?
  8. 麒麟信安面向场景化创新,赋能openEuler商业验证
  9. GitHub 标星 2.5K+,U^2-Net 跨界肖像画,完美复刻人物细节!
  10. 研发团队来了高颜值的妹子,这结局万万没想到 | 每日趣闻