codeforces 474F Ant colony(线段树)

来源：互联网发布：没货源怎么开淘宝店编辑：程序博客网时间：2024/05/18 13:48

对于每个查询的区间，能活下来的蚂蚁的个数就是这段区间中等于这段区间gcd的s的个数，而gcd一定小于等于这段区间的最小值，所以可用线段树来维护每个区间的gcd，最小值和最小值出现的次数，如果gcd等于最小值，那么这段区间中能活下来的蚂蚁的个数就是最小值出现的次数。

#include <iostream>#include<cstdio>#include<cstring>using namespace std;int s[100005];int gcd(int a,int b){    if(!b) return a;    return gcd(b,a%b);}struct A{    int l,r,d,mn,times;}a[4*100005];void build(int k,int l,int r){    a[k].l=l;    a[k].r=r;    if(l==r) {        a[k].mn=s[l];        a[k].times=1;        a[k].d=s[l];        return ;    }    int mid=(l+r)>>1;    build(k<<1,l,mid);    build(k<<1|1,mid+1,r);    a[k].d=gcd(a[k<<1].d,a[k<<1|1].d);    if(a[k<<1].mn<a[k<<1|1].mn){        a[k].mn=a[k<<1].mn;        a[k].times=a[k<<1].times;    }    else if(a[k<<1].mn>a[k<<1|1].mn){        a[k].mn=a[k<<1|1].mn;        a[k].times=a[k<<1|1].times;    }    else{        a[k].mn=a[k<<1].mn;        a[k].times=a[k<<1].times+a[k<<1|1].times;    }}int ans;struct B{    int d,mn,times;};B query(int k,int ql,int qr){    B tmp;    if(ql<=a[k].l && a[k].r<=qr){        tmp.d=a[k].d;        tmp.mn=a[k].mn;        tmp.times=a[k].times;        return tmp;    }    int mid=(a[k].l+a[k].r)>>1;    B left,right;    if(qr<=mid) return query(k<<1,ql,qr);    else if(mid<ql) return query(k<<1|1,ql,qr);    else{        left=query(k<<1,ql,qr);        right=query(k<<1|1,ql,qr);        tmp.d=gcd(left.d,right.d);        if(left.mn<right.mn){            tmp.mn=left.mn;            tmp.times=left.times;        }        else if(left.mn>right.mn){            tmp.mn=right.mn;            tmp.times=right.times;        }        else{            tmp.mn=left.mn;            tmp.times=left.times+right.times;        }        return tmp;    }}int check(B x){    if(x.d==x.mn) return x.times;    else return 0;}int main(){    int n,t,i;    scanf("%d",&n);    for(i=1;i<=n;i++){        scanf("%d",s+i);    }    build(1,1,n);    scanf("%d",&t);    for(i=1;i<=t;i++){        int l,r;        scanf("%d%d",&l,&r);        printf("%d\n",(r-l+1)-check( query(1,l,r) ));    }    return 0;}

0 0