【题目链接】 http://www.lydsy.com/JudgeOnline/problem.php?id=3809

【题目大意】

  给定一个长度为n(1<=n<=100000)的正整数序列s(1<=si<=n),
  对于m(1<=m<=1000000)次询问“l,r,a,b”,
  每次输出sl...sr中,权值∈[a,b]的权值的种类数。

【题解】

  用莫队维护,对于修改操作用树状数组,
  发现复杂度为O(msqrtnlogn),难以接受,
  考虑用分块处理整数序列,修改操作O(1),查询操作O(sqrt(n))
  总复杂度O(msqrt(n)+mlog(m))

【代码】

#include <cstdio>
#include <algorithm>
#include <cmath>
const int N=100001,M=1000001;
using namespace std;
typedef long long LL;
int size[N],ans[M],pos[N],color[N],c[400],n,m,limit,bl[400],br[400];
struct Q{
int l,r,id,a,b;
friend bool operator < (const Q &a,const Q &b){
return pos[a.l]<pos[b.l]||(pos[a.l]==pos[b.l]&&a.r<b.r);
}
}ask[M];
bool cmp(const Q &a,const Q &b){return a.id<b.id;}
void read(int&a){
char ch;while(!((ch=getchar())>='0')&&(ch<='9'));
a=ch-'0';while(((ch=getchar())>='0')&&(ch<='9'))a*=10,a+=ch-'0';
}
int query(int x,int y){
int res=0;
int L=pos[x],R=pos[y];
for(int i=L+1;i<R;i++)res+=c[i];
if(L==R){for(int i=x;i<=y;i++)if(size[i])res++;}
else{
for(int i=x;i<=br[L];i++)if(size[i])res++;
for(int i=bl[R];i<=y;i++)if(size[i])res++;
}return res;
}
void modify(int u,int x){
if(size[color[u]]==1&&x==-1)c[pos[color[u]]]--;
if(size[color[u]]==0&&x==1)c[pos[color[u]]]++;
size[color[u]]+=x;
}
int main(){
read(n);read(m);
limit=(int)sqrt(n+0.5);
for(int i=1;i<=n;i++){read(color[i]);pos[i]=(i-1)/limit+1;}
for(int i=1;i<=n;i++){br[pos[i]]=i;if(!bl[pos[i]])bl[pos[i]]=i;}
for(int i=1;i<=m;i++){read(ask[i].l);read(ask[i].r);read(ask[i].a);read(ask[i].b);ask[i].id=i;}
sort(ask+1,ask+m+1);
for(int i=1,l=1,r=0;i<=m;i++){
while(r<ask[i].r)modify(++r,1);
while(r>ask[i].r)modify(r--,-1);
while(l<ask[i].l)modify(l++,-1);
while(l>ask[i].l)modify(--l,1);
ans[ask[i].id]=query(ask[i].a,ask[i].b);
}for(int i=1;i<=m;i++)printf("%d\n",ans[i]);
return 0;
}

最新文章

  1. 【Summary】ANSYS TRANSIENT ANALYSIS
  2. 覆盖的面积(HDU 1255 线段树)
  3. 用 mCustomScrollbar 滚动条插件实现滚动更新添加数据
  4. PHP redis负载均衡代码
  5. uboot完全手册---14
  6. storm系统架构学习
  7. java实现身份证校验
  8. Oracle 介绍 (未完待续)
  9. MSBuild 教程(2)
  10. Maven属性、profile和资源过滤
  11. Python文件处理之文件指针(四)
  12. Git之路——Git的使用
  13. PHP学习笔记十七【面向对象定义类】
  14. startssl证书firefox支持配置
  15. Dubbo与Zookeeper、SpringMVC整合和使用
  16. 手把手教你全家桶之React(一)
  17. 在Java中使用redisTemplate操作缓存
  18. 微信小程序记账本进度七
  19. Java全栈程序员之09:IDEA+GitHub
  20. LinkedHashMap

热门文章

  1. Join vs merge vs lookup
  2. SourceTree 过期,注册导入许可证
  3. 自动化测试===uiautomator2类似appium
  4. 【bzoj4448】SCOI2015 情报传递
  5. HDU 6119 小小粉丝度度熊 双指针
  6. 【软件设计】UML类图怎么看
  7. Iptables基础整理
  8. python代码实现dijkstra算法
  9. 多路复用I/O模型epoll() 模型 代码实现
  10. 深入理解python多进程编程