BZOJ3236[Ahoi2013]作业——莫队+树状数组/莫队+分块
题目描述
输入
输出
样例输入
1 2 2
1 2 1 3
1 2 1 1
1 3 1 3
2 3 2 3
样例输出
1 1
3 2
2 1
提示
N=100000,M=1000000
莫队+树状数组:
先考虑每次询问没有权值区间限制的情况,将询问离线排序,用一个数组记录答案,莫队即可。
但现在每次询问有了查询的权值区间,显然一个数组无法记录答案,我们用树状数组来记录答案。
对于第一问直接用树状数组即可,对于第二问先用数组记录每个权值是否出现过再用树状数组维护即可。
莫队时间复杂度是$O(mlog_{n}\sqrt{n})$,查询时间复杂度是$O(mlog_{n})$。
#include<set>
#include<map>
#include<queue>
#include<stack>
#include<cmath>
#include<cstdio>
#include<vector>
#include<bitset>
#include<cstring>
#include<iostream>
#include<algorithm>
#define ll long long
using namespace std;
struct lty
{
int l,r;
int a,b;
int num;
}q[1000010];
int ans1[1000010];
int ans2[1000010];
int v1[100010];
int v2[100010];
int t[100010];
int s[100010];
int n,m;
int L,R;
int block;
inline bool cmp(lty x,lty y)
{
return (x.l/block)==(y.l/block)?x.r<y.r:x.l<y.l;
}
inline void add1(int x,int val)
{
for(int i=x;i<=100000;i+=i&-i)
{
v1[i]+=val;
}
}
inline void add2(int x,int val)
{
for(int i=x;i<=100000;i+=i&-i)
{
v2[i]+=val;
}
}
inline int ask1(int x)
{
int res=0;
for(int i=x;i;i-=i&-i)
{
res+=v1[i];
}
return res;
}
inline int ask2(int x)
{
int res=0;
for(int i=x;i;i-=i&-i)
{
res+=v2[i];
}
return res;
}
int main()
{
scanf("%d%d",&n,&m);
block=(int)sqrt(n+0.5);
block-=rand()%30;
for(int i=1;i<=n;i++)
{
scanf("%d",&s[i]);
}
for(int i=1;i<=m;i++)
{
scanf("%d%d%d%d",&q[i].l,&q[i].r,&q[i].a,&q[i].b);
q[i].num=i;
}
sort(q+1,q+1+m,cmp);
L=1;
for(int i=1;i<=m;i++)
{
while(L>q[i].l)
{
L--;
if(!t[s[L]])
{
add2(s[L],1);
}
t[s[L]]++;
add1(s[L],1);
}
while(R<q[i].r)
{
R++;
if(!t[s[R]])
{
add2(s[R],1);
}
t[s[R]]++;
add1(s[R],1);
}
while(L<q[i].l)
{
if(t[s[L]]==1)
{
add2(s[L],-1);
}
t[s[L]]--;
add1(s[L],-1);
L++;
}
while(R>q[i].r)
{
if(t[s[R]]==1)
{
add2(s[R],-1);
}
t[s[R]]--;
add1(s[R],-1);
R--;
}
ans1[q[i].num]=ask1(q[i].b)-ask1(q[i].a-1);
ans2[q[i].num]=ask2(q[i].b)-ask2(q[i].a-1);
}
for(int i=1;i<=m;i++)
{
printf("%d %d\n",ans1[i],ans2[i]);
}
}
莫队+分块:
可以发现莫队+树状数组的做法时间复杂度主要在双指针移动时对答案的维护。
我们将树状数组改成分块,那么单次移动指针时间复杂度为$O(1)$,而单次查询的时间复杂度是$O(\sqrt{n})$。
这样莫队的时间复杂度是$O(m\sqrt{n})$,查询的时间复杂度为$O(m\sqrt{n})$。
#include<set>
#include<map>
#include<queue>
#include<stack>
#include<cmath>
#include<cstdio>
#include<vector>
#include<bitset>
#include<cstring>
#include<iostream>
#include<algorithm>
#define ll long long
using namespace std;
struct lty
{
int l,r;
int a,b;
int num;
}q[1000010];
int ans1[1000010];
int ans2[1000010];
int sum1[100010];
int sum2[100010];
int v1[100010];
int v2[100010];
int t[100010];
int s[100010];
int n,m;
int L,R;
int block;
inline bool cmp(lty x,lty y)
{
return (x.l/block)==(y.l/block)?x.r<y.r:x.l<y.l;
}
inline int get(int x)
{
return (x-1)/200+1;
}
inline void add1(int x,int val)
{
sum1[get(x)]+=val;
v1[x]+=val;
}
inline void add2(int x,int val)
{
sum2[get(x)]+=val;
v2[x]+=val;
}
inline int ask1(int l,int r)
{
int res=0;
for(int i=get(l)+1;i<=get(r)-1;i++)
{
res+=sum1[i];
}
if(get(l)==get(r))
{
for(int i=l;i<=r;i++)
{
res+=v1[i];
}
return res;
}
for(int i=l;i<=min(n,get(l)*200);i++)
{
res+=v1[i];
}
for(int i=(get(r)-1)*200+1;i<=r;i++)
{
res+=v1[i];
}
return res;
}
inline int ask2(int l,int r)
{
int res=0;
for(int i=get(l)+1;i<=get(r)-1;i++)
{
res+=sum2[i];
}
if(get(l)==get(r))
{
for(int i=l;i<=r;i++)
{
res+=v2[i];
}
return res;
}
for(int i=l;i<=min(n,get(l)*200);i++)
{
res+=v2[i];
}
for(int i=(get(r)-1)*200+1;i<=r;i++)
{
res+=v2[i];
}
return res;
}
int main()
{
scanf("%d%d",&n,&m);
block=(int)sqrt(n+0.5);
for(int i=1;i<=n;i++)
{
scanf("%d",&s[i]);
}
for(int i=1;i<=m;i++)
{
scanf("%d%d%d%d",&q[i].l,&q[i].r,&q[i].a,&q[i].b);
q[i].num=i;
}
sort(q+1,q+1+m,cmp);
L=1;
for(int i=1;i<=m;i++)
{
while(L>q[i].l)
{
L--;
if(!t[s[L]])
{
add2(s[L],1);
}
t[s[L]]++;
add1(s[L],1);
}
while(R<q[i].r)
{
R++;
if(!t[s[R]])
{
add2(s[R],1);
}
t[s[R]]++;
add1(s[R],1);
}
while(L<q[i].l)
{
if(t[s[L]]==1)
{
add2(s[L],-1);
}
t[s[L]]--;
add1(s[L],-1);
L++;
}
while(R>q[i].r)
{
if(t[s[R]]==1)
{
add2(s[R],-1);
}
t[s[R]]--;
add1(s[R],-1);
R--;
}
ans1[q[i].num]=ask1(q[i].a,q[i].b);
ans2[q[i].num]=ask2(q[i].a,q[i].b);
}
for(int i=1;i<=m;i++)
{
printf("%d %d\n",ans1[i],ans2[i]);
}
}
BZOJ3236[Ahoi2013]作业——莫队+树状数组/莫队+分块的更多相关文章
- COGS.1822.[AHOI2013]作业(莫队 树状数组/分块)
题目链接: COGS.BZOJ3236 Upd: 树状数组实现的是单点加 区间求和,采用值域分块可以\(O(1)\)修改\(O(sqrt(n))\)查询.同BZOJ3809. 莫队为\(O(n^{1. ...
- bzoj3236 作业 莫队+树状数组
莫队+树状数组 #include<cstdio> #include<cstring> #include<iostream> #include<algorith ...
- BZOJ 3236: [Ahoi2013]作业(莫队+树状数组)
传送门 解题思路 莫队+树状数组.把求\([a,b]\)搞成前缀和形式,剩下的比较裸吧,用\(cnt\)记一下数字出现次数.时间复杂度\(O(msqrt(n)log(n)\),莫名其妙过了. 代码 # ...
- BZOJ_3289_Mato的文件管理_莫队+树状数组
BZOJ_3289_Mato的文件管理_莫队+树状数组 Description Mato同学从各路神犇以各种方式(你们懂的)收集了许多资料,这些资料一共有n份,每份有一个大小和一个编号 .为了防止他人 ...
- bzoj 3289: Mato的文件管理 莫队+树状数组
3289: Mato的文件管理 Time Limit: 40 Sec Memory Limit: 128 MB[Submit][Status][Discuss] Description Mato同学 ...
- 51nod 1290 Counting Diff Pairs | 莫队 树状数组
51nod 1290 Counting Diff Pairs | 莫队 树状数组 题面 一个长度为N的正整数数组A,给出一个数K以及Q个查询,每个查询包含2个数l和r,对于每个查询输出从A[i]到A[ ...
- 【BZOJ3460】Jc的宿舍(树上莫队+树状数组)
点此看题面 大致题意: 一棵树,每个节点有一个人,他打水需要\(T_i\)的时间,每次询问两点之间所有人去打水的最小等待时间. 伪·强制在线 这题看似强制在线,但实际上,\(pre\ mod\ 2\) ...
- HihoCoder 1488 : 排队接水(莫队+树状数组)
描述 有n个小朋友需要接水,其中第i个小朋友接水需要ai分钟. 由于水龙头有限,小Hi需要知道如果为第l个到第r个小朋友分配一个水龙头,如何安排他们的接水顺序才能使得他们等待加接水的时间总和最小. 小 ...
- BZOJ 3236 莫队+树状数组
思路: 莫队+树状数组 (据说此题卡常数) yzy写了一天(偷笑) 复杂度有点儿爆炸 O(msqrt(n)logn) //By SiriusRen #include <cmath> #in ...
随机推荐
- Linux journalctl命令
在Systemd出现之前,Linux系统及各应用的日志都是分别管理的,Systemd开始统一管理了所有Unit的启动日志,这样带来的好处就是可以只用一个 journalctl命令,查看所有内核和应用的 ...
- Spring Aop底层原理详解
Spring Aop底层原理详解(来源于csdn:https://blog.csdn.net/baomw)
- Java web的一些面试题
1.Tomcat 的优化经验 答:去掉对 web.xml 的监视,把 jsp 提前编辑成 Servlet. 有富余物理内存的情况,加大 tomcat 使用的 jvm 的内存 2.HTTP 请求的 GE ...
- HeadFirst设计模式读书笔记之工厂模式
1. 简单工厂 1. 你开了一家披萨店,点披萨的方法可能是这样: public Pizza orderPizza(String type) { Pizza pizza; if (type.equals ...
- 关于Java中基类构造器的调用问题
在<Java编程思想>第7章复用类中有这样一段话,值得深思.当子类继承了父类时,就涉及到了基类和导出类(子类)这两个类.从外部来看,导出类就像是一个与基类具有相同接口的新类,或许还会有一些 ...
- mysql基本操作(1)
1.mysql数据库客户端安装 brew install mysql-client 2.mysql 连接数据库 mysql -h <数据库地址> -P <端口> -u < ...
- FIFO深度计算
本文设计思想采用明德扬至简设计法.在使用FPGA设计系统时,常需要利用FIFO进行数据缓存保证数据不丢失,因此计算FIFO深度是至关重要的.FIFO的深度主要取决于“最恶劣”的情况,以下对于两种最常见 ...
- About A Scam
事件起因 本篇记录一个我遇到一个诈骗故事. 这两年我陆续有收到邮件,内容为有一大笔遗产我可以继承,让我提供银行卡号,身份证号相关信息. 后面邮件的内容就变为,有一笔公益款项,可以用我名义去管理,让我提 ...
- .Net Core 在Linux服务器下部署程序--(1). Windows 连接 Linux服务器
下载Linux服务器连接软件,市面上有Putty,FinalShell等,我以FinalShell为例,下载地址为 :http://www.hostbuf.com/t/988.html,软件安装结束后 ...
- java加载properties文件的六中基本方式实现
java加载properties文件的方式主要分为两大类:一种是通过import java.util.Properties类中的load(InputStream in)方法加载: 另一种是通过impo ...