http://acm.hust.edu.cn/vjudge/problem/14689

三个操作

  1. [a,b]覆盖为0
  2. [a,b]覆盖为1
  3. [a,b]颠倒每项

两个查询

  1. [a,b]间1数量
  2. [a,b]间最长连续1数量

其实是区间覆盖,区间颠倒和区间合并的结合,可以先看这几道题

区间合并:http://www.cnblogs.com/qlky/p/5745065.html

区间更新(两种更新方式):http://www.cnblogs.com/qlky/p/5737676.html

可以分为两项做:

1的数量很简单,用len[]保存每个区间数量即可

[a,b]区间内最长连续其实之前没做过,具体做法是:

在区间合并的基础上,在query时取左子树,右子树以及min(mid-L+1,rsum[ls])+min(R-mid,lsum[rs])的最小值

这里的min(mid-L+1,rsum[ls])+min(R-mid,lsum[rs])其实就是左子树的右端和右子树左端结合,前面那项是为了限定左右的最大值,因为要在选定的区间内找连续值

#include <iostream>
#include <string>
#include <cstring>
#include <cstdlib>
#include <cstdio>
#include <cmath>
#include <algorithm>
#include <stack>
#include <queue>
#include <cctype>
#include <vector>
#include <iterator>
#include <set>
#include <map>
#include <sstream>
using namespace std; #define mem(a,b) memset(a,b,sizeof(a))
#define pf printf
#define sf scanf
#define spf sprintf
#define pb push_back
#define debug printf("!\n")
#define MAXN 100000+5
#define MAX(a,b) a>b?a:b
#define blank pf("\n")
#define LL long long
#define ALL(x) x.begin(),x.end()
#define INS(x) inserter(x,x.begin())
#define pqueue priority_queue
#define INF 0x3f3f3f3f #define ls (rt<<1)
#define rs (rt<<1|1) int n,m; int col[MAXN<<],a[MAXN<<],len[MAXN<<]; int sum[MAXN<<],lsum[MAXN<<],rsum[MAXN<<],mz[MAXN<<],lz[MAXN<<],rz[MAXN<<]; void cg(int &a,int &b)
{
int tmp = a;
a = b;
b = tmp;
} void FXOR(int rt,int k)
{
//pf("xor%d %d %d %d t%d %d %d\n",rt,sum[rt],lsum[rt],rsum[rt],mz[rt],lz[rt],rz[rt]);
//统计1数量
len[rt] = k-len[rt]; //最长连续1
cg(sum[rt],mz[rt]);
cg(lsum[rt],lz[rt]);
cg(rsum[rt],rz[rt]);
//pf("xor%d %d %d %d t%d %d %d\n",rt,sum[rt],lsum[rt],rsum[rt],mz[rt],lz[rt],rz[rt]);
if(a[rt]!= -)
{
a[rt]^=;
}
else col[rt]^=;
} void PushDown(int rt,int l,int r)
{
if(l==r) return;
int mid = (l+r)>>;
if(a[rt]!=-)
{
//统计1数量
a[ls] = a[rt]?a[rt]:;
a[rs] = a[rt]?a[rt]:;
len[rs] = a[rt]?(r-mid):;
len[ls] = a[rt]?(mid-l+):;
//最长连续1
sum[rs] = lsum[rs] = rsum[rs] = a[rt]?(r-mid):;
sum[ls] = lsum[ls] = rsum[ls] = a[rt]?(mid-l+):;
mz[rs] = lz[rs] = rz[rs] = a[rt]?:(r-mid);
mz[ls] = lz[ls] = rz[ls] = a[rt]?:(mid-l+); col[rt] = ;
a[rt] = -;
}
if(col[rt])
{
if(l!=r)
{
FXOR(ls,mid-l+);
FXOR(rs,r-mid);
}
col[rt] = ;
}
} void PushUp(int rt,int k)
{
if(k==) return;
//最长连续1
lsum[rt] = lsum[ls];
lz[rt] = lz[ls];
rsum[rt] = rsum[rs];
rz[rt] = rz[rs];
if(lsum[rt]==k-(k>>)) lsum[rt]+=lsum[rs];
if(rsum[rt]==k>>) rsum[rt]+=rsum[ls];
if(lz[rt]==k-(k>>)) lz[rt]+=lz[rs];
if(rz[rt]==k>>) rz[rt]+=rz[ls];
sum[rt] = max(lsum[rs]+rsum[ls],max(sum[rs],sum[ls]));
mz[rt] = max(lz[rs]+rz[ls],max(mz[rs],mz[ls]));
//统计1数量
len[rt] = len[rs]+len[ls];
} void build(int l,int r,int rt)
{
a[rt] = -;
sum[rt] = lsum[rt] = rsum[rt] = len[rt] = col[rt] = ;
mz[rt] = lz[rt] = rz[rt] = r-l+;
if(l==r) return;
int mid = (l+r)>>;
build(l,mid,ls);
build(mid+,r,rs);
} void update(int val,int L,int R,int l,int r,int rt)
{
if(L<=l && r<=R)
{
if(val==)
{
//统计1数量
len[rt] = a[rt] = col[rt] = ; //统计最长连续1
sum[rt] = rsum[rt] = lsum[rt] = ;
mz[rt] = rz[rt] = lz[rt] = r-l+;
}
else if(val==)
{
//统计1数量
len[rt] = r-l+;
a[rt] = ;
col[rt] = ; //统计最长连续1
sum[rt] = rsum[rt] = lsum[rt] = r-l+;
mz[rt] = rz[rt] = lz[rt] = ;
}
else
{
FXOR(rt,r-l+);
}
return;
}
PushDown(rt,l,r);
int mid = (l+r)>>;
if(L <= mid) update(val,L,R,l,mid,ls);
if(R > mid) update(val,L,R,mid+,r,rs);
PushUp(rt,r-l+);
} LL ans = ; //统计1数量
void query(int L,int R,int l,int r,int rt)
{
//pf("%d %d %d %d\n",L,R,l,r);
if(L <= l && r <= R)
{
ans+=len[rt];
return;
}
PushDown(rt,l,r);
int mid = (l+r)>>;
if(L<=mid) query(L,R,l,mid,ls);
if(R>mid) query(L,R,mid+,r,rs);
} //最长连续1
int query2(int L,int R,int l,int r,int rt)
{
//pf("%d %d %d %d\n",L,R,l,r);
if(L <= l && r <= R)
{
return sum[rt];
}
PushDown(rt,l,r);
int res = ;
int mid = (l+r)>>;
if(L<=mid)
{
int tmp = query2(L,R,l,mid,ls);
res = max(tmp,res);
}
if(R>mid)
{
int tmp = query2(L,R,mid+,r,rs);
res = max(tmp,res);
}
res = max(res,min(mid-L+,rsum[ls])+min(R-mid,lsum[rs]));
return res;
} int main()
{
int n,i,j,t,kase=;
sf("%d",&t);
while(t--)
{
sf("%d%d",&n,&m);
build(,n,);
for(i=;i<=n;i++)
{
int tmp;
sf("%d",&tmp);
update(tmp,i,i,,n,);
}
//for(i=1;i<=18;i++) pf("t%d %d %d %d\n",i,sum[i],lsum[i],rsum[i]);
//for(i=1;i<=18;i++) pf("t%d %d %d %d\n",i,mz[i],lz[i],rz[i]);
for(i=;i<=m;i++)
{
int tmp,c,d;
sf("%d%d%d",&tmp,&c,&d);
if(tmp==) update(,c+,d+,,n,);
else if (tmp==) update(,c+,d+,,n,);
else if(tmp==) update(,c+,d+,,n,);
else if(tmp==)
{
ans = ;
query(c+,d+,,n,);
pf("%I64d\n",ans);
}
else
{
pf("%d\n",query2(c+,d+,,n,));
}
//for(j=1;j<=18;j++) pf("t%d %d %d %d t%d %d\n",j,sum[j],lsum[j],rsum[j],a[j],col[j]);
//for(j=1;j<=18;j++) pf("tt%d %d %d %d\n",j,mz[j],lz[j],rz[j]);
}
}
return ;
}
/*
1
10 1000
0 0 0 1 1 0 1 0 1 1
3 0 0
3 1 1
3 2 2
3 3 3
3 4 4
3 5 5
3 6 6
3 7 7
3 8 8
3 9 9
3 0 0
3 0 1
3 0 2
3 0 3
3 0 4
3 0 5
3 0 6
3 0 7
3 0 8
3 0 9
*/

