2019徐州网络赛I

2019徐州网络赛I

题意

给1个长度为n的数组p,然后给出m次区间询问,求满足区间内满足$min(p_i,p_j)=gcd(p_i,p_j)$的对数.

思路

由$min(p_i,p_j)=gcd(p_i,p_j)$可以看出一个数必然是另一个数的倍数,预处理出对于每个$p_i$以内的约数,然后用树状数组维护(二维偏序查询问题).

代码

#include<bits/stdc++.h>
using namespace std;
#define  lowbits(x)  x&(-x)
const int maxn=1e6+10;
const int MAXN=1e5+10;
struct  node{
    int l,r,type;
    bool operator < (const node& a)const {
        return l==a.l?type<a.type:l>a.l;
    }
}node[MAXN*20];
int a[MAXN];
int pos[MAXN];
int sum[MAXN];
void add(int x){
    for(;x<MAXN;x+=lowbits(x)) sum[x]+=1;
}
int query(int x){
    int res=0;
    for(;x;x-=lowbits(x)) res+=sum[x];
    return res;
}
int ans[MAXN];
int main(){
    int n,m;scanf("%d%d",&n,&m);
    for(int i=1;i<=n;i++) scanf("%d",&a[i]),pos[a[i]]=i;
    int tot=0;
    for(int i=1;i<=n;i++){
        for(int j=2;j*a[i]<=n;j++){
            int l=i,r=pos[j*a[i]];
            if(l==0||r==0) continue;
            if(r<l) swap(l,r);
            node[++tot]={l,r,0};
        }
    }
    for(int i=1;i<=m;i++){
        int l,r;scanf("%d%d",&l,&r);
        node[++tot]={l,r,i};
    }
    sort(node+1,node+tot+1);
    for(int i=1;i<=tot;i++){
        if(node[i].type==0) add(node[i].r);
        else ans[node[i].type]=query(node[i].r);
    }
    for(int i=1;i<=m;i++) printf("%d\n",ans[i]);
    // system("pause");
}