原文链接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的更多相关文章

  1. BZOJ2821 作诗(Poetize) 【分块】

    BZOJ2821 作诗(Poetize) Description 神犇SJY虐完HEOI之后给傻×LYD出了一题: SHY是T国的公主,平时的一大爱好是作诗. 由于时间紧迫,SHY作完诗之后还要虐OI ...

  2. 【分块】BZOJ2821 作诗(Poetize)

    2821: 作诗(Poetize) Time Limit: 50 Sec  Memory Limit: 128 MBSubmit: 3265  Solved: 951[Submit][Status][ ...

  3. BZOJ2821 作诗(Poetize) 分块

    题意 算法 经验总结 代码 题意 不带修改,查询数列[1,n]中[l,r]内的出现正偶数次的数的个数, 数列中的数 <= 1e5, n <= 1e5, 强制在线 算法 ​ 查询的内容: 区 ...

  4. bzoj2821: 作诗(Poetize)

    分块 分sqrt(n)块 F[i][j]表示块i到块j的答案 s[i][j]表示数字i在前j块内出现了几次 #include <iostream> #include <cstdio& ...

  5. 2018.09.30 bzoj2821: 作诗(Poetize)(分块)

    传送门 分块经典题目. 先将数列分块. 然后预处理出每两个块之间有多少个数出现了正偶数次. 这样查询的时候对于中间的完整块直接用预处理出的数组搞定. 剩下的暴力枚举求解. 代码: #include&l ...

  6. 2821: 作诗(Poetize)

    2821: 作诗(Poetize) Time Limit: 50 Sec  Memory Limit: 128 MBSubmit: 1078  Solved: 348[Submit][Status] ...

  7. 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 ...

  8. BZOJ_2821_作诗(Poetize)_分块

    BZOJ_2821_作诗(Poetize)_分块 Description 神犇SJY虐完HEOI之后给傻×LYD出了一题:SHY是T国的公主,平时的一大爱好是作诗.由于时间紧迫,SHY作完诗 之后还要 ...

  9. 【BZOJ2821】作诗(Poetize) 分块

    Description 神犇SJY虐完HEOI之后给傻×LYD出了一题:SHY是T国的公主,平时的一大爱好是作诗.由于时间紧迫,SHY作完诗之后还要虐OI,于是SHY找来一篇长度为N的文章,阅读M次, ...

随机推荐

  1. Win7开机卡在Windows Update 35%的解决办法

    一台Win7老机器,多年未清理,用DISM++清理后,开机重启一直卡在Windows Update 35%转圈圈数小时,无法进入系统. 强制按关机键,F8进入安全模式依然同样现象. 查阅MSDN后,有 ...

  2. Laravel 5.2控制器--为什么命令行创建控制器的时候会报错?

    今天我遇到过一种情况:就是之前创建的控制器名称想更改(仅仅只是大小写问题而已),然后用命令行重新创建新的名称的控制器. 结果就报错了? [ReflectionException] Class App\ ...

  3. ICO和IPO

    ICO 和 IPOIPO:在经过种子轮,天使轮,x轮之后已经非常成熟的状态下进行的投资,基本上是盈利的,上市的是股份,现金流ICO:早期有一个idea,然后又很多的小散户进行的几千万-几亿的投资,风险 ...

  4. 14)django-模板(计数器)

    模块中for循环自带计数器. 使用场景:表格数据相增加序列号. 问:可以使用表的自增加序列做为序列号,但是这个存在个问题,即表中数据被删除,则会出现序列不连续. 1)每个循环都有6个公共方法,双循环有 ...

  5. HTML之Position用法

    在此,先做声明,本篇仅是摘录自互联网,个人认为这篇文章讲的很不错,附于此地与大家共同欣赏. position的四个属性值: 1.relative2.absolute3.fixed4.static下面分 ...

  6. spring各版本jar包和源码

    spring各版本jar包和源码 spring历史版本源码:https://github.com/spring-projects/spring-framework/tags spring历史jar包和 ...

  7. vue的多选框

  8. verilog-产生axis数据流

    首先这是产生aixs数据流的代码 `timescale 1ps/1ps `default_nettype none module axis_switch_0_example_master #( ) ( ...

  9. python并发编程之多线程1

    一多线程的概念介绍 threading模块介绍 threading模块和multiprocessing模块在使用层面,有很大的相似性. 二.开启多线程的两种方式 1.创建线程的开销比创建进程的开销小, ...

  10. 最短路径之Bellman-Ford算法

    第一行为源点个数,边的个数m 接下来m行为a->b和权值 5 52 3 21 2 -31 5 54 5 23 4 3 Bellman-Ford算法(1): #include<iostrea ...