BZOJ2821 作诗(Poetize) 主席树 bitset
原文链接https://www.lydsy.com/JudgeOnline/problem.php?id=2821
题目传送门 - BZOJ2821
题意
$n$ 个数,$m$ 组询问,每次问 $[l,r]$ 中有多少个数出现正偶数次。
$1\leq n,m,a_i\leq 10^5$
题解
这题的标算是一个分块。但是我不想写分块怎么办?
bitset 大法好!
bitset 大法好!
bitset 大法好!
发现我们可以每隔 $\sqrt n$ 保存一下前缀序列中每一个数字的出现情况。
然后询问的时候相当于找到长度 $L-1$ 的前缀序列的数的出现情况和长度 $R$ 的前缀序列的数的出现情况,异或起来,数 0 的个数即可。这个显然可以通过预处理的东西暴力调整,然后手写 bitset 再加一个预处理可以使统计 0 的个数的复杂度也变成 $O(nm/32)$ 。所以总复杂度 $O(nm/32)$ 。这里注意一下调整 $L-1$ 和 $R$ 的 bitset 不能是同一个!。
然后发现样例萎掉了。
发现一个数没有出现过不能算进去。
那怎么办呢?
补集转化一下:答案变成 “区间不同数值的种数 - 统计异或得到的 bitset 的 1 的个数” 。
显然“区间不同数值的种数”是可以通过转化成二维数点问题的,直接拖主席树板子。
然后就 AC 啦。
时间复杂度 $O(n\log n+n\sqrt n+n^2/32)$ (由于 $n,m,c$ 同阶,所以这里都看作 $n$ )。
代码
#include <bits/stdc++.h>
#define y1 __zzd001
using namespace std;
typedef unsigned uint;
const int N=100405;
int cnt1[65536];
void init(){
for (int i=0;i<65536;i++)
cnt1[i]=cnt1[i>>1]+(i&1);
}
struct BitSet{
uint v[3200];
void clear(){
memset(v,0,sizeof v);
}
uint XOR(int x){v[x>>5]^=1<<(x&31);}
};
int calc(uint v){
return cnt1[v>>16]+cnt1[v&65535];
}
struct Ptree{
static const int S=N*40;
int n;
int root[N],sum[S],ls[S],rs[S],tot;
void build(int &rt,int L,int R){
sum[rt=++tot]=0;
if (L==R)
return;
int mid=(L+R)>>1;
build(ls[rt],L,mid);
build(rs[rt],mid+1,R);
}
void update(int prt,int &rt,int L,int R,int x){
if (!rt||rt==prt)
sum[rt=++tot]=sum[prt];
sum[rt]++;
if (L==R)
return;
if (!ls[rt])
ls[rt]=ls[prt];
if (!rs[rt])
rs[rt]=rs[prt];
int mid=(L+R)>>1;
if (x<=mid)
update(ls[prt],ls[rt],L,mid,x);
else
update(rs[prt],rs[rt],mid+1,R,x);
}
int query(int rt,int L,int R,int xL,int xR){
if (!rt||R<xL||L>xR)
return 0;
if (xL<=L&&R<=xR)
return sum[rt];
int mid=(L+R)>>1;
return query(ls[rt],L,mid,xL,xR)+query(rs[rt],mid+1,R,xL,xR);
}
int Query(int x1,int x2,int y1,int y2){
if (x1>x2||y1>y2)
return 0;
return query(root[x2],0,n,y1,y2)-query(root[x1-1],0,n,y1,y2);
}
void init(int _n){
tot=0;
n=_n;
build(root[0],0,n);
}
void insert(int x,int y){
update(root[x-1],root[x],0,n,y);
}
}pt;
int read(){
int x=0;
char ch=getchar();
while (!isdigit(ch))
ch=getchar();
while (isdigit(ch))
x=(x<<1)+(x<<3)+(ch^48),ch=getchar();
return x;
}
int n,m,a[N],Pre[N],c;
BitSet v1[330],v2[330];
int main(){
init();
n=read(),c=read(),m=read();
pt.init(n);
memset(Pre,0,sizeof Pre);
for (int i=1;i<=n;i++){
a[i]=read();
pt.insert(i,Pre[a[i]]);
Pre[a[i]]=i;
}
v1[0].clear(),v2[0]=v1[0];
for (int i=1;i*320<=n;i++){
v1[i]=v1[i-1];
for (int j=(i-1)*320+1;j<=i*320;j++)
v1[i].XOR(a[j]);
v2[i]=v1[i];
}
int ans=0;
while (m--){
int L=(read()+ans)%n+1,R=(read()+ans)%n+1;
if (L>R)
swap(L,R);
int all=pt.Query(L,R,0,L-1);
L--;
int l=L/320,r=R/320;
BitSet &sl=v1[l],&sr=v2[r];
l*=320,r*=320;
while (l<L)
sl.XOR(a[++l]);
while (r<R)
sr.XOR(a[++r]);
int tot=0;
for (int i=0;i<3200;i++)
tot+=calc(sl.v[i]^sr.v[i]);
printf("%d\n",ans=all-tot);
int ls=L/320*320,rs=R/320*320;
while (l>ls)
sl.XOR(a[l--]);
while (r>rs)
sr.XOR(a[r--]);
}
return 0;
}
BZOJ2821 作诗(Poetize) 主席树 bitset的更多相关文章
- BZOJ2821 作诗(Poetize) 【分块】
BZOJ2821 作诗(Poetize) Description 神犇SJY虐完HEOI之后给傻×LYD出了一题: SHY是T国的公主,平时的一大爱好是作诗. 由于时间紧迫,SHY作完诗之后还要虐OI ...
- 【分块】BZOJ2821 作诗(Poetize)
2821: 作诗(Poetize) Time Limit: 50 Sec Memory Limit: 128 MBSubmit: 3265 Solved: 951[Submit][Status][ ...
- BZOJ2821 作诗(Poetize) 分块
题意 算法 经验总结 代码 题意 不带修改,查询数列[1,n]中[l,r]内的出现正偶数次的数的个数, 数列中的数 <= 1e5, n <= 1e5, 强制在线 算法 查询的内容: 区 ...
- bzoj2821: 作诗(Poetize)
分块 分sqrt(n)块 F[i][j]表示块i到块j的答案 s[i][j]表示数字i在前j块内出现了几次 #include <iostream> #include <cstdio& ...
- 2018.09.30 bzoj2821: 作诗(Poetize)(分块)
传送门 分块经典题目. 先将数列分块. 然后预处理出每两个块之间有多少个数出现了正偶数次. 这样查询的时候对于中间的完整块直接用预处理出的数组搞定. 剩下的暴力枚举求解. 代码: #include&l ...
- 2821: 作诗(Poetize)
2821: 作诗(Poetize) Time Limit: 50 Sec Memory Limit: 128 MBSubmit: 1078 Solved: 348[Submit][Status] ...
- BZOJ 2821: 作诗(Poetize)( 分块 )
分块,分成N^0.5块.O(N^1.5)预处理出sm[i][j]表示前i块中j的出现次数, ans[i][j]表示第i~j块的答案. 然后就可以O(N^0.5)回答询问了.总复杂度O((N+Q)N^0 ...
- BZOJ_2821_作诗(Poetize)_分块
BZOJ_2821_作诗(Poetize)_分块 Description 神犇SJY虐完HEOI之后给傻×LYD出了一题:SHY是T国的公主,平时的一大爱好是作诗.由于时间紧迫,SHY作完诗 之后还要 ...
- 【BZOJ2821】作诗(Poetize) 分块
Description 神犇SJY虐完HEOI之后给傻×LYD出了一题:SHY是T国的公主,平时的一大爱好是作诗.由于时间紧迫,SHY作完诗之后还要虐OI,于是SHY找来一篇长度为N的文章,阅读M次, ...
随机推荐
- LIght OJ 1179
题意: 约瑟夫环问题, 给你N 个人, 没K个出队, 问最后剩下的人的编号. 思路: 直接模拟会T, 对于N个人 , 是一个约瑟夫环问题, 当第一个人出队后, (标号一定为 k % n -1) 剩下的 ...
- 10)django-ORM(创建,字段类型,字段参数)
一:ORM关系对象映射(Object Relational Mapping,简称ORM) ORM分两种: DB first 先在数据库中创建数据库表等 Code first 先创建类,然后根据类创建数 ...
- 彻底搞透OAuth 2.0
OAuth是一个关于授权(authorization)的开放网络标准,在全世界得到广泛应用,目前的版本是2.0版. 本文对OAuth 2.0的设计思路和运行流程,做一个简明通俗的解释,主要参考材料为R ...
- ssh 登录报错 packet_write_wait: Connection to x.x.x.x port 22: Broken pipe
问题 更新个人博客文章时遇到:Error: packet_write_wait: Connection to 192.30.253.113 port 22: Broken pipe packet_wr ...
- Confluence 6 整合到支持的附件存储选项
如果你现在正在存储附件到 WebDav 或者你的数据库中.你可以整合附件的存储到文件系统中.当你的附件从数据库中被合并到文件系统后,你存储在数据库中的附件数据就可以从数据库中删除了. 当附件合并进行的 ...
- HTML5-长按事件
<!DOCTYPE html> <html> <head> <meta charset="utf-8"/> <title> ...
- JPA核心类与使用
点击访问:JPA环境配置(一) Persistence: Persistence用于获取EntityManagerFactory实例,这个类中包含一个名为createEntityManagerFact ...
- 四.awk、sde深度讲解
###sed### 查询 1创建测试文件 cat>person.txt<<EOF> 101,oldboy,CEO> 102,zhangyao,CTO> 103,Al ...
- java Properties (属性集)
加载Properties Properties downloadLog = new Properties(); try { //加载logFile文件 downloadLog.load(new Fil ...
- poj3417lca+树上差分
/* 给定n个点的树,在其中加入m条新边(称为非树边) 现在可以割断一条树边,一条非树边,使图分裂成两个联通块,请问有几种切割方式 对树边进行分情况讨论 如果树边不处在环中,则割断这条树边后可以割断任 ...