LOJ#3094. 「BJOI2019」删数

之前做atcoder做到过这个结论结果我忘了。。。

em,就是\([1,n]\)之间每个数\(i\),然后\([i - cnt[i] + 1,i]\)可以放一条线段,没被线段放的地方就是需要改的数的总和

之后我们线段树维护区间最小值以及个数

我们要注意如果+1后使得一个本来在\([1,N]\)的点越出了范围,那么就要把这个区间给删掉,-1同理,要加进来

值域开成\(N + 2M\)也就是\(4.5*10^{5}\)即可

#include <bits/stdc++.h>
#define fi first
#define se second
#define pii pair<int,int>
#define mp make_pair
#define pb push_back
#define space putchar(' ')
#define enter putchar('\n')
#define eps 1e-10
#define MAXN 500005
#define ba 47
//#define ivorysi
using namespace std;
typedef long long int64;
typedef unsigned int u32;
typedef double db;
template<class T>
void read(T &res) {
res = 0;T f = 1;char c = getchar();
while(c < '0' || c > '9') {
if(c == '-') f = -1;
c = getchar();
}
while(c >= '0' && c <= '9') {
res = res * 10 +c - '0';
c = getchar();
}
res *= f;
}
template<class T>
void out(T x) {
if(x < 0) {x = -x;putchar('-');}
if(x >= 10) {
out(x / 10);
}
putchar('0' + x % 10);
}
const int MAXV = 150000,LEN = 450000;
int N,M;
int a[MAXN],d;
int cnt[LEN + 5];
int getpos(int x) {
return x - d + MAXV;
}
struct node {
int l,r,val,cnt,lz;
}tr[LEN * 4 + 5];
void update(int u) {
tr[u].val = min(tr[u << 1].val,tr[u << 1 | 1].val);
tr[u].cnt = 0;
if(tr[u].val == tr[u << 1].val) tr[u].cnt += tr[u << 1].cnt;
if(tr[u].val == tr[u << 1 | 1].val) tr[u].cnt += tr[u << 1 | 1].cnt;
}
void build(int u,int l,int r) {
tr[u].l = l;tr[u].r = r;
if(l == r) {tr[u].cnt = 1;return;}
int mid = (l + r) >> 1;
build(u << 1,l,mid);
build(u << 1 | 1,mid + 1,r);
update(u);
}
void addlz(int u,int v) {
tr[u].val += v;tr[u].lz += v;
}
void pushdown(int u) {
if(tr[u].lz) {
addlz(u << 1,tr[u].lz);
addlz(u << 1 | 1,tr[u].lz);
tr[u].lz = 0;
}
}
void add(int u,int l,int r,int v) {
if(tr[u].l == l && tr[u].r == r) {
addlz(u,v);
return;
}
pushdown(u);
int mid = (tr[u].l + tr[u].r) >> 1;
if(r <= mid) add(u << 1,l,r,v);
else if(l > mid) add(u << 1 | 1,l,r,v);
else {add(u << 1,l,mid,v);add(u << 1 | 1,mid + 1,r,v);}
update(u);
}
pii Query(int u,int l,int r) {
if(tr[u].l == l && tr[u].r == r) return mp(tr[u].val,tr[u].cnt);
pushdown(u);
int mid = (tr[u].l + tr[u].r) >> 1;
if(r <= mid) return Query(u << 1,l,r);
else if(l > mid) return Query(u << 1 | 1,l,r);
else {
pii a = Query(u << 1,l,mid),b = Query(u << 1 | 1,mid + 1,r);
if(a.fi > b.fi) swap(a,b);
if(a.fi == b.fi) a.se += b.se;
return a;
}
}
void Solve() {
read(N);read(M);
build(1,1,LEN);
for(int i = 1 ; i <= N ; ++i) {
read(a[i]);
a[i] += MAXV;
add(1,a[i] - cnt[a[i]],a[i] - cnt[a[i]],1);
cnt[a[i]]++;
}
int p,x;
for(int i = 1 ; i <= M ; ++i) {
read(p);read(x);
if(p == 0) {
if(x == 1) {
if(cnt[getpos(N)]) {add(1,getpos(N) - cnt[getpos(N)] + 1,getpos(N),-1);}
}
else {
if(cnt[getpos(N + 1)]) {add(1,getpos(N + 1) - cnt[getpos(N + 1)] + 1,getpos(N + 1),1);}
}
d += x;
}
else {
if(a[p] <= getpos(N)) {
add(1,a[p] - cnt[a[p]] + 1,a[p] - cnt[a[p]] + 1,-1);
}
cnt[a[p]]--;
a[p] = x - d + MAXV;
if(a[p] <= getpos(N)) {
add(1,a[p] - cnt[a[p]],a[p] - cnt[a[p]],1);
}
cnt[a[p]]++;
}
pii res = Query(1,getpos(1),getpos(N));
int ans = 0;
if(res.fi == 0) ans = res.se;
out(ans);enter;
}
}
int main() {
#ifdef ivorysi
freopen("f1.in","r",stdin);
#endif
Solve();
}

最新文章

  1. 汽车驱动之家 http://bmw360.cn
  2. 转:安装MySQL遇到MySQL Server Instance Configuration Wizard未响应的解决办法
  3. JS中基本window.document对象操作以及常用事件!
  4. bzoj3135: [Baltic2013]pipesd
  5. MVC4升级MVC5 异常处理
  6. Python开发【第二十篇】:缓存
  7. 为什么c程序里一定要写main函数
  8. Python之路第十二天,高级(6)-paramiko
  9. erlang学习笔记(2)
  10. LINUX下getsockopt和setsockopt函数
  11. IM 融云 之 列表及封装
  12. JAVAEE学习路线分享
  13. 错误:You can&#39;t specify target table &#39;xxx&#39; for update in FROM clause的解决
  14. 【Tyvj 1728】普通平衡树
  15. PAT1002:A+B for Polynomials
  16. 你不知道的JavaScript--Item21 漂移的this
  17. redis 在 php 中的应用(string篇)
  18. Dos 批处理 Shutdown
  19. JavaScript--浅谈DOM操作
  20. hdu1199 线段树

热门文章

  1. Java线程优先级及守护线程(二)
  2. Linux常用目录及目录作用说明
  3. Innodb内存结构
  4. Leetcode题目48.旋转图像(中等)
  5. Final——Nishang
  6. 【转载】详解CI、CD相关概念
  7. JS高级_变量提升和函数提升
  8. Node安装配置
  9. php的时区修改
  10. React Native常用的第三方开源库