【BZOJ3809】Gty的二逼妹子序列

Description

Autumn和Bakser又在研究Gty的妹子序列了!但他们遇到了一个难题。
对于一段妹子们,他们想让你帮忙求出这之内美丽度∈[a,b]的妹子的美丽度的种类数。
为了方便,我们规定妹子们的美丽度全都在[1,n]中。
给定一个长度为n(1<=n<=100000)的正整数序列s(1<=si<=n),对于m(1<=m<=1000000)次询问“l,r,a,b”,每次输出sl...sr中,权值∈[a,b]的权值的种类数

Input

第一行包括两个整数n,m(1<=n<=100000,1<=m<=1000000),表示数列s中的元素数和询问数。
第二行包括n个整数s1...sn(1<=si<=n)。
接下来m行,每行包括4个整数l,r,a,b(1<=l<=r<=n,1<=a<=b<=n),意义见题目描述。
保证涉及的所有数在C++的int内。
保证输入合法。

Output

对每个询问,单独输出一行,表示sl...sr中权值∈[a,b]的权值的种类数。

Sample Input

10 10
4 4 5 1 4 1 5 1 2 1
5 9 1 2
3 4 7 9
4 4 2 5
2 3 4 7
5 10 4 4
3 9 1 1
1 4 5 9
8 9 3 3
2 2 1 6
8 9 1 4

Sample Output

2
0
0
2
1
1
1
0
1
2

HINT

样例的部分解释:
5 9 1 2
子序列为4 1 5 1 2
在[1,2]里的权值有1,1,2,有2种,因此答案为2。
3 4 7 9
子序列为5 1
在[7,9]里的权值有5,有1种,因此答案为1。
4 4 2 5
子序列为1
没有权值在[2,5]中的,因此答案为0。
2 3 4 7
子序列为4 5
权值在[4,7]中的有4,5,因此答案为2。
建议使用输入/输出优化。

题解:一看到题第一感觉仍然是莫队+树状数组,但是时间复杂度O(m*sqrt(n)*log(n)),承受不起啊,但是我们可以分块

对于原来的算法,修改时O(m*sqrt(n)*log(n))的,但是查询却是O(m*log(n))的,我们用分块相当于牺牲一点查询的时间,使修改更快一点

言归正传,我们只需要将权值分块,维护每个块内不同权值的种类数以及区间中每个权值的出现次数,然后查询时先查询[a,b]中间的块的种类数,在暴力统计两边的块内的出现次数,于是修改和查询都是O(m*sqrt(n))的了

别忘了特判a,b在一个块内的情况

从1开始的分块真的很别扭啊~

#include <cstdio>
#include <cstring>
#include <iostream>
#include <cmath>
#include <algorithm>
using namespace std;
int n,m,siz;
struct node
{
int qa,qb,ql,qr,org;
}q[1000010];
int v[100010],sk[100010],s[100010],ans[1000010];
int rd()
{
int ret=0; char gc=getchar();
while(gc<'0'||gc>'9') gc=getchar();
while(gc>='0'&&gc<='9') ret=ret*10+gc-'0',gc=getchar();
return ret;
}
bool cmp(node a,node b)
{
if((a.ql-1)/siz==(b.ql-1)/siz) return a.qr<b.qr;
return (a.ql-1)/siz<(b.ql-1)/siz;
}
int main()
{
n=rd(),m=rd();
siz=(int)sqrt((double)n);
int i,j;
for(i=1;i<=n;i++) v[i]=rd();
for(i=1;i<=m;i++) q[i].ql=rd(),q[i].qr=rd(),q[i].qa=rd(),q[i].qb=rd(),q[i].org=i;
sort(q+1,q+m+1,cmp);
int l=1,r=0;
for(i=1;i<=m;i++)
{
while(r<q[i].qr) r++,sk[(v[r]-1)/siz]+=(s[v[r]]==0),s[v[r]]++;
while(r>q[i].qr) s[v[r]]--,sk[(v[r]-1)/siz]-=(s[v[r]]==0),r--;
while(l>q[i].ql) l--,sk[(v[l]-1)/siz]+=(s[v[l]]==0),s[v[l]]++;
while(l<q[i].ql) s[v[l]]--,sk[(v[l]-1)/siz]-=(s[v[l]]==0),l++;
if((q[i].qa-1)/siz==(q[i].qb-1)/siz)
{
for(j=q[i].qa;j<=q[i].qb;j++) ans[q[i].org]+=(s[j]>0);
continue;
}
for(j=q[i].qa;j<=(q[i].qa-1)/siz*siz+siz&&j<=n;j++) ans[q[i].org]+=(s[j]>0);
for(j=(q[i].qb-1)/siz*siz+1;j<=q[i].qb;j++) ans[q[i].org]+=(s[j]>0);
for(j=(q[i].qa-1)/siz+1;j<(q[i].qb-1)/siz;j++) ans[q[i].org]+=sk[j];
}
for(i=1;i<=m;i++) printf("%d\n",ans[i]);
return 0;
}

