[codeforces 301D]：Yaroslav and Divisors

传送门
题目大意：
给定一个[1,n]的排列，有m次询问：
每次给定$Li,Ri$，询问有多少对$q,w\in [Li,Ri]$使得$P_q|P_w$
题解：
这个题太妙辣
大概就是把每一个询问拆成两个，具体题解能够看这里：
https://www.cnblogs.com/oyking/archive/2013/08/14/3256922.html
极限时间复杂度：$O(2472113*logn)$，那个2472113是我跑出来的对数
代码：(404ms)（滑稽）

#include<cstdio>
#include<cstdlib>
#include<cmath>
#include<cstring>
#include<iostream>
#include<algorithm>
#define ll long long
#define max(a,b) a>b?a:b
#define min(a,b) a<b?a:b
using namespace std;
inline int read(){
    int x=0,f=1;char ch=' ';
    while(ch<'0'||ch>'9'){if(ch=='-')f=-1;ch=getchar();}
    while(ch>='0'&&ch<='9')x=(x<<3)+(x<<1)+(ch^48),ch=getchar();
    return x*f;
}
const int N=2e5+5;
struct Q{
    int id,l,r;
}q1[N],q2[N];
inline bool cmp1(const Q& a,const Q& b){return a.l==b.l?a.r<b.r:a.l<b.l;}
inline bool cmp2(const Q& a,const Q& b){return a.r==b.r?a.l<b.l:a.r<b.r;}
int n,m;
int a[N],t[N],ans[N],pos[N];
inline void update(int x){
    while(x<=n){
        t[x]++;
        x+=x&-x;
    }
}
inline int query(int x){
    int ans=0;
    while(x){
        ans+=t[x];
        x-=x&-x;
    }
    return ans;
}
int main(){
    n=read();m=read();
    for(int i=1;i<=n;i++)a[i]=read(),pos[a[i]]=i;
    for(int i=1;i<=m;i++){
        q1[i].l=read(),q1[i].r=read();q1[i].id=i;
        q2[i].l=q1[i].l,q2[i].r=q1[i].r;q2[i].id=i;
    }
    sort(q1+1,q1+m+1,cmp1);
    sort(q2+1,q2+m+1,cmp2);
    for(int i=1,j=1,k=1;i<=n;i++){
        while(j<=m&&q1[j].l==i){
            ans[q1[j].id]-=query(q1[j].r)-query(q1[j].l-1);
            j++;
        }
        for(int x=a[i];x<=n;x+=a[i])update(pos[x]);
        while(k<=m&&q2[k].r==i){
            ans[q2[k].id]+=query(q2[k].r)-query(q2[k].l-1);
            k++;
        }
    }
    for(int i=1;i<=m;i++)printf("%d\n",ans[i]);
    return 0;
}