2019icpc徐州ネットワークゲーム_I_query

問題の意味

配列与えられた、繰り返し間隔尋ね\（[L、R] \ ） を満足\（MIN（[I] 、[J]）== GCD（[I]、[J]）を\） 数\（（i、j）は\ ） 番号。

分析

• 実際には、数に必要な間隔の数倍もあります。
• すべての数だけの数倍を持って、シーケンスがすべて配置されているので、\（nlogn \）ヶ月、あなたはいくつかのすべての暴力を見つけ、その後、二つの部分ために、オフラインで質問をすることができますので、フェンウィックツリーの使用は解決しますそれはすることができます。
• フェンウィックツリーは、実際に求めて逆しようとしている（I <J \＆\ \ \＆A [i]の> [J]）を二次元の半順序であり、このタイトルの需要がされる（Lを\ [I] <L [J] \＆\ ＆R [I]> R [J] \） 二次元の半順序の\（L [i]は、Rが [I] \） 、第一種の寸法を求めています反転数を求める方法は、ツリー・アレイのように計算することができます。

コード

#include <bits/stdc++.h>
using namespace std;
const int N=1e5+50;
int n,m,l,r,c[N],a[N],p[N],ans[N];
struct node{
    int id,l,r;
    bool operator<(const node& rhs)const{
        if(r==rhs.r){
            return l<rhs.l;
        }else{
            return r<rhs.r;
        }
    }
};
vector<node> ns,ps;
int lowbit(int x){
    return x&(-x);
}
void add(int i,int x){
    while(i<=n){
        c[i]+=x;
        i+=lowbit(i);
    }
}
int sum(int i){
    int ans=0;
    while(i){
        ans+=c[i];
        i-=lowbit(i);
    }
    return ans;
}
int main(){
    // freopen("in.txt","r",stdin);
    scanf("%d%d",&n,&m);
    for(int i=1;i<=n;i++){
        scanf("%d",&a[i]);
        p[a[i]]=i;
    }
    for(int i=1;i<=n;i++){
        for(int j=i+i;j<=n;j+=i){
            int a=p[i],b=p[j];
            if(a>b){
                swap(a,b);
            }
            ps.push_back({0,a,b});
        }
    }
    for(int i=1;i<=m;i++){
        scanf("%d%d",&l,&r);
        ns.push_back(node{i,l,r});
    }
    sort(ns.begin(),ns.end());
    sort(ps.begin(),ps.end());
    int j=0;
    int sz=ps.size();
    for(int i=0;i<m;i++){
        while(j<sz && ps[j].r<=ns[i].r){
            add(ps[j].l,1);
            j++;
        }
        ans[ns[i].id]=j-sum(ns[i].l-1);
    }
    for(int i=1;i<=m;i++){
        printf("%d\n",ans[i]);
    }
    return 0;
}