codeforces #271(div2) F. Ant colony

题意：一个数列，q个询问，问某个区间，如果某个数是其他数的因子，则s++,问r-l+1-s是多少

思路：一个数是其他数的因子，这个数肯定是最小数，且为这个区间的gcd，所以RMQ求区间GCD，区间最小值，在用线段树求这个区间最小值个数

#include<bits/stdc++.h>
using namespace std;
typedef long long ll;
const int N=1e5+10;

ll dp1[N][20],dp2[N][20],a[N];
int n;

void init(){
    memset(dp2,127,sizeof(dp2));
    for(int i=1;i<=n;i++)
        dp1[i][0]=dp2[i][0]=a[i];
    for(int j=1;(1<<j)<=n;j++){
        for(int i=1;i+(1<<j)-1<=n;i++)
        {
            dp1[i][j]=__gcd(dp1[i][j-1],dp1[i+(1<<(j-1))][j-1]);
            dp2[i][j]=min(dp2[i][j-1],dp2[i+(1<<(j-1))][j-1]);
        }
    }
}
int  check(int x,int y){
    int kk=0;
    while((1<<(kk+1))<=y-x+1) kk++;
    int Min=min(dp2[x][kk],dp2[y-(1<<kk)+1][kk]);
    return Min;
}
int  check1(int x,int y){
    int kk=0;
    while((1<<(kk+1))<=y-x+1) kk++;
    int Min=__gcd(dp1[x][kk],dp1[y-(1<<kk)+1][kk]);
    return Min;
}
struct node{
    int l,r,M,num;
}e[N*4];

void build(int x,int L,int R){
    e[x].l=L;e[x].r=R;

    if(L==R){
        e[x].M=a[L];e[x].num=1;return ;
    }
    int mid=(L+R)>>1;
    build(2*x,L,mid);
    build(2*x+1,mid+1,R);
    if(e[2*x].M==e[2*x+1].M){
        e[x].M=e[2*x].M;
        e[x].num=e[x*2].num+e[2*x+1].num;
    }
    else {
        if(e[2*x].M<e[2*x+1].M){
            e[x].M=e[2*x].M;e[x].num=e[2*x].num;
        }
        else
        {
             e[x].M=e[2*x+1].M;e[x].num=e[2*x+1].num;
        }
    }
}

int sum;
void  query(int x,int L,int R,int y){
    if(e[x].l>=L&&R>=e[x].r){
        if(e[x].M==y) {
            sum+=e[x].num;
        }
        return ;
    }
    int mid=(e[x].l+e[x].r)>>1;
    if(mid>=R)  query(2*x,L,R,y);
    else if(L>mid) query(2*x+1,L,R,y);
    else {
        query(2*x,L,mid,y);
        query(2*x+1,mid+1,R,y);
    }
}
int main(){
    scanf("%d",&n);
    for(int i=1;i<=n;i++){
        scanf("%d",&a[i]);
    }
    build(1,1,n);
    init();
    int q;

    cin>>q;
    int x,y;
    while(q--){
        scanf("%d%d",&x,&y);
        ll z=check(x,y);
        ll zz=check1(x,y);
       // cout<<z<<" "<<zz<<endl;
        if(z==zz){
            sum=0;
            query(1,x,y,zz);
            printf("%d\n",y-x+1-sum);
        }
        else {
            printf("%d\n",y-x+1);
        }
    }
}