BZOJ 3236 AHOI 2013 作业 莫队+树状数组
BZOJ 3236 AHOI 2013 作业
题目大意:
此时己是凌晨两点,刚刚做了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 作业 莫队+树状数组的更多相关文章
- BZOJ 3236 AHOI 2013 作业 莫队算法
题目大意:给出一些数,问在一个区间中不同的数值有多少种,和在一个区间中不同的数值有多少个. 思路:因为没有改动,所以就想到了莫队算法.然后我写了5K+的曼哈顿距离最小生成树,然后果断T了.(100s的 ...
- BZOJ 3236: [Ahoi2013]作业(莫队+树状数组)
传送门 解题思路 莫队+树状数组.把求\([a,b]\)搞成前缀和形式,剩下的比较裸吧,用\(cnt\)记一下数字出现次数.时间复杂度\(O(msqrt(n)log(n)\),莫名其妙过了. 代码 # ...
- bzoj3236 作业 莫队+树状数组
莫队+树状数组 #include<cstdio> #include<cstring> #include<iostream> #include<algorith ...
- 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 ...
- COGS.1822.[AHOI2013]作业(莫队 树状数组/分块)
题目链接: COGS.BZOJ3236 Upd: 树状数组实现的是单点加 区间求和,采用值域分块可以\(O(1)\)修改\(O(sqrt(n))\)查询.同BZOJ3809. 莫队为\(O(n^{1. ...
- bzoj 3289: Mato的文件管理 莫队+树状数组
3289: Mato的文件管理 Time Limit: 40 Sec Memory Limit: 128 MB[Submit][Status][Discuss] Description Mato同学 ...
- Bzoj 3289: Mato的文件管理 莫队,树状数组,逆序对,离散化,分块
3289: Mato的文件管理 Time Limit: 40 Sec Memory Limit: 128 MBSubmit: 1539 Solved: 665[Submit][Status][Di ...
- bzoj 3289 : Mato的文件管理 (莫队+树状数组)
题目链接: https://www.lydsy.com/JudgeOnline/problem.php?id=3289 思路: 求区间最小交换的次数将区间变成一个不降序列其实就是求区间逆序对的数量,这 ...
- [AHOI2013]作业 莫队 树状数组
#include<cmath> #include<cstdio> #include<algorithm> #include<string> #inclu ...
随机推荐
- SpringBatch批处理框架
1.前言:本博客是对于刘相SpringBatch批处理框架的学习 1.1.参考网站:https://docs.spring.io/spring-batch/4.2.x/reference/html/i ...
- HTML --- 简单的标签
HTML --- 简单的标签 html概述和基本结构 html概述 HTML是 HyperText Mark-up Language 的首字母简写,意思是超文本标记语言,超文本指的是超链接,标记指的是 ...
- LUOGU P4159 [SCOI2009]迷路(矩阵乘法)
传送门 解题思路 以前bpw讲过的一道题,顺便复习一下矩阵乘法.做法就是拆点,把每个点拆成\(9\)个点,然后挨个连边.之后若\(i\)与\(j\)之间的边长度为\(x\),就让\(i\)的第\(x\ ...
- error C2712: Cannot use __try in functions that require object unwinding
转自VC错误:http://www.vcerror.com/?p=52 问题描述: error C2712: Cannot use __try in functions that require ob ...
- Vue学习笔记——Vue-router
转载:https://blog.csdn.net/guanxiaoyu002/article/details/81116616 第1节:Vue-router入门 .解读router/index.js文 ...
- laravel请求处理管道例子
例子: <?php interface Middleware{ public static function handle (Closure $next);} class VerifyCsrfT ...
- javaSpring学习总结day_02
使用注解注入: 1.用于创建bean对象 @Component: 作用:相当于配置了一个bean标签 位置:类上面 属性:value,含义是bean的id,当不写时,有默认值,默认值是当前类的短名,首 ...
- python基础语法(变量与数据类型)
python基础语法(变量与数据类型) 一.python变量 python中的变量不需要声明.每个变量在使用钱都需要赋值,变量赋值以后,该变量才会被创建 在python中,变量就是变量,它没有类型,我 ...
- Linux 容器 vs 虚拟机——谁更胜一筹
自从Linux上的容器变得流行以来,了解Linux容器和虚拟机之间的区别变得更加棘手.本文将向您提供详细信息,以了解Linux容器和虚拟机之间的差异. Linux容器vs虚拟机 - 应用程序与操作系统 ...
- 2019-8-31-dotnet-特性-DynamicallyInvokable-是用来做什么的
title author date CreateTime categories dotnet 特性 DynamicallyInvokable 是用来做什么的 lindexi 2019-08-31 16 ...