【BZOJ3236】[Ahoi2013]作业

别的和上题都一样,就是新增一个求[l,r]中数值∈[a,b]的数的个数,这个怎么搞都可以吧~

#include <cstdio>
#include <cstring>
#include <iostream>
#include <cmath>
#include <algorithm>
using namespace std;
int n,m,siz;
struct node
{
int qa,qb,ql,qr,org;
}q[1000010];
int v[100010],sk[100010],sv[100010],s[100010],ans[1000010],sum[1000010];
int rd()
{
int ret=0; char gc=getchar();
while(gc<'0'||gc>'9') gc=getchar();
while(gc>='0'&&gc<='9') ret=ret*10+gc-'0',gc=getchar();
return ret;
}
bool cmp(node a,node b)
{
if((a.ql-1)/siz==(b.ql-1)/siz) return a.qr<b.qr;
return (a.ql-1)/siz<(b.ql-1)/siz;
}
int main()
{
n=rd(),m=rd();
siz=(int)sqrt((double)n);
int i,j;
for(i=1;i<=n;i++) v[i]=rd();
for(i=1;i<=m;i++) q[i].ql=rd(),q[i].qr=rd(),q[i].qa=rd(),q[i].qb=rd(),q[i].org=i;
sort(q+1,q+m+1,cmp);
int l=1,r=0;
for(i=1;i<=m;i++)
{
while(r<q[i].qr) r++,sk[(v[r]-1)/siz]+=(s[v[r]]==0),s[v[r]]++,sv[(v[r]-1)/siz]++;
while(r>q[i].qr) s[v[r]]--,sv[(v[r]-1)/siz]--,sk[(v[r]-1)/siz]-=(s[v[r]]==0),r--;
while(l>q[i].ql) l--,sk[(v[l]-1)/siz]+=(s[v[l]]==0),s[v[l]]++,sv[(v[l]-1)/siz]++;
while(l<q[i].ql) s[v[l]]--,sv[(v[l]-1)/siz]--,sk[(v[l]-1)/siz]-=(s[v[l]]==0),l++;
if((q[i].qa-1)/siz==(q[i].qb-1)/siz)
{
for(j=q[i].qa;j<=q[i].qb;j++) ans[q[i].org]+=(s[j]>0),sum[q[i].org]+=s[j];
continue;
}
for(j=q[i].qa;j<=(q[i].qa-1)/siz*siz+siz&&j<=n;j++) ans[q[i].org]+=(s[j]>0),sum[q[i].org]+=s[j];
for(j=(q[i].qb-1)/siz*siz+1;j<=q[i].qb;j++) ans[q[i].org]+=(s[j]>0),sum[q[i].org]+=s[j];
for(j=(q[i].qa-1)/siz+1;j<(q[i].qb-1)/siz;j++) ans[q[i].org]+=sk[j],sum[q[i].org]+=sv[j];
}
for(i=1;i<=m;i++) printf("%d %d\n",sum[i],ans[i]);
return 0;
}

