2019-ACM-ICPC 徐州网络赛 I .query 二维偏序问题

题目链接：https://nanti.jisuanke.com/t/41391
题目大意：给出一个长度为n的全排列，q次询问，每次询问区间[l,r]内有多少对数字(p,q)，满足p%q==0。
数据范围n,q同阶10^5

数据范围保证了序列是一个全排列，那么可以预处理出所有的满足条件的数字对，这样的数字对数不会超过n∗logn

#include<bits/stdc++.h>
#define LL long long
using namespace std;

int n, x, sum[100005];

void add(int p, int y)
{
    while(p<=n){
        sum[p]+=y;p+=(p&-p);
    }
}

int cx(int p)
{
    int ans=0;
    while(p){
        ans+=sum[p];p-=(p&-p);
    }
    return ans;
}

vector<int> G[100005];
vector<pair<int, int> > q[100005];
int pos[100005];
int ans[100005]={0};
int main()
{
    int m;
    memset(sum, 0, sizeof(sum));
    scanf("%d%d",&n,&m);
    for(int i=1; i<=n; i++){
        int x;
        scanf("%d", &x);
        pos[x]=i;
    }
    for(int i=1; i<=n; i++){
        for(int j=i*2; j<=n; j+=i){

            int x=pos[i], y=pos[j];
            if(x<y){
                swap(x, y);
            }
            cout<<x<<" "<<y<<endl;
            G[x].push_back(y);//保存a[x]的所有坐标比a[x]小的倍数和因数
        }
    }

    for(int i=1; i<=m; i++){
        int L, R;
        scanf("%d%d", &L,&R);
        q[R].push_back({L, i});
    }

    for(int i=1; i<=n; i++){
        for(int j=0; j<G[i].size(); j++){
            add(G[i][j], 1);
        }
        for(int j=0; j<q[i].size(); j++){
            ans[q[i][j].second]=cx(i)-cx(q[i][j].first-1);
        }
    }

    for(int i=1; i<=m; i++){
        printf("%d\n", ans[i]);
    }

    return 0;
}