I. query
题意:给你一个长度为n的排列,询问m次区间,问区间内有多少对 i,j 使得 gcd(pi,pj)= min(pi,pj)。
题解:gcd(pi,pj)= min(pi,pj)实际上就是pi,pj成倍数关系。首先可以预处理出所有的成倍数关系的pi,pj对。
大概只有(1 + 1/2 + 1/3 + 1/4 + 1/5 + .....)* 1e5 = 1166750 这么多对,命名为贡献区间。
现在问题就是如何快速求出询问区间包含了多少预处理出的贡献区间。
考虑数据结构,并且使用一种比较巧妙的方法。先将查询保存,然后边插入边查询。
将查询区间与贡献区间混合并排序。
排序方法是,区间左端点从大到小的排序,如果左端点相等,那么贡献区间先排前面,然后才是查询区间。
这时候左端点就没有用了。。。使用线段树维护区间和。
遇到贡献区间时,就将其插入到线段树中。注意此时不是区间更新,而是点更新,更新到其右端点。
遇到查询区间时,查询区间 1 ~ r 区间和。
查询时,之前插入的贡献区间左端点都比查询区间要大,所以只要查询区间内有多少个贡献区间的右端点就行了。
#include<bits/stdc++.h>
using namespace std;
#define mid (l + r) / 2
#define ls o * 2
#define rs o * 2 + 1
const int N = 1e5 + 10;
struct node{
int l, r, id, f;//区间左右端点,查询顺序,区分查询区间贡献区间
bool friend operator < (node c, node d){
if(c.l != d.l)
return c.l > d.l;
return c.f > d.f;
}
}e[N << 6];
int ans[N], s[N * 4], a[N], pos[N], cnt;
int n, m;
void up(int o, int l, int r, int pos){
if(l == r){
s[o]++;
return;
}
if(pos <= mid)
up(ls, l, mid, pos);
else
up(rs, mid + 1, r, pos);
s[o] = s[ls] + s[rs];
}
int qu(int o, int l, int r, int ql, int qr){
if(ql <= l && qr >= r)
return s[o];
int res = 0;
if(ql <= mid)
res += qu(ls, l, mid, ql, qr);
if(qr > mid)
res += qu(rs, mid + 1, r, ql, qr);
return res;
}
void gkd()
{
for(int i = 1; i <= cnt; ++i){
if(e[i].f)
up(1, 1, n, e[i].r);
else
ans[e[i].id] = qu(1, 1, n, 1, e[i].r);
}
for(int i = 1; i <= m; ++i)
printf("%d\n", ans[i]);
}
int main()
{
scanf("%d%d", &n, &m);
for(int i = 1; i <= n; ++i)
scanf("%d", &a[i]), pos[a[i]] = i;
int l, r;
for(int i = 1; i <= m; ++i){//记录查询区间
scanf("%d%d", &l, &r);
e[++cnt] = node{l, r, i, 0};
}
for(int i = 1; i <= n; ++i)//记录贡献区间
for(int j = 2 * a[i]; j <= n; j += a[i])
e[++cnt] = node{min(i, pos[j]), max(i, pos[j]), 0, 1};
sort(e + 1, e + 1 + cnt);//混合
gkd();
}