【JZOJ4792】【NOIP2016提高A组模拟9.21】整除
题目描述
麦克雷有一个1~n的排列,他想知道对于一些区间,有多少对区间内的数(x,y),满足x能被y整除。
输入
第一行包含2个正整数n,m。表示有n个数,m个询问。
接下来一行包含n个正整数,表示麦克雷有的数列。
接下来m行每行包含2个正整数l,r。表示询问区间[l,r]。
输出
共 m 行,每行一个整数,表示满足条件的对数。
样例输入
10 9
1 2 3 4 5 6 7 8 9 10
1 10
2 9
3 8
4 7
5 6
2 2
9 10
5 10
4 10
样例输出
27
14
8
4
2
1
2
7
9
数据范围
30%:1<=n,m<=100
100%:1<=n,m<=2*10^5,1<=pi<=n
解法
对于一个区间[l,r],它的答案=区间[1,r]中的总贡献-区间[1,r]中[1,l-1]的贡献。
证明:
显然,区间[l,r]的答案=区间[1,r]的答案-区间[1,l]的答案-所有(a[i],a[j])的贡献。(其中i∈[1,l],j∈[1,r])
而又,区间[1,l]的答案+所有(a[i],a[j])的贡献=区间[1,r]中[1,l-1]的贡献;
所以区间[l,r]的答案=区间[1,r]中的总贡献-区间[1,r]中[1,l-1]的贡献。
如果把区间按右端点升序排序,那么使用扫描线即可将问题转化为:
区间[1,r]中,任意区间贡献和。
用树状数组维护即可。
代码
#include<iostream>
#include<stdio.h>
#include<math.h>
#include<string.h>
#include<algorithm>
#define ll long long
#define ln(x,y) ll(log(x)/log(y))
#define sqr(x) ((x)*(x))
using namespace std;
const char* fin="aP3.in";
const char* fout="aP3.out";
const ll inf=0x7fffffff;
const ll maxn=200007,maxt=maxn*5;
ll n,m,i,j,k,l,r,ks,cnt;
ll tmp=0;
ll a[maxn],b[maxn],c[maxn];
ll ans[maxn];
ll hh[maxn];
bool bz[maxn];
struct q{
ll l,r,id;
}qu[maxn];
bool cmp(q a,q b){
return a.r<b.r;
}
void change(ll v,ll v1){
if (!bz[v]) return;
v=b[v];
cnt+=v1;
for (;v<=n;v+=v&-v) c[v]+=v1;
}
ll getsum(ll v){
ll x=0;
for (;v;v-=v&-v) x+=c[v];
return x;
}
ll getsum(ll l,ll r){
return getsum(r)-getsum(l-1);
}
void join(ll v){
ll i,j,k;
bz[a[v]]=true;
for (i=1;i<=(ll)(sqrt(a[v]));i++){
if (a[v]%i==0){
change(a[v]/i,1);
if (i*i!=a[v]) change(i,1);
}
}
for (i=a[v]*2;i<=n;i+=a[v]) change(i,1);
}
int main(){
scanf("%d%d",&n,&m);
ks=sqrt(n);
for (i=1;i<=n;i++){
scanf("%d",&a[i]);
b[a[i]]=i;
}
for (i=1;i<=m;i++){
scanf("%d%d",&qu[i].l,&qu[i].r);
qu[i].id=i;
}
sort(qu+1,qu+m+1,cmp);
r=0;
for (i=1;i<=m;i++){
while (r<qu[i].r) {
join(++r);
}
ans[qu[i].id]=cnt-getsum(qu[i].l-1);
}
for (i=1;i<=m;i++) printf("%lld\n",ans[i]);
return 0;
}
启发
离线的区间问题归结为:
莫队算法
范围
对于已知区间,给其加入元素或删除元素以较低复杂度完成。
则适用莫队算法。
扩展
如果加入元素简单,删除元素难,那么就可以使用只增莫队算法。
参考这一题。
分块
坑待填
区间拆分扫描线
范围
当[1,i]容易计算时,且区间[l,r]可拆分成[1,r]和[1,l-1];
则适用区间拆分扫描线。
参考这一题。
扩展
如果[l,r]拆分后满足ans([l,r])=ans([1,r])-ans([1,l-1])-ans(分属于[1,l-1],[l,r]);
那么是可以将ans([l,r])=ans([1,r])-[1,r]下[1,l-1]的贡献。
参考本题。
【JZOJ4792】【NOIP2016提高A组模拟9.21】整除的更多相关文章
- [jzoj 4722] [NOIP2016提高A组模拟8.21] 跳楼机 解题报告 (spfa+同余)
题目链接: http://172.16.0.132/senior/#main/show/4722 题目: DJL为了避免成为一只咸鱼,来找srwudi学习压代码的技巧.Srwudi的家是一幢h层的摩天 ...
- 【JZOJ4791】【NOIP2016提高A组模拟9.21】矩阵
题目描述 在麦克雷的面前出现了一个有n*m个格子的矩阵,每个格子用"."或"#"表示,"."表示这个格子可以放东西,"#" ...
- 【JZOJ4790】【NOIP2016提高A组模拟9.21】选数问题
题目描述 在麦克雷的面前有N个数,以及一个R*C的矩阵.现在他的任务是从N个数中取出R*C个,并填入这个矩阵中.矩阵每一行的法值为本行最大值与最小值的差,而整个矩阵的法值为每一行的法值的最大值.现在, ...
- JZOJ 4732. 【NOIP2016提高A组模拟8.23】函数
4732. [NOIP2016提高A组模拟8.23]函数 (Standard IO) Time Limits: 1500 ms Memory Limits: 262144 KB Detailed ...
- 【NOIP2016提高A组模拟9.14】数列编辑器
题目 分析 比赛上,没有注意到询问只询问光标前面,于是只打了个暴力. 因为询问只询问光标前面,首先,当光标向后每移动到一个位置,顺便将这个位置的前缀和,和最大前缀和求出来. 总之,模拟 #includ ...
- 【NOIP2016提高A组模拟9.24】总结
第一题纯模拟,结果那个出题人脑子似乎进水了,空间限制开了1G!!! 导致我捉摸了半天为什么空间要开那么大,最后只能得出上面的结论. 第二题是个矩阵快速幂,比赛上我没把递推式求出来,但是根据各种乱搞,得 ...
- 【JZOJ4711】【NOIP2016提高A组模拟8.17】Binary
题目描述 输入 输出 样例输入 6 6 8 9 1 13 9 3 1 4 5 2 6 9 1 3 7 2 7 7 1 6 1 2 11 13 样例输出 45 19 21 数据范围 解法 40%暴力即可 ...
- 【JZOJ4747】【NOIP2016提高A组模拟9.3】被粉碎的线段树
题目描述 输入 第一行包括两个正整数,N ,M ,分别表示线段树的宽以及询问次数. 以下N-1 行以先序遍历(dfs深搜顺序)描述一个小R线段树,每行一个正整数表示当前非叶子节点的 mid,保证每个节 ...
- 【JZOJ4746】【NOIP2016提高A组模拟9.3】树塔狂想曲
题目描述 相信大家都在长训班学过树塔问题,题目很简单求最大化一个三角形数塔从上往下走的路径和.走的规则是:(i,j)号点只能走向(i+1,j)或者(i+1,j+1).如下图是一个数塔,映射到该数塔上行 ...
随机推荐
- HZOI2019 B. 那一天她离我而去 最小环
题目大意:https://www.cnblogs.com/Juve/articles/11219089.html 那一天,我们......行啦,不要帮出题人脑补画面了,我们来正经的题解 我们发现我们可 ...
- 系统性能信息模块psutil
目录 前言 获取系统性能信息 CPU 内存 磁盘 网络信息 其他系统信息 系统进程管理方法 进程信息 popen类 查看系统硬件的小脚本 前言 psutil 是一个跨平台库,能够轻松实现获取系统运行的 ...
- C语言学习笔记 函数式宏
不学C光搞PHP不知道还有这种东西-函数式宏,宏前面学过了Macro,编译器在对代码进行编译时会对宏表达式进行展开替换,这样宏就起到了全局变量的作用,这里函数式宏也是类似,编译器进行编译时按函数表达是 ...
- 深入浅出 Java Concurrency (23): 并发容器 part 8 可阻塞的BlockingQueue (3)[转]
在Set中有一个排序的集合SortedSet,用来保存按照自然顺序排列的对象.Queue中同样引入了一个支持排序的FIFO模型. 并发队列与Queue简介 中介绍了,PriorityQueue和Pri ...
- SpringMVC配置顺序的问题
1:web.xml:web应用一经加载,先来找他 1):指明applicationContext的位置 2):引入spring监听,ContextLoaderListe ...
- fork 与 vfork
fork 函数复制父进程(包括父进程的地址空间)产生子进程 在父进程返回子进程ID,在子进程本身返回0. fork一般有两个用处: 1.网络服务进程等待请求,新请求到来,fork一个子进程处理,父进程 ...
- Leetcode551.Student Attendance Record I学生出勤记录1
给定一个字符串来代表一个学生的出勤纪录,这个纪录仅包含以下三个字符: 'A' : Absent,缺勤 'L' : Late,迟到 'P' : Present,到场 如果一个学生的出勤纪录中不超过一个' ...
- 获取MMSQL数据库表信息
SELECT 表名 then d.name else '' end, 表说明 then isnull(f.value,'') else '' end, 字段序号=a.colorder, 字段名=a.n ...
- mybatis深入理解(二)-----Mybatis数据源与连接池
对于ORM框架而言,数据源的组织是一个非常重要的一部分,这直接影响到框架的性能问题.本文将通过对MyBatis框架的数据源结构进行详尽的分析,并且深入解析MyBatis的连接池.本文首先会讲述MyBa ...
- Spring线程安全的实现机制--ThreadLocal
转载: http://blog.csdn.net/lufeng20/article/details/24314381