最新文章

  1. SVN基本指令
  2. 有一个团队协同工具,叫Worktile
  3. LB负载均衡之Nginx-Proxy
  4. javascript 中的事件机制
  5. css-文本超出后显示省略号
  6. eclipse环境下tomcat远程调试方法
  7. 基于TLS的反调试技术
  8. (转载)自己实现spring
  9. C++ Primer 5th 第9章 顺序容器
  10. 玩Linux桌面发现一个最佳的组合配置
  11. 显示器 Linux 性能 18 (一个命令行工具传递)
  12. 判断数字 字母 isDigit(), isalpha()
  13. SQL Server 备份所有数据库代码
  14. 记一次 联想杨天A4600K 安装操作系统遇到的尴尬
  15. [SDOI 2015]约数个数和
  16. centos7 安装 smplayer
  17. 【洛谷P1507 NASA的食物计划】
  18. 机器学习&amp;深度学习基础(目录)
  19. 为什么页面设计宽度要控制在960px
  20. 转载如何实现portlet之间的传递参数

热门文章

  1. luoguP2982 [USACO10FEB]慢下来Slowing down
  2. 简述ref与out区别
  3. opencv学习笔记2
  4. Reviewing notes 1.1 of Analytic geometry
  5. [I/O]一览图
  6. SDUT OJ 数据结构实验之二叉树三:统计叶子数
  7. 扩展jQuery---选中指定索引的文本
  8. [Microsoft] 微软技术平台的Cloud Building平台AppVeyor
  9. Call requires API level 11 (current min is 8)报错
  10. Android 文字转语音(TTS)