计算L到R有多少对,本来想的是用1~R的对数减去1~L-1的对数,结果发现这样算的结果包含了
一个属于1~L-1 另一个属于L~R的合法对。
于是进一步思考,令x = 1~R的对数减去1~L-1的对数,y = 一个属于1~L-1另一个属于L~R的合法对数,答案ans = x - y。
采用离线处理。
i从1到n循环,对所有的query.L=i,减去sum[R] - sum[L-1],即上面所说的y(此时L~R的合法对还没被计算出来)。然后找出所有i的倍数,加到sum里面。再对所有的query.R=i,加上sum[R] - sum[L-1],即上面所说的x(此时L~R的合法对已经计算出来了)。
树状数组,单点更新,区间求值。
由于求的是不同的两个数组成的数对,还需要对每个查询减去自身与自身组成的数对
即ans[i] = - (r[i] - l[i] + 1);
#include <bits/stdc++.h>
using namespace std;
#define N 100500
int a[N],p[N],b[N],l[N],r[N],ans[N];
vector<int> L[N], R[N];
int n,m;
inline int lowbit(int x){
return x & (-x);
}
void update(int x){
while(x <= n) {
b[x]++;
x += lowbit(x);
}
}
int query(int x){
int res = 0;
while(x >= 1) {
res += b[x];
x -= lowbit(x);
}
return res;
}
int main (){
scanf("%d%d", &n, &m);
for(int i = 1; i <= n; ++i) {
scanf("%d", &a[i]);
p[a[i]] = i;
}
for(int i=1;i <= m;++i) {
scanf("%d%d",&l[i],&r[i]);
ans[i] = - (r[i] - l[i] + 1);
L[l[i]].push_back(i);
R[r[i]].push_back(i);
}
for(int i=1;i<=n;++i) {
for(int j=0;j<L[i].size();++j) ans[L[i][j]]-=query(r[L[i][j]])-query(i-1);
for(int j=1;j*a[i]<=n;++j) update(p[j*a[i]]);
for(int j=0;j<R[i].size();++j) ans[R[i][j]]+=query(i)-query(l[R[i][j]]-1);
}
for(int i=1;i<=m;++i) printf("%d\n",ans[i]);
return 0;
}