hdu-4777Rabbit Kingdom 树状数组

http://acm.split.hdu.edu.cn/showproblem.php?pid=4777

题意：给一段序列，多个查询，查询某段区间内和其他数都互质的数的个数

接触过类似的题目，但训练时竟没有想到树状数组，搜了发题解，看到树状数组四个字，立马就把题解关掉了，想了想，好像不难做

必然是离线处理，套路类似spoj  D-query

需要修改的地方就是，以因子为索引寻找前驱，

w【i】原数组

pre【i】记录当前扫描 因子包含 i 的最左边的数的下标

pre1【i】记录 i 左边的最右边的和w【i】不互质数的下标

首先对于这种贡献，直接对树状数组某一个点加1是不够的，还需要对左边最近的那个不互质的数减1

从左到右一个一个数处理，通过对w【i】质因子分解，得到因子x，pre【x】找到前面和w【i】不互质的数，这个时候他们是不能有贡献的，要还原，加的减回来，减的加回来，这也是记录pre1【i】的用处，最后找到一个 i 左边最近的和w【i】不互质的数（下标为j），树状数组里 i对应位置加1，j对应位置减1，处理右端点为i的查询

然后T了。。。。。

仔细分析了一波，改了改，还是T

看了一下别人代码，顿时醒悟，质因子分解取模的地方的确有可能T，需要类似筛法预处理所有数的素因子 1466msAC

关键地方全不是自己想的QAQ

#include<bits/stdc++.h>
using namespace std;
const int maxn=200005;
struct node
{
    int le,ri,id;
    bool operator < (const node u)const
    {
        if(ri==u.ri)
            return le<u.le;
        return ri<u.ri;
    }
}Q[maxn];
int n,m;
int w[maxn],f[maxn],pre[maxn],pre1[maxn],book[maxn];
vector<int>factor[maxn];
void init()
{
    for(int i=2;i<=200000;i++)
    {
        if(factor[i].empty())
        {
            for(int j=i;j<=200000;j+=i)
            {
                factor[j].push_back(i);
            }
        }
    }
}
int sum(int x)
{
    int ans=0;
    while(x)
    {
        ans+=f[x];
        x-=x&-x;
    }
    return ans;
}
void add(int x,int k)
{
    if(x==0)
        return ;
    while(x<=n)
    {
        f[x]+=k;
        x+=x&-x;
    }
    return ;
}
void update(int x,int pp)
{
    if(book[pre[pp]])
    {
        add(pre[pp],-1);
        book[pre[pp]]=0;
        add(pre1[pre[pp]],1);
    }
    pre[pp]=x;
}
int main()
{
    int i,j;
    int ans[maxn];
    init();
    while(cin>>n>>m,n+m)
    {
        for(i=1;i<=n;i++)
        {
            scanf("%d",&w[i]);
        }
        for(i=0;i<m;i++)
        {
            scanf("%d %d",&Q[i].le,&Q[i].ri);
            Q[i].id=i;
        }
        sort(Q,Q+m);
        memset(f,0,sizeof(f));
        memset(pre,0,sizeof(pre));
        memset(pre1,0,sizeof(pre1));
        memset(book,0,sizeof(book));
        for(j=0,i=1;i<=n&&j<m;i++)
        {
            int temp=w[i],least=0;
            for(int k=0;k<factor[temp].size();k++)
            {
                int now=factor[temp][k];
                least=max(least,pre[now]);
                update(i,now);
            }
            add(i,1);
            add(least,-1);
            book[i]=1;
            pre1[i]=least;

            while(Q[j].ri==i&&j<m)
            {
                ans[Q[j].id]=sum(i)-sum(Q[j].le-1);
                j++;
            }
        }
        for(i=0;i<m;i++)
            printf("%d\n",ans[i]);
    }
    return 0;
}