【BZOJ3809/3236】Gty的二逼妹子序列 [Ahoi2013]作业 莫队算法+分块的更多相关文章

  1. 【bzoj3809/bzoj3236】Gty的二逼妹子序列/[Ahoi2013]作业 莫队算法+分块

    原文地址:http://www.cnblogs.com/GXZlegend/p/6805252.html bzoj3809 题目描述 Autumn和Bakser又在研究Gty的妹子序列了!但他们遇到了 ...

  2. 【BZOJ-3809】Gty的二逼妹子序列 分块 + 莫队算法

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

  3. 【bzoj3809】Gty的二逼妹子序列

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

  4. 【BZOJ3809】Gty的二逼妹子序列 莫队 分块

    题目描述 给你一个长度为\(n\)的数列,还有\(m\)个询问,对于每个询问\((l,r,a,b)\),输出区间\([l,r]\)有多少范围在\([a,b]\)的权值. \(n\leq 100000, ...

  5. 莫队p2 【bzoj3809】Gty的二逼妹子序列

    发现一篇已经够长了...所以就放在这里吧... http://hzwer.com/5749.html ↑依然是看大牛题解过的   袜子那道题太简单了.... 然后被这道题超时卡了一段时间....... ...

  6. BZOJ3809:Gty的二逼妹子序列

    浅谈莫队:https://www.cnblogs.com/AKMer/p/10374756.html 题目传送门:https://lydsy.com/JudgeOnline/problem.php?i ...

  7. [bzoj3809]Gty的二逼妹子序列_莫队_分块

    Gty的二逼妹子序列 bzoj-3809 题目大意:给定一个n个正整数的序列,m次询问.每次询问一个区间$l_i$到$r_i$中,权值在$a_i$到$b_i$之间的数有多少个. 注释:$1\le n\ ...

  8. [bzoj3809]Gty的二逼妹子序列/[bzoj3236][Ahoi2013]作业

    [bzoj3809]Gty的二逼妹子序列/[bzoj3236][Ahoi2013]作业 bzoj   bzoj 题目大意:一个序列,m个询问在$[l,r]$区间的$[x,y]$范围内的数的个数/种类. ...

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

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

随机推荐

  1. sim800L调试问题

    SIM800L默认上电开机,若此时没有把rst和pwk引脚提前设置好,SIM800l会使stm32进入硬件中断(这可能是因为方面电源的原因导致的),同时sim800L开机后需要一定的时间稳定下来,建议 ...

  2. struts2 转发、重定向概述

    转发等参数传递的注解方式: @Action(value = "operatorRoleAction", results = { @Result(name = "view& ...

  3. AutoFac文档5(转载)

    目录 开始 Registering components 控制范围和生命周期 用模块结构化Autofac xml配置 与.net集成 深入理解Autofac 指导 关于 词汇表 扫描 autofac可 ...

  4. 如何使用VIM的列编辑模式 [转]

    如何使用VIM的列编辑模式? * windows 我使用的VIM FOR WINDOWS,一直都听说VIM有列编辑模式,一直没有使用过,试了几次都失败了.今天又因为工作需要,到网上查了一下,经过不断的 ...

  5. &lt;LeetCode OJ&gt; 234. Palindrome Linked List

    Total Accepted: 40445 Total Submissions: 148124 Difficulty: Easy Given a singly linked list, determi ...

  6. 一定要使用-server参数来调试并发程序

    在阅读JCIP的时候想手工测试一下,结果发现总是没有出现书中描述的并发问题 后来我琢磨,以前记得书上说过,在debugging的环境下,JVM是低并发的,一定要在server的环境下测试,让JVM在高 ...

  7. js 静态方法 静态变量 实例方法 实例变量

    1.静态方法的定义 Js代码   var BaseClass = function() {}; // var BaseClass=new Function(); BaseClass.f1 = func ...

  8. hdu1078(记忆化搜索)

    题意:给出n*n的格子,每个各自里面有些食物,问一只老鼠每次走最多k步所能吃到的最多的食物 这道题目,值得我记住它,re了n次,以前写搜索没有注意的一个小地方,导致re这么多次的 ac代码: #inc ...

  9. QT4.8.5 源码编译记录

    今天想将以前的虚拟机的 QT4.8.5 集成到一个虚拟机里面,所以就重新编译了一次 QT4.8.5的源码 走了一点点小弯路,特此记录. 一.交叉编译器,不能直接从原来的虚拟机里面拷贝,必须使用官网的交 ...

  10. ip: either "dev" is duplicate, or "type" is garbage

    Driver installation Nothing has to be done for the APF51Dev, the APF28Dev and the APF6Dev. For the A ...