分块:

先预处理,将原序列分成长度为len的许多块,计算从第i块到第j块的答案,(可以做到O(n*n/len))。

每次询问时,将询问的区间分成三部分,:左边,中间,右边,中间是尽量大的一个块区间,其答案已经计算得到,左右两边加起来最多有2*len个元素,暴力计算其对答案的影响。O(q*len*f(n)),f(n)是暴力加入一个元素的代价。

这道题f(n)是log(n)

总的复杂度:f(n) = O( n*n/len + q*len*log(n) ),

当len = n*(q*log(n))-1/2时取最小值.

 /**************************************************************
Problem: 2821
User: idy002
Language: C++
Result: Accepted
Time:24352 ms
Memory:13728 kb
****************************************************************/ #include <cstdio>
#include <cstring>
#include <cmath>
#include <algorithm>
#define maxn 100010
#define maxs 1510
#define isok(a) ((a)>0&&!((a)&1))
using namespace std; int n, c, m;
int a[maxn], f[maxs][maxs];
int lx[maxn], rx[maxn], mno[maxn], stot;
int cnt[maxn], cur_ans;
int vs[maxn], vt[maxn], vp[maxn];
bool met[maxn]; struct Pair {
int w, pos;
Pair(){}
Pair( int w, int pos ) : w(w), pos(pos) {}
};
Pair vod[maxn];
bool cmpwp( const Pair & a, const Pair & b ) {
return a.w<b.w || (a.w==b.w && a.pos<b.pos);
}
bool cmpp( const Pair & a, const Pair & b ) {
return a.pos<b.pos;
} void part() {
int len = (int)(double)sqrt((double)n/(log((double)n)/log()))+;
stot = n/len + (n%len!=);
for( int i=; i<=stot; i++ ) {
lx[i] = rx[i-]+;
rx[i] = rx[i-]+len;
if( i==stot ) rx[i] = n;
for( int j=lx[i]; j<=rx[i]; j++ )
mno[j] = i;
}
}
inline void update( int oldv, int newv ) {
bool od, nw;
od = isok(oldv);
nw = isok(newv);
if( od && !nw ) cur_ans--;
else if( !od && nw ) cur_ans++;
} void prep() {
for( int i=; i<=n; i++ )
vod[i] = Pair(a[i],i);
sort( vod+, vod++n, cmpwp );
for( int i=; i<=n; i++ ) {
if( vod[i].w!=vod[i-].w )
vs[vod[i].w] = vt[vod[i-].w] = i;
}
vt[vod[n].w] = n+;
for( int i=; i<=n; i++ ) vp[i] = vod[i].pos; for( int i=; i<=stot; i++ ) {
for( int j=i; j<=stot; j++ ) {
for( int k=lx[j]; k<=rx[j]; k++ ) {
update( cnt[a[k]], cnt[a[k]]+ );
cnt[a[k]]++;
}
f[i][j] = cur_ans;
}
cur_ans = ;
memset( cnt, , sizeof(cnt) );
}
} int qu_cnt( int w, int l, int r ) {
return upper_bound( vp+vs[w], vp+vt[w], r )
- lower_bound( vp+vs[w], vp+vt[w], l );
}
int query( int l, int r ) {
int ml=mno[l], mr=mno[r];
int rt;
if( ml==mr ) {
for( int j=l; j<=r; j++ ) {
update( cnt[a[j]], cnt[a[j]]+ );
cnt[a[j]]+=;
}
rt = cur_ans;
cur_ans = ;
for( int j=l; j<=r; j++ )
cnt[a[j]]--;
return rt;
}
if( mno[l]==mno[l-] ) ml++;
if( mno[r]==mno[r+] ) mr--;
cur_ans = f[ml][mr];
for( int j=l; j<lx[ml]; j++ ) {
if( met[a[j]] ) continue;
met[a[j]] = true;
int t1 = qu_cnt(a[j],l,r);
int t2 = qu_cnt(a[j],lx[ml],rx[mr]);
update( t2, t1 );
}
for( int j=r; j>rx[mr]; j-- ) {
if( met[a[j]] ) continue;
met[a[j]] = true;
int t1 = qu_cnt(a[j],l,r);
int t2 = qu_cnt(a[j],lx[ml],rx[mr]);
update( t2, t1 );
}
rt = cur_ans; cur_ans = ;
for( int j=l; j<lx[ml]; j++ ) met[a[j]]=false;
for( int j=r; j>rx[mr]; j-- ) met[a[j]]=false;
return rt;
} int main() {
scanf( "%d%d%d", &n, &c, &m );
for( int i=; i<=n; i++ ) scanf( "%d", a+i );
part();
prep();
for( int i=,x=,l,r; i<=m; i++ ) {
scanf( "%d%d", &l, &r );
l = (l+x)%n+;
r = (r+x)%n+;
if( l>r ) swap(l,r);
printf( "%d\n", x=query(l,r) );
}
}

