2013年 ACMICPC 杭州赛区H题

思路：树状数组统计。待验证，不知道是否对。 

#include<cstdio>
#include<cstring>
#include<cmath>
#include<iostream>
#include<algorithm>
#define lowbit(x) (x&(-x))
#define Maxn 200010
using namespace std;
int C[Maxn],vi[Maxn],pre[Maxn],n,num[Maxn],ans[Maxn],ov[Maxn];
struct Qu{
    int l,r,i;
}q[Maxn];
int cmp(Qu a,Qu b)
{
    return a.r<b.r;
}
void update(int pos,int val)
{
    while(pos){
        C[pos]+=val;
        pos-=lowbit(pos);
    }
}
int sum(int pos)
{
    int s=0;
    while(pos<=n){
        s+=C[pos];
        pos+=lowbit(pos);
    }
    return s;
}
int main()
{
    int m,i,j;
    while(scanf("%d%d",&n,&m),n||m){
        memset(C,0,sizeof(C));
        memset(vi,0,sizeof(vi));
        memset(pre,0,sizeof(pre));
        memset(ov,0,sizeof(ov));
        for(i=1;i<=n;i++)
            scanf("%d",num+i);
        for(i=1;i<=m;i++){
            scanf("%d%d",&q[i].l,&q[i].r);
            q[i].i=i;
        }
        sort(q+1,q+1+m,cmp);
        int r=1,temp,pos;
        for(i=1;i<=n;i++){
            update(i,1);
            pos=0;
            for(j=2;j*j<=num[i];j++){
                if(num[i]%j) continue;
                if(!vi[pre[j]]&&pre[j]) {
                    update(pre[j],-1);
                    if(ov[pre[j]])
                        update(ov[pre[j]],1);
                    vi[pre[j]]=1;
                }
                pos=max(pos,pre[j]);
                pre[j]=i;
                if(j*j==num[i]) continue;
                temp=num[i]/j;
                if(!vi[pre[temp]]&&pre[temp]){
                    update(pre[temp],-1);
                    if(ov[pre[temp]])
                        update(ov[pre[temp]],1);
                    vi[pre[temp]]=1;
                }
                pos=max(pos,pre[temp]);
                pre[temp]=i;
            }
            if(!vi[pre[num[i]]]&&pre[num[i]]){
                update(pre[num[i]],-1);
                if(ov[pre[num[i]]])
                        update(ov[pre[num[i]]],1);
                vi[pre[num[i]]]=1;
            }
            pos=max(pos,pre[num[i]]);
            pre[num[i]]=i;
            if(pos){
                update(pos,-1);
                ov[i]=pos;
            }
            while(r<=m&&q[r].r==i){
                ans[q[r].i]=sum(q[r].l);
                r++;
            }
        }
        for(i=1;i<=m;i++)
            printf("%d
",ans[i]);
    }
    return 0;
}