\(\\\)

Description


给出一个长为 \(n\) 的数列 \(A\) 和 \(k\),多次询问:

对于一个区间 \([L_i,R_i]\),问区间内有多少个数在 \([a_i,b_i]\) 内,以及这些数共有多少个不同的值。

  • \(n\le 10^5,m\le 10^6\)

By wangyisong1996加强数据

\(\\\)

Solution


看到最后一行心都凉了......

真的佩服松松松的速度 不知道比我高到哪里去了

卡常卡到想吐(见代码部分吧)`

\(\\\)

首先肯定莫队,然后考虑第一问。

直接离散化之后权值树状数组,每次新加进来一个就在对应权值处 \(+1\) ,删除 \(-1\) 。

对于第二问,我们无法确定当前区间里有多少个是个问题。

于是直接再开一个辅助桶,以及另一个询问用的权值树状数组。

加入时若以前没有(桶为空),则在这一权值处 \(+1\) ,删除时若桶清空成 \(0\) ,则在对应权值处 \(-1\) 。

回答就直接区间减法即可。

注意 lower_boundupper_bound 的时候可能越界,所以要加上哨兵。

\(\\\)

Code


还是说一下卡常用了点啥吧.....

  • BZOJ 专用的 int 优化

  • 读入 & 输出优化

  • 莫队对询问排序时的奇偶性讨论

  • 很迷的块的大小,实测 \(\frac{N}{\sqrt M}\) 最快

    关于这个粘一个洛谷日报上的证明

    我们设块长度为 \(S\) ,那么对于任意多个在同一块内的询问,挪动的距离就是 \(n\),一共\(\frac{n}{S}\) 个块,移动的总次数就是\(\frac{n^2}{S}\),移动可能跨越块,所以还要加上一个 \(mS\) 的复杂度,总复杂度为 \(O(\frac{n^2}{S}+mS)\) ,我们要让这个值尽量小,\(S\) 取 \(\frac{n}{\sqrt{m}}\) 是最优的,此时复杂度为

    \[O(\frac{n^2}{\frac{n}{\sqrt{m}}}+m(\frac{n}{\sqrt{m}}))=O(n\sqrt{m})
    \]

#include<cmath>
#include<cstdio>
#include<cctype>
#include<cstdlib>
#include<cstring>
#include<iostream>
#include<algorithm>
#define N 100010
#define M 1000010
#define Rg register
#define gc getchar
using namespace std; inline int rd(){
int x=0; bool f=0; char c=gc();
while(!isdigit(c)){if(c=='-')f=1;c=gc();}
while(isdigit(c)){x=(x<<1)+(x<<3)+(c^48);c=gc();}
return f?-x:x;
} inline void print(int x){
Rg int y=10,len=1;
while(x>=y){y=(y<<1)+(y<<3);++len;}
while(len--){y/=10;putchar('0'+x/y);x%=y;}
} int n,m,tot,ans,ans1[M],ans2[M],s[N],bl[N],cnt[N],tmp[N]; struct Q{int l,r,L,R,id;}q[M]; inline bool cmp(Q x,Q y){
if(bl[x.l]!=bl[y.l]) return bl[x.l]<bl[y.l];
return bl[x.l]&1?x.r<y.r:x.r<y.r;
} struct BIT{
int c[N];
inline int lowbit(int x){return x&-x;}
inline void add(int p,int x){
for(;p<=n;p+=lowbit(p)) c[p]+=x;
}
inline int query(int p){
int res=0;
for(;p;p-=lowbit(p)) res+=c[p];
return res;
}
}bitcnt,bitsum; inline void add(int p){
++cnt[s[p]];
bitsum.add(s[p],1);
if(cnt[s[p]]==1) bitcnt.add(s[p],1);
} inline void del(int p){
--cnt[s[p]];
bitsum.add(s[p],-1);
if(!cnt[s[p]]) bitcnt.add(s[p],-1);
} int main(){
n=rd(); m=rd();
int t=n/sqrt(m);
for(Rg int i=1;i<=n;++i){
tmp[i]=s[i]=rd();
bl[i]=i/t+1;
}
sort(tmp+1,tmp+1+n);
for(Rg int i=1;i<=n;++i){
tmp[++tot]=tmp[i];
while(tmp[i+1]==tmp[i]) ++i;
}
tmp[++tot]=2000000000;
for(Rg int i=1;i<=n;++i) s[i]=lower_bound(tmp+1,tmp+1+tot,s[i])-tmp;
for(Rg int i=1;i<=m;++i){
q[i].l=rd(); q[i].r=rd(); q[i].id=i;
q[i].L=lower_bound(tmp+1,tmp+1+tot,rd())-tmp;
q[i].R=upper_bound(tmp+1,tmp+1+tot,rd())-tmp-1;
}
sort(q+1,q+1+m,cmp);
int l=1,r=1;
cnt[s[1]]=1;
bitcnt.add(s[1],1);
bitsum.add(s[1],1);
for(Rg int i=1;i<=m;++i){
if(q[i].L>q[i].R){
ans1[q[i].id]=ans2[q[i].id]=0;
continue;
}
while(l<q[i].l){del(l);++l;}
while(l>q[i].l){--l;add(l);}
while(r>q[i].r){del(r);--r;}
while(r<q[i].r){++r;add(r);}
ans1[q[i].id]=bitsum.query(q[i].R)-bitsum.query(q[i].L-1);
ans2[q[i].id]=bitcnt.query(q[i].R)-bitcnt.query(q[i].L-1);
}
for(Rg int i=1;i<=m;++i){
print(ans1[i]); putchar(' ');
print(ans2[i]); putchar('\n');
}
return 0;
}

[ AHOI 2013 ] 作业 & [ BZOJ 3809 ] Gty的二逼妹子序列的更多相关文章

  1. BZOJ 3809: Gty的二逼妹子序列

    3809: Gty的二逼妹子序列 Time Limit: 80 Sec  Memory Limit: 28 MBSubmit: 1387  Solved: 400[Submit][Status][Di ...

  2. Bzoj 3809: Gty的二逼妹子序列 莫队,分块

    3809: Gty的二逼妹子序列 Time Limit: 35 Sec  Memory Limit: 28 MBSubmit: 868  Solved: 234[Submit][Status][Dis ...

  3. BZOJ 3809 Gty的二逼妹子序列 莫队算法+分块

    Description Autumn和Bakser又在研究Gty的妹子序列了!但他们遇到了一个难题. 对于一段妹子们,他们想让你帮忙求出这之内美丽度∈[a,b]的妹子的美丽度的种类数. 为了方便,我们 ...

  4. BZOJ 3809: Gty的二逼妹子序列 & 3236: [Ahoi2013]作业 [莫队]

    题意: 询问区间权值在$[a,b]$范围内种类数和个数 莫队 权值分块维护种类数和个数$O(1)-O(\sqrt{N})$ #include <iostream> #include < ...

  5. BZOJ 3809 Gty的二逼妹子序列(莫队+分块)

    [题目链接] http://www.lydsy.com/JudgeOnline/problem.php?id=3809 [题目大意] 给定一个长度为n(1<=n<=100000)的正整数序 ...

  6. bzoj 3809 Gty的二逼妹子序列——莫队+分块

    题目:https://www.lydsy.com/JudgeOnline/problem.php?id=3809 容易想到树状数组维护值域.但修改和查询都是 log 太慢. 考虑有 nsqrt(n) ...

  7. bzoj 3809 Gty的二逼妹子序列 —— 莫队+分块

    题目:https://www.lydsy.com/JudgeOnline/problem.php?id=3809 据说一开始应该想到莫队+树状数组,然而我想的却是莫队+权值线段树... 如果用权值线段 ...

  8. bzoj 3809 Gty的二逼妹子序列(莫队算法,块状链表)

    [题意] 回答若干个询问,(l,r,a,b):区间[l,r]内权值在[a,b]的数有多少[种]. [思路] 考虑使用块状链表实现莫队算法中的插入与删除. 因为权值处于1..n之间,所以我们可以建一个基 ...

  9. BZOJ.3809.Gty的二逼妹子序列(分块 莫队)

    题目链接 /* 25832 kb 26964 ms 莫队+树状数组:增加/删除/查询 都是O(logn)的,总时间复杂度O(m*sqrt(n)*logn),卡不过 莫队+分块:这样查询虽然变成了sqr ...

随机推荐

  1. oracle:程序后台提示Io异常: The Network Adapter could not establish the connection)

    今天要用tomcat部署一个系统,windows环境下,oracle和程序在一台机器上. 在配置jndi时,直接写成 jdbc:oracle:thin:@127.0.0.1:1521:orcl.启动程 ...

  2. wukong引擎源码分析之搜索——docid有序的数组里二分归并求交集,如果用跳表的话,在插入索引时会更快

    searcher.Search(types.SearchRequest{Text: "百度中国"}) // 查找满足搜索条件的文档,此函数线程安全 func (engine *En ...

  3. LA-3716(sort的神用)

    题意: 给出两条长度均为n的NDA链A和B,找出一段最长的字串[l,r]使得该区域的突变位置不超过p%; 思路: sum[i]表示[1,i]中不相同的个数,可得表达式(sum[i]-sum[j])/( ...

  4. 【矩阵---求A的1到N次幂之和】

    引例: Matrix Power Series: 题目大意,给定矩阵A,求A^+A^+A^+...A^N. 题解:已知X=a,可以通过以下矩阵求出ans=a^+a^+...a^=矩阵^(n+)后右上格 ...

  5. 理解Objective-C Runtime(三)消息转发机制

    消息转发机制概述 上一篇博客消息传递机制中讲解了Objective-C中对象的「消息传递机制」.本文需要讲解另外一个重要问题:当对象受到无法处理的消息之后会发生什么情况? 显然,若想令类能理解某条消息 ...

  6. 利用PyCharm的Profile工具进行Python性能分析

    Profile:PyCharm提供了性能分析工具Run->Profile,如下图所示.利用Profile工具可以对代码进行性能分析,找出瓶颈所在. 测试:下面以一段测试代码来说明如何使用pych ...

  7. MongoDB 用户名密码登录

    Mongodb enable authentication MongoDB 默认直接连接,无须身份验证,如果当前机器可以公网访问,且不注意Mongodb 端口(默认 27017)的开放状态,那么Mon ...

  8. MYSQL数据库学习----MYSQL数据类型

    一切数据在计算中都是以二进制形式存储,而8位二进制数就表示1个字节. 通常我们说一种数据类型占多少字节,就是说这种数据类型可以表示多少位的二进制数,同时二进制数可以转换为十进制数,进而得到这种数据类型 ...

  9. bzoj 4289 TAX —— 点边转化

    题目:https://www.lydsy.com/JudgeOnline/problem.php?id=4289 把边转化成点,同一个原有点相连的边中,边权小的向大的连差值的边,大的向小的连0的边: ...

  10. ASP.NET Core:WebAppCoreApi

    ylbtech-ASP.NET Core:WebAppCoreApi 1.返回顶部 1. 2. 3.           4. 2. Controllers返回顶部 1.ValuesControlle ...