bzoj 2821 分块的更多相关文章

  1. bzoj 2821 分块处理

    大题思路就是分块,将n个数分成sqrt(n)个块,然后 处理出一个w数组,w[i,j]代表第i个块到第j个块的答案 那么对于每组询问l,r如果l,r在同一个块中,直接暴力做就行了 如果不在同一个块中, ...

  2. BZOJ 2821 分块+二分

    题意: N个数,M组询问,每次问[l,r]中有多少个数出现正偶数次. 思路: 把N个数分成sqrt(n)块,预处理d[i][j]表示第i块起点到第j块末尾的答案 枚举起点i,并维护一个数组记录每个数到 ...

  3. [BZOJ 2821] 作诗(Poetize) 【分块】

    题目链接:BZOJ - 2821 题目分析 因为强制在线了,所以无法用莫队..可以使用分块来做. 做法是,将 n 个数分成 n/x 个块,每个块大小为 x .先预处理出 f[i][j] ,表示从第 i ...

  4. [BZOJ 2821] 作诗

    Link: BZOJ 2821 传送门 Solution: 一道类似区间众数的经典分块 由于个数为偶数这样的条件不能支持快速合并 因此要先$O(n*sqrt(n))$预处理出$pre[i][j]$表示 ...

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

  6. BZOJ 2821作诗(Poetize) 分块

    Description 有一个长度为n的序列,序列每个元素的范围[1,c],有m个询问x y,表示区间[x,y]中出现正偶数次的数的种类数. Solution 大力分块解决问题. 把序列分块,f[i] ...

  7. bzoj 2821 作诗 分块

    基本思路和蒲公英一样 还是预处理出每两个块间的答案 询问时暴力跑两边的贡献 #include<cstdio> #include<cstring> #include<ios ...

  8. 【BZOJ 2821】作诗

    [题目链接] https://www.lydsy.com/JudgeOnline/problem.php?id=2821 [算法] 如果不强制在线,显然莫队是可以解决此题的,那么,强制在线怎么办呢? ...

  9. bzoj 2741 分块+可持久化trie

    多个询问l,r,求所有子区间异或和中最大是多少 强制在线 做法: 分块+可持久化trie 1.对于每块的左端点i,预处理出i到任意一个j,()i,j)间所有子区间异或和中最大为多少,复杂度O(\(n\ ...

随机推荐

  1. Coursera在线学习---第五节.Logistic Regression

    一.假设函数与决策边界 二.求解代价函数 这样推导后最后发现,逻辑回归参数更新公式跟线性回归参数更新方式一摸一样. 为什么线性回归采用最小二乘法作为求解代价函数,而逻辑回归却用极大似然估计求解? 解答 ...

  2. 远程工具(SSH Secure)连接Centos出现中文乱码问题的解决办法

    问题原因 使用远程工具进行连接时,如果linux有中文文件或目录,显示时会出现乱码,原因是linux编码是UTF-8,而远程工具默认是当前系统本地编码即GBK.所以解决方案是统一两者编码就OK了,但是 ...

  3. docker之设置开机自启动(二)

    docker的自启动 通过sysv-rc-conf等管理 启动脚本 # docker.service #!/bin/sh sudo systemctl enable docker sudo syste ...

  4. HttpURLConnection传json

    private static String sendToWangTing(DataRow dataRow) throws IOException{ String ip = Configuration. ...

  5. js实现图片上传预览

    <!DOCTYPE html> <html> <head> <meta http-equiv="Content-Type" content ...

  6. URAL题解一

    URAL题解一 URAL 1002 题目描述:一种记住手机号的方法就是将字母与数字对应,如图.这样就可以只记住一些单词,而不用记住数字.给出一个数字串和n个单词,用最少的单词数来代替数字串,输出对应的 ...

  7. Ubuntu或者Ubuntu server重新设置IP地址

    1.打开终端输入: sudo vi /etc/network/interfaces 2.进入编辑页面 改一处,添加5行内容,如下图: 3.修改好后esc    :wq进行保存 4.输入: sudo / ...

  8. python模块之itertools

    在循环对象和函数对象中,我们了解了循环器(iterator)的功能.循环器是对象的容器,包含有多个对象.通过调用循环器的next()方法 (__next__()方法,在Python 3.x中),循环器 ...

  9. linux下查看资源使用情况

    //查看占用内存最多的前K的程序ps aux | sort -k4nr | head -K //查看占用CPU最多的前K的程序 ps aux | sort -k3nr | head -K

  10. Perl 连接Oracle 出现OCI missing的问题及解决

    问题描述 新申请了一个虚拟机操作系统: Win Server 2008, 64位 , 8核, 16G Memory 上 http://www.activestate.com/activeperl 下载 ...