19 icpc徐州网络赛 query //树状数组离线+素数筛

query

我没有好好善待树状数组啊,树状数组在考虑离线+区域求和的时候蛮好用

题意:

给一个排列 p p p,给 q q q个询问,每次询问 [ l , r ] [l,r] [l,r]里满足 g c d ( i , j ) = = m i n ( i , j ) ( i < j ) gcd(i,j)==min(i,j)(i<j) gcd(i,j)==min(i,j)(i<j) p a i r pair pair的数量。

思路:

将询问离线,从左到右遍历数组 p p p,遍历到 p [ i ] p[i] p[i],统计新出现的 p a i r pair pair,左端点加入到树状数组。
然后处理询问右端点为 i i i的询问。
预处理每一个 p a i r pair pair,不超过 N l o g N NlogN NlogN个。

/*   Author : Rshs
 *   Data : 2019-09-18-19.28
 */
#include<bits/stdc++.h>
using namespace std;
#define FI first
#define SE second
#define LL long long
#define MP make_pair
#define PII pair<int,int>
#define SZ(a) (int)a.size()
const double pai = acos(-1);
const double eps = 1e-10;
const LL mod = 1e9+7;
const int MX = 1e6+5;
int a[MX];
vector<PII>qu[MX];
int id[MX];
int bit[MX];
int ans[MX];
vector<int>v[MX];
void upd(int p,int c,int nn){
    while(p<=nn) bit[p]+=c,p=p+(p&-p);
}
int qur(int p){
    int re=0;
    while(p>0) re+=bit[p],p=p-(p&-p);
    return re;
}
int main(){
    int n,m;cin>>n>>m;
    for(int i=1;i<=n;i++)scanf("%d",&a[i]);
    for(int i=1;i<=m;i++){
        int sl,sr;scanf("%d%d",&sl,&sr);
        qu[sr].push_back(MP(sl,i));
    }
    for(int i=1;i<=n;i++){
        for(int j=i*2;j<=n;j=j+i){
            v[j].push_back(i);//j的因子
            v[i].push_back(j);//j的倍数
        }
    }
    for(int i=1;i<=n;i++){
        for(auto j:v[a[i]])if(id[j])upd(id[j],1,n);
        id[a[i]]=i;
        for(auto z:qu[i]){
            ans[z.SE]=( qur(i)-qur(z.FI-1) );
        }
    }
    for(int i=1;i<=m;i++)cout<<ans[i]<<'\n';
    return 0;
}

评论
添加红包

请填写红包祝福语或标题

红包个数最小为10个

红包金额最低5元

当前余额3.43前往充值 >
需支付:10.00
成就一亿技术人!
领取后你会自动成为博主和红包主的粉丝 规则
hope_wisdom
发出的红包
实付
使用余额支付
点击重新获取
扫码支付
钱包余额 0

抵扣说明:

1.余额是钱包充值的虚拟货币,按照1:1的比例进行支付金额的抵扣。
2.余额无法直接购买下载,可以购买VIP、付费专栏及课程。

余额充值