BZOJ 3236 AHOI 2013 作业

内存限制:512 MiB 时间限制:10000 ms 标准输入输出
 
 
题目类型:传统 评测方式:文本比较

题目大意:

此时己是凌晨两点,刚刚做了Codeforces的小A掏出了英语试卷。英语作业其实不算多,一个小时刚好可以做完。然后是一个小时可以做完的数学作业,接下来是分别都是一个小时可以做完的化学,物理,语文......小A压力巨大。woc他竟然一个小时搞完语文卷子

这是小A碰见了一道非常恶心的数学题,给定了一个长度为n的数列和若干个询问,每个询问是关于数列的区间表示数列的第l个数到第r个数),首先你要统计该区间内大于等于a,小于等于b的数的个数,其次是所有大于等于a,小于等于b的,且在该区间中出现过的数值的个数。

小A望着那数万的数据规模几乎绝望,只能向大神您求救,请您帮帮他吧。

这么说其实不太清楚,说白了就是一个是可以记录重复的而一个不能

输入格式

第一行n,m

接下来n个数表示数列

接下来m行,每行四个数l,r,a,b

输出格式

输出m行,分别对应每个询问,输出两个数,分别为在l到r这段区间中大小在[a,b]中的数的个数,以及大于等于a,小于等于b的,且在该区间中出现过的数值的个数(具体可以参考样例)。

样例

样例输入

3 4
1 2 2
1 2 1 3
1 2 1 1
1 3 1 3
2 3 2 3

样例输出

2 2
1 1
3 2
2 1

数据范围与提示

N=100000,M=1000000

题解:

我们采用莫队的算法,将询问的区间排序

维护两个权值树状数组,一个是可以重复的树状数组A,一个是不能重复的树状数组B,

我们在删除操作时,对于那个可以有重复的树状数组A,我们无条件删除它

而那个不能重复的树状数组B,我们需要判断能否删去

开一个数组,cnt[i]表示权值i出现的次数

如果我们在删除一个i后cnt[i]变成了0,说明现在区间中没有权值i出现,那么我们就可以删去它

对于加的操作,其实是一样的,对于树状数组A,我们无条件加上它

对于树状数组B,如果当前cnt[i]==0,那么说明i是一个新出现的值,需要在树状数组B中加入i

最后一遍询问[a,b]就可得出答案

#include<iostream>
#include<cstdio>
#include<cstring>
#include<algorithm>
#include<string>
#include<vector>
#include<cmath>
#define MAXN 100005
#define MAXM 1000005
#define ll long long
#define re register
using namespace std;
ll n,m,a[MAXN],block[MAXM],blo,l=1,r=0,ans1[MAXM],ans2[MAXM],cnt[MAXN];
struct node{
ll l,r,a,b,id;
friend inline bool operator < (node a,node b){
return block[a.l]==block[b.l]?a.r<b.r:a.l<b.l;
}
}ask[MAXM];
struct BIT{
ll c[MAXN];
ll lowbit(re ll x){return x&(-x);}
inline void update(re ll pos,re ll data){
while(pos<=n){
c[pos]+=data;
pos+=lowbit(pos);
}
}
inline ll query(re ll pos){
re ll sum=0;
while(pos>0){
sum+=c[pos];
pos-=lowbit(pos);
}
return sum;
}
}tree1,tree2;
inline void add(re ll x){
if(x==0) return ;
tree1.update(x,1);
if(cnt[x]==0) tree2.update(x,1);
cnt[x]++;
}
inline void del(re ll x){
if(x==0) return ;
tree1.update(x,-1);
cnt[x]--;
if(cnt[x]==0) tree2.update(x,-1);
}
signed main(){
scanf("%lld%lld",&n,&m);
blo=(ll)sqrt(n);
for(re ll i=1;i<=n;i++)
scanf("%lld",&a[i]);
for(re ll i=1;i<=m;i++){
scanf("%lld%lld%lld%lld",&ask[i].l,&ask[i].r,&ask[i].a,&ask[i].b);
ask[i].id=i;
block[i]=i/blo+1;
}
sort(ask+1,ask+m+1);
for(re ll i=1;i<=m;i++){
while(l<ask[i].l) del(a[l++]);
while(l>ask[i].l) add(a[--l]);
while(r>ask[i].r) del(a[r--]);
while(r<ask[i].r) add(a[++r]);
ans1[ask[i].id]=tree1.query(ask[i].b)-tree1.query(ask[i].a-1);
ans2[ask[i].id]=tree2.query(ask[i].b)-tree2.query(ask[i].a-1);
}
for(re ll i=1;i<=m;i++)
printf("%lld %lld\n",ans1[i],ans2[i]);
return 0;
}

