hdu4777——Rabbit Kingdom

题目链接
http://acm.hdu.edu.cn/showproblem.php?pid=4777
题目大意
长度为n的序列，m次询问，每次询问 [l,r] 中与其它数不互质的数有几个。
解题思路
离线树状数组。将询问按右端点 r 从小到达排序。预处理出每个数左端与其最近的互质数位置，以及右端与其最近的互质数位置。然后维护一个树状数组，表示不合法的数的数量。最后答案就是当前询问长度 r-l+1 减去 当前区间不合法的数的数量。
树状数组维护思路，每次到一个点 i ，将其 l[i] 位置树状数组 +1，因为只有l[i]在询问区间里，i这个点才算是不合法的数。然后将每个 j (r[j]=i) 位置树状数组 +1，并且 l[j] 树状数组 -1，因为询问区间是按右端点 r 从小到大排序的，所以其 l[j] 位置肯定比 r[j] 位置先不在询问区间内。
AC代码

#include <bits/stdc++.h>
using namespace std;
const int N=2e5;
int n,m;
int a[N+10];
bool vis[N+10];
int tmp[N+10];
int l[N+10];
int r[N+10];
int tree[N+10];
int cnt;
int su[N];
vector<int>mp[N+10];
vector<int>vc[N+10];
struct node
{
    int l,r,id;
}q[N+10];
int ans[N+10];
bool cmp(node a,node b)
{
    return a.r<b.r;
}
void get_prime()
{
    for(int i=2;i<=N;++i)
    {
        if(!vis[i])
            su[++cnt]=i;
        for(int j=1;j<=cnt&&su[j]*i<=N;++j)
            vis[su[j]*i]=1;
    }
    for(int i=1;i<=cnt;++i)
        for(int j=1;j*su[i]<=N;++j)
            mp[j*su[i]].push_back(i);
}
void init()
{
    for(int i=0;i<=n;++i)
    {
        l[i]=0;
        r[i]=n+1;
        vc[i].clear();
    }
    memset(tree,0,sizeof(tree));
    vc[n+1].clear();
}
int lowbit(int i)
{
    return i&(-i);
}
void add(int k,int x)
{
    if(!k)
        return ;
    while(k<=n)
    {
        tree[k]+=x;
        k+=lowbit(k);
    }
}
int ask(int k)
{
    int res=0;
    while(k)
    {
        res+=tree[k];
        k-=lowbit(k);
    }
    return res;
}
int main()
{
    get_prime();
    while(~scanf("%d%d",&n,&m) && (n||m))
    {
        for(int i=1;i<=n;++i)
            scanf("%d",&a[i]);
        for(int i=1;i<=m;++i)
        {
            scanf("%d%d",&q[i].l,&q[i].r);
            q[i].id=i;
        }
        sort(q+1,q+1+m,cmp);
        init();
//        puts("1");
        memset(tmp,0,sizeof(tmp));
        for(int i=1;i<=n;++i)
        {
            for(int j=0;j<mp[a[i]].size();++j)
            {
                l[i]=max(l[i],tmp[mp[a[i]][j]]);
                tmp[mp[a[i]][j]]=i;
            }
        }
        for(int i=1;i<=cnt;++i)
            tmp[i]=n+1;
        for(int i=n;i>=1;i--)
        {
            for(int j=0;j<mp[a[i]].size();++j)
            {
                r[i]=min(r[i],tmp[mp[a[i]][j]]);
                tmp[mp[a[i]][j]]=i;
            }
        }
        for(int i=1;i<=n;++i)
            vc[r[i]].push_back(i);
        int now=1;
//        puts("1");
        for(int i=1;i<=m;++i)
        {
            while(now<=q[i].r)
            {
                add(l[now],1);
//                puts("1");
                for(int j=0;j<vc[now].size();++j)
                {
                    add(vc[now][j],1);
                    add(l[vc[now][j]],-1);
                }
                now++;
            }
//            puts("1");
            ans[q[i].id]=q[i].r-q[i].l+1-(ask(q[i].r)-ask(q[i].l-1));
        }
        for(int i=1;i<=m;++i)
            printf("%d\n",ans[i]);
    }
    return 0;
}