2019徐州网络赛I题

题目

已知数列p是1～n的排列，现在给定m组询问，每组询问给出左右端点l和r，求出[l,r]区间中有多少对<i,j>满足 p[i]%p[j]=0 或 p[j]%p[i]=0,要求i<j。

题解

首先维护一个前缀和sum，其中sum[i]维护的是[1,i]中满足条件的对数。
那么对于询问[l,r]，结果就是sum[r]-sum[l-1]。
但是这看上去有些问题，因为存在i在[1,l-1]，j却在[l,r]上的情况。
于是我们用一个树状数组来维护这种情况，在查询完一个左端点l后，将它所有的倍数（除了它自己）全都更新，每次查询结果减去ask®-ask(l-1)。
举个例子，样例为

5 2
1 2 3 4 5
1 5
3 5

第一次查询完成后，树状数组会在[2,5]区间上都+1，表示这些是1的倍数。在去查[3,5]区间时，就会减去这些值。

代码

代码写的有些丑陋，因为之前读错题了套了一个莫队的板子，多了一个f结构体数组，其实用不到。

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

const int maxn = 100005;
typedef long long ll;

template<class T>
void read(T &x) {
    T f=1;x=0;char s=getchar();
    while(s<'0'||s>'9'){if(s=='-')f=-1;s=getchar();}
    while(s>='0'&&s<='9'){x=x*10+s-'0';s=getchar();}
    x *= f;
}

struct node{
    int r,id;
    node(int r,int id):r(r),id(id){}
    bool operator<(const node &b)const {
        return r<b.r;
    }
};
int block;
struct node2{
    int l, r, id;
    bool operator<(const node2 &b) const
    {
        if (l / block == b.l / block)
            return r < b.r;
        return l < b.l;
    }
}f[maxn]; 
int c[maxn],a[maxn],ans2[maxn],n,m,res,pos[maxn],sum[maxn],cnt[maxn],cnt[maxn],ans[maxn];
vector<int> e[maxn];
vector<node> q[maxn];

int lowbit(int x) {return x&-x;}

void update(int x) {
    for(int i=x;i<=n;i+=lowbit(i))
        c[i]++;
}

int ask(int x) {
    int res = 0;
    for(int i=x;i>0;i-=lowbit(i))
        res += c[i];
    return res;
}

int main() {
    read(n), read(m);
    for(int i=1;i<=n;i++) read(a[i]),pos[a[i]]=i;;
    for(int i=1,l,r;i<=m;i++) {
        read(l), read(r);
        f[i].l = l, f[i].r = r;
        f[i].id = i;
        q[l].push_back(node(r,i));
    }
    for(int i=1;i<=n;i++) {
        for(int j=i;j<=n;j+=i) {
            e[pos[i]].push_back(pos[j]);
            e[pos[j]].push_back(pos[i]);
            int temp = max(pos[j],pos[i]);
            cnt[temp]++;
        }
    }
    for(int i=1;i<=n;i++) sort(e[i].begin(),e[i].end());
    for(int i=1;i<=n;i++) sum[i] = sum[i-1] + cnt[i];
    for(int i=1;i<=n;i++) {
        for(int j=0;j<q[i].size();j++) {
            ans[q[i][j].id] = sum[q[i][j].r] - sum[i-1] - ask(q[i][j].r) + ask(i-1) ;
        }
        int temp = upper_bound(e[i].begin(),e[i].end(),i) - e[i].begin();
        for(int j=temp;j<e[i].size();j++) update(e[i][j]);
    }
    for(int i=1;i<=m;i++) printf("%d\n",ans[i]-(f[i].r-f[i].l+1));
    // for(int i=1;i<=m;i++) cout << ans2[i] << ' ';
    return 0;
}