BZOJ 3236 AHOI 2013 作业 莫队+树状数组的更多相关文章

  1. BZOJ 3236 AHOI 2013 作业 莫队算法

    题目大意:给出一些数,问在一个区间中不同的数值有多少种,和在一个区间中不同的数值有多少个. 思路:因为没有改动,所以就想到了莫队算法.然后我写了5K+的曼哈顿距离最小生成树,然后果断T了.(100s的 ...

  2. BZOJ 3236: [Ahoi2013]作业(莫队+树状数组)

    传送门 解题思路 莫队+树状数组.把求\([a,b]\)搞成前缀和形式,剩下的比较裸吧,用\(cnt\)记一下数字出现次数.时间复杂度\(O(msqrt(n)log(n)\),莫名其妙过了. 代码 # ...

  3. bzoj3236 作业 莫队+树状数组

    莫队+树状数组 #include<cstdio> #include<cstring> #include<iostream> #include<algorith ...

  4. BZOJ3236[Ahoi2013]作业——莫队+树状数组/莫队+分块

    题目描述 输入 输出 样例输入 3 4 1 2 2 1 2 1 3 1 2 1 1 1 3 1 3 2 3 2 3 样例输出 2 2 1 1 3 2 2 1 提示 N=100000,M=1000000 ...

  5. COGS.1822.[AHOI2013]作业(莫队 树状数组/分块)

    题目链接: COGS.BZOJ3236 Upd: 树状数组实现的是单点加 区间求和,采用值域分块可以\(O(1)\)修改\(O(sqrt(n))\)查询.同BZOJ3809. 莫队为\(O(n^{1. ...

  6. bzoj 3289: Mato的文件管理 莫队+树状数组

    3289: Mato的文件管理 Time Limit: 40 Sec  Memory Limit: 128 MB[Submit][Status][Discuss] Description Mato同学 ...

  7. Bzoj 3289: Mato的文件管理 莫队,树状数组,逆序对,离散化,分块

    3289: Mato的文件管理 Time Limit: 40 Sec  Memory Limit: 128 MBSubmit: 1539  Solved: 665[Submit][Status][Di ...

  8. bzoj 3289 : Mato的文件管理 (莫队+树状数组)

    题目链接: https://www.lydsy.com/JudgeOnline/problem.php?id=3289 思路: 求区间最小交换的次数将区间变成一个不降序列其实就是求区间逆序对的数量,这 ...

  9. [AHOI2013]作业 莫队 树状数组

    #include<cmath> #include<cstdio> #include<algorithm> #include<string> #inclu ...

随机推荐

  1. 【JZOJ6275】小L的数列

    description analysis 考虑矩阵乘法 设初始\(m×m\)矩阵上\(i\)行\(j\)列的数字表示该矩阵第\(j\)位上\(f[i]\)的指数 那么一开始表示\(f[1..k]\)的 ...

  2. springboot启动报 A child container failed during start 错误解决过程

    启动结果如下: "C:\Program Files\Java\jdk1.8.0_201\bin\java.exe" -agentlib:jdwp=transport=dt_sock ...

  3. (转)详解HttpURLConnection

    请求响应流程 设置连接参数的方法 setAllowUserInteraction setDoInput setDoOutput setIfModifiedSince setUseCaches setD ...

  4. django2 rest api

    版本django 2.2.3,我用的社区版 1.安装django:pycharm——Settings——Project Interpreter——选择你要的Python解释器版本——点击右侧的加号,输 ...

  5. Vuex听说很难?

    Vuex 是什么? Vuex 是一个专为 Vue.js 应用程序开发的状态管理模式.它采用集中式存储管理应用的所有组件状态,并以相应的规则保证状态以一种可预测的方式发生变化.   什么鬼东西 看完这段 ...

  6. Pycharm新建文档的模板设置

    下图演示的是关于python的文档的模板设置! 这样,以后的每一个新建的python的py文件,开头都会有下图中的两句:解释器路径与编码方式 步聚5的第二行内容打错了,应该是utf: #!/usr/b ...

  7. MATLAB 中自定义函数的使用

    MATLAB在文件内部(在函数内部)定义函数,但文件名以开头函数来命名,与Java中每个文件只能有一个公开类,但在文件内部还是可以定义其他非公开类一个道理. 无参函数 do.m function do ...

  8. 《代码整洁之道》ch5~ch9读书笔记 PB16110698(~3.15) 第二周

    <代码整洁之道>ch5~ch9读书笔记 本周我阅读了本书的第5~9章节,进一步了解整洁代码需要注意的几个方面:格式.对象与数据结构.错误处理.边界测试.单元测试和类的规范.以下我将分别记录 ...

  9. UMP系统功能 资源隔离

  10. 类欧几里德算法(洛谷 P5170

    #include <iostream> #include <cstdio> #include <queue> #include <algorithm> ...