hdu 4777 Rabbit Kingdom 离线树状数组

         杭州现场赛的H，现场的时候最后一个小时都在想这题，到底也没搞出来...当时搞出来怎么也Ag了...题意挺简单，给n个数m个查询，每个查询【L,R】，返回这个区间内有多少个数与其他任何数互质。先预处理一下，对每个数a【i】处理一个区间，使得在这个区间中a【i】和其他的数都互质，这个预处理挺简单，提前把1--20W每个数的因子存起来，然后正推一边反推一边就好了。之后从前到后遍历每个数，遍历到a[i]的时候，将l[i]-1的位置-1，r[i]+1位置的也应该要-1，但这个时候减掉的话会影响查询结果，所以把这个操作替换成“遍历到r[i]+1的时候，把i位置-1”,这个操作也用个vector数组 go[ ]记一下，等遍历到那个位置-1就好。因为有了这个操作，所以在遍历到a[i]的时候，第一件事应该是把go[i]当中存的要-1的操作执行掉，并且把l[go[i][j]]-1的位置+1 ，否则会导致同个位置重复被-1，导致结果偏小..遍历到一个右端R的时候，这段的长度再减掉这个区间里所有的-1就是这个查询的结果啦...

        

#include <iostream>
#include <cstdio>
#include <algorithm>
#include <cmath>
#include <string>
#include <cstring>
#include <queue>
#include <vector>
using namespace std;
typedef long long ll;
const int maxn=200000+20;
const int inf=(1<<30);
int n,m,p,q;
struct BIT
{
    int dt[maxn];
    void init()
    {
        memset(dt,0,sizeof dt);
    }
    int lowbit(int x)
    {
        return x&(-x);
    }

    void modify(int pos,int c)
    {
        if (pos<=0 || pos>n) return;

        for (int x=pos; x<=n; x+=lowbit(x))
        dt[x]+=c;
    }
    int query(int pos)
    {
        int res=0;
        for (int x=pos; x>0; x-=lowbit(x))
        res+=dt[x];
        return res;
    }
}bit;
struct node
{
    int x,y;
    int id;
    node()
    {

    }
    node(int a,int b)
    {
        x=a;
        y=b;
    }
    bool operator<(const node &p)const
    {
        if (y!=p.y) return y<p.y;
        return x<p.x;
    }
}que[maxn];
vector<int>go[maxn];
vector<int>fj[maxn];
int ans[maxn];
int a[maxn],l[maxn],r[maxn];
int last[maxn];

int main()
{
//    freopen("in.txt","r",stdin);
//    freopen("out.txt","w",stdout);
       for(int i=2;i<=200000;++i)
       for(int j=i;j<=200000;j+=i)
        fj[j].push_back(i);
    while(scanf("%d%d",&n,&m) && n && m)
    {
        for (int i=1; i<=n; i++)
        go[i].clear();
        memset(ans,0,sizeof ans);
        memset(last,0,sizeof last);
        memset(l,0,sizeof l);
        memset(r,0,sizeof r);
        for (int i=1; i<=n; i++)
        scanf("%d",&a[i]);

        for (int i=1; i<=n; i++)
        {

            int maxx=0;
            for (int j=0; j<fj[a[i]].size(); j++)
            {
                 maxx=max(last[fj[a[i]][j]],maxx);
                 last[fj[a[i]][j]]=i;
            }
            l[i]=maxx+1;
        }
        memset(last,0x3f3f,sizeof last);
        for (int i=n; i>=1; i--)
        {
            int minn=inf;
            for (int j=0; j<fj[a[i]].size(); j++)
            {
                 minn=min(last[fj[a[i]][j]],minn);
                 last[fj[a[i]][j]]=i;
            }
            if (minn>n) r[i]=n;
            else r[i]=minn-1;
        }
//        for (int i=1; i<=n; i++)
//        cout<<l[i]<<" "<r[i]<<endl;
        for (int i=0; i<m; i++)
        {
            scanf("%d%d",&que[i].x,&que[i].y);
            que[i].id=i;
        }
        sort(que,que+m);
        int num=0;
        bit.init();
        for (int i=1; i<=n; i++)
        {
            for (int j=0; j<go[i].size(); j++)
            {
                bit.modify(go[i][j],-1);
                bit.modify(l[go[i][j]]-1,1);
            }
            bit.modify(l[i]-1,-1);

            while(que[num].y==i && num<m)
            {
                int sub=bit.query(que[num].y)-bit.query(que[num].x-1);
                ans[que[num].id]=que[num].y-que[num].x+sub+1;
                num++;
            }
            if (num>=m) break;

            int tgt=r[i]+1;
            if (tgt<=n)
            {
                go[tgt].push_back(i);
            }
        }
        for (int i=0; i<m; i++)
        printf("%d\n",ans[i]);
//        tt--;
//        if (!tt)  break;
    }
    return 0;
}