[ AHOI 2013 ] 作业 & [ BZOJ 3809 ] Gty的二逼妹子序列
\(\\\)
Description
给出一个长为 \(n\) 的数列 \(A\) 和 \(k\),多次询问:
对于一个区间 \([L_i,R_i]\),问区间内有多少个数在 \([a_i,b_i]\) 内,以及这些数共有多少个不同的值。
- \(n\le 10^5,m\le 10^6\)
\(\\\)
Solution
看到最后一行心都凉了......
真的佩服松松松的速度 不知道比我高到哪里去了
卡常卡到想吐(见代码部分吧)`
\(\\\)
首先肯定莫队,然后考虑第一问。
直接离散化之后权值树状数组,每次新加进来一个就在对应权值处 \(+1\) ,删除 \(-1\) 。
对于第二问,我们无法确定当前区间里有多少个是个问题。
于是直接再开一个辅助桶,以及另一个询问用的权值树状数组。
加入时若以前没有(桶为空),则在这一权值处 \(+1\) ,删除时若桶清空成 \(0\) ,则在对应权值处 \(-1\) 。
回答就直接区间减法即可。
注意 lower_bound
和upper_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的二逼妹子序列的更多相关文章
- BZOJ 3809: Gty的二逼妹子序列
3809: Gty的二逼妹子序列 Time Limit: 80 Sec Memory Limit: 28 MBSubmit: 1387 Solved: 400[Submit][Status][Di ...
- Bzoj 3809: Gty的二逼妹子序列 莫队,分块
3809: Gty的二逼妹子序列 Time Limit: 35 Sec Memory Limit: 28 MBSubmit: 868 Solved: 234[Submit][Status][Dis ...
- BZOJ 3809 Gty的二逼妹子序列 莫队算法+分块
Description Autumn和Bakser又在研究Gty的妹子序列了!但他们遇到了一个难题. 对于一段妹子们,他们想让你帮忙求出这之内美丽度∈[a,b]的妹子的美丽度的种类数. 为了方便,我们 ...
- BZOJ 3809: Gty的二逼妹子序列 & 3236: [Ahoi2013]作业 [莫队]
题意: 询问区间权值在$[a,b]$范围内种类数和个数 莫队 权值分块维护种类数和个数$O(1)-O(\sqrt{N})$ #include <iostream> #include < ...
- BZOJ 3809 Gty的二逼妹子序列(莫队+分块)
[题目链接] http://www.lydsy.com/JudgeOnline/problem.php?id=3809 [题目大意] 给定一个长度为n(1<=n<=100000)的正整数序 ...
- bzoj 3809 Gty的二逼妹子序列——莫队+分块
题目:https://www.lydsy.com/JudgeOnline/problem.php?id=3809 容易想到树状数组维护值域.但修改和查询都是 log 太慢. 考虑有 nsqrt(n) ...
- bzoj 3809 Gty的二逼妹子序列 —— 莫队+分块
题目:https://www.lydsy.com/JudgeOnline/problem.php?id=3809 据说一开始应该想到莫队+树状数组,然而我想的却是莫队+权值线段树... 如果用权值线段 ...
- bzoj 3809 Gty的二逼妹子序列(莫队算法,块状链表)
[题意] 回答若干个询问,(l,r,a,b):区间[l,r]内权值在[a,b]的数有多少[种]. [思路] 考虑使用块状链表实现莫队算法中的插入与删除. 因为权值处于1..n之间,所以我们可以建一个基 ...
- BZOJ.3809.Gty的二逼妹子序列(分块 莫队)
题目链接 /* 25832 kb 26964 ms 莫队+树状数组:增加/删除/查询 都是O(logn)的,总时间复杂度O(m*sqrt(n)*logn),卡不过 莫队+分块:这样查询虽然变成了sqr ...
随机推荐
- oracle:程序后台提示Io异常: The Network Adapter could not establish the connection)
今天要用tomcat部署一个系统,windows环境下,oracle和程序在一台机器上. 在配置jndi时,直接写成 jdbc:oracle:thin:@127.0.0.1:1521:orcl.启动程 ...
- wukong引擎源码分析之搜索——docid有序的数组里二分归并求交集,如果用跳表的话,在插入索引时会更快
searcher.Search(types.SearchRequest{Text: "百度中国"}) // 查找满足搜索条件的文档,此函数线程安全 func (engine *En ...
- LA-3716(sort的神用)
题意: 给出两条长度均为n的NDA链A和B,找出一段最长的字串[l,r]使得该区域的突变位置不超过p%; 思路: sum[i]表示[1,i]中不相同的个数,可得表达式(sum[i]-sum[j])/( ...
- 【矩阵---求A的1到N次幂之和】
引例: Matrix Power Series: 题目大意,给定矩阵A,求A^+A^+A^+...A^N. 题解:已知X=a,可以通过以下矩阵求出ans=a^+a^+...a^=矩阵^(n+)后右上格 ...
- 理解Objective-C Runtime(三)消息转发机制
消息转发机制概述 上一篇博客消息传递机制中讲解了Objective-C中对象的「消息传递机制」.本文需要讲解另外一个重要问题:当对象受到无法处理的消息之后会发生什么情况? 显然,若想令类能理解某条消息 ...
- 利用PyCharm的Profile工具进行Python性能分析
Profile:PyCharm提供了性能分析工具Run->Profile,如下图所示.利用Profile工具可以对代码进行性能分析,找出瓶颈所在. 测试:下面以一段测试代码来说明如何使用pych ...
- MongoDB 用户名密码登录
Mongodb enable authentication MongoDB 默认直接连接,无须身份验证,如果当前机器可以公网访问,且不注意Mongodb 端口(默认 27017)的开放状态,那么Mon ...
- MYSQL数据库学习----MYSQL数据类型
一切数据在计算中都是以二进制形式存储,而8位二进制数就表示1个字节. 通常我们说一种数据类型占多少字节,就是说这种数据类型可以表示多少位的二进制数,同时二进制数可以转换为十进制数,进而得到这种数据类型 ...
- bzoj 4289 TAX —— 点边转化
题目:https://www.lydsy.com/JudgeOnline/problem.php?id=4289 把边转化成点,同一个原有点相连的边中,边权小的向大的连差值的边,大的向小的连0的边: ...
- ASP.NET Core:WebAppCoreApi
ylbtech-ASP.NET Core:WebAppCoreApi 1.返回顶部 1. 2. 3. 4. 2. Controllers返回顶部 1.ValuesControlle ...