给你一个1~n的排列A，询问一些区间内有多少对x,y使得A[x]|A[y]
我们考虑将所有的询问按照r排序，让后分成两部分计算：
1.A[i]|A[j] 且i>j
2.A[i]|A[j] 且i<j
这里只考虑第一种，第二种将所有询问对称翻转即可
我们建立一个树状数组，我们令i=1~n
每次将考虑A[i]的倍数，若kA[i]的位置j在i之前，我们就在树状数组中把s[j]+1
那么显然，对于一个询问，若其右端点为i，那么显然这个询问的答案就是树状数组中[l,i]的和
这样就完成了一部分的计算，另一部分可以翻转后按照同样的方法计算
最后答案记得加上r-l+1（区间长度）
#include<stdio.h>
#include<string.h>
#include<algorithm>
#define N 200020
using namespace std;
inline int c(int x){ return x&-x; }
struct E{ int l,r,k; } s[N];
inline bool c1(E a,E b){ return a.r<b.r; }
int n,m,t=0,b[N],d[N],A[N],v[N];
inline void sp(int& a,int& b){ a^=b; b^=a; a^=b; }
inline void add(int x){ for(;x<=n;x+=c(x)) b[x]++; }
inline int sum(int x,int k=0){ for(;x;x^=c(x)) k+=b[x]; return k; }
void KM(){
memset(b,0,sizeof b);
for(int i=1,j;i<=m;++i){
for(j=s[i-1].r+1;j<=s[i].r;++j)
for(int k=2,z=n/d[j];k<=z;++k)
if(v[k*d[j]]<j) add(v[k*d[j]]);
A[s[i].k]+=sum(s[i].r)-sum(s[i].l-1);
}
}
int main(){
scanf("%d%d",&n,&m);
for(int i=1;i<=n;++i) scanf("%d",d+i),v[d[i]]=i;
for(int i=1;i<=m;s[i].k=i,++i) scanf("%d%d",&s[i].l,&s[i].r),A[i]=s[i].r-s[i].l+1;
sort(s+1,s+1+m,c1); KM();
reverse(d+1,d+1+n);
for(int i=1;i<=n;++i) v[i]=n-v[i]+1;
for(int i=1;i<=m;++i) {
sp(s[i].l,s[i].r);
s[i].l=n-s[i].l+1;
s[i].r=n-s[i].r+1;
}
sort(s+1,s+1+m,c1); KM();
for(int i=1;i<=m;++i) printf("%d
",A[i]);
}