Jzoj4792 整除

给你一个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 200020using 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\n",A[i]);}