Codeforces Round #271 (Div. 2) F Ant colony

来源：互联网发布：怎么安装spss13软件编辑：程序博客网时间：2024/05/21 09:43

【题意】

给了n个数，然后给了m个查询，每个查询有个左区间l和右区间r。在这个区间里求出每一个数能别其他多少个数整除。如果，这个个数是r-l这个数就会被删除，要求个各区间里面不被删除的数的个数！

【解题方法】

如果区间[l,r]内所有数的gcd！=该区间内最小的数minv，则对于区间内的每个数都至少有一个数不是其倍数，那么每个数的v值都不可能等于r-l，所以所有的蚂蚁都要被吃掉（统统吃掉），答案是r-l+1。现在如果gcd=minv，那么对于所有的等于minv的数，所有的数都是他的倍数，对于所有的不等于minv的数，至少有一个不是他的倍数（比如minv），所以答案就是r-l+1-minv的个数。

【AC code】

//CF.271 F//Ant colony#include <cstdio>#include <cstring>#include <iostream>#include <algorithm>using namespace std;const int maxn=1e5+5;const int inf=1e9+10;struct node{    int l,r;    int minn,num,gc;//维护区间的最小值，最小值的个数，区间的最大公约数}Tree[maxn<<2];int gcd(int a,int b){    return b==0?a:gcd(b,a%b);}void PushUp(int rt){    Tree[rt].gc=gcd(Tree[rt*2].gc,Tree[rt*2+1].gc);}void Build(int l,int r,int rt){    Tree[rt].l=l,Tree[rt].r=r;    if(l==r)    {        scanf("%d",&Tree[rt].minn);        Tree[rt].num=1;        Tree[rt].gc=Tree[rt].minn;        return ;    }    int mid=(l+r)/2;    Build(l,mid,rt*2);    Build(mid+1,r,rt*2+1);    if(Tree[rt*2].minn<Tree[rt*2+1].minn)    {        Tree[rt].minn=Tree[rt*2].minn;        Tree[rt].num=Tree[rt*2].num;    }    else if(Tree[rt*2].minn>Tree[rt*2+1].minn)    {        Tree[rt].minn=Tree[rt*2+1].minn;        Tree[rt].num=Tree[rt*2+1].num;    }    else    {        Tree[rt].minn=Tree[rt*2].minn;        Tree[rt].num=Tree[rt*2].num+Tree[rt*2+1].num;    }    PushUp(rt);}int minv,minnum,mingcd;//记录[l,r]区间的答案void queryans(int L,int R,int rt){    if(L<=Tree[rt].l&&Tree[rt].r<=R)    {        if(minv>Tree[rt].minn)        {            minv=Tree[rt].minn;            minnum=Tree[rt].num;        }        else if(minv==Tree[rt].minn)        {//            minv=Tree[rt].minn;            minnum+=Tree[rt].num;        }        mingcd=gcd(mingcd,Tree[rt].gc);        return ;    }    int mid=(Tree[rt].l+Tree[rt].r)/2;    if(L<=mid) queryans(L,R,rt*2);    if(mid<R) queryans(L,R,rt*2+1);}int main(){    int n;    cin>>n;    Build(1,n,1);    //puts("success!");    int t,l,r;    cin>>t;    while(t--)    {        cin>>l>>r;        minv=inf;        minnum=0;        mingcd=0;        queryans(l,r,1);        //cout<<minv<<" "<<mingcd<<" "<<minnum<<endl;        if(mingcd==minv)        {            cout<<r-l+1-minnum<<endl;        }        else        {            cout<<r-l+1<<endl;        }    }}

0 0