hdu4777 树状数组

题意:给了n个数，然后又m次查询，询问[L,R] 内有多少个数与其他的数不互质。

解:

我们首先可以通过处理得出每个数的有效区间，LR 就是 左边L位置上的数 和他不互质， 右边R位置上的数和不互质，

我们对于询问排序，R小的排前面，枚举每个R，在loc位置就将第loc个点在loc的位置加上一个1在loc的L(左不互质点)减一个1,再将枚举到该位的时候对于有在这个位置上R的点 在loc位置减1

在loc的L位置加1

#include <iostream>
#include <algorithm>
#include <cstdio>
#include <string.h>
#include <vector>
#include <cmath>
using namespace std;
const int maxn=;
typedef long long LL;
bool vis[maxn];
int yinzi[maxn][],numofYZ[maxn];
void sieve()
{

     memset(vis,false,sizeof(vis));
     memset(numofYZ,false,sizeof(numofYZ));
    for(LL i=; i<=; i++)
    {
        if(vis[i])continue;
        yinzi[i][numofYZ[i]++]=i;
         for(LL j=i+i; j<=; j+=i)
         {
            vis[j]=true;
            yinzi[j][numofYZ[j]++]=i;
         }
    }
}
struct point{
   int L,R,id;
   bool operator <(const point &rhs)const {
     return R<rhs.R||( R == rhs.R && L<rhs.L);
   }
}wLR[maxn],Q[maxn];
int Loc[maxn];
int w[maxn];
vector<int>G[maxn];
void init(int n)
{
    for(int i=; i<=n+; i++)
        G[i].clear();
}
int n,m;
int C[maxn];
int lowbit(int x)
{
    return x&(-x);
}
void add(int x,int v)
{
    if(x<=)return ;
    while(x<=n)
    {
       C[x]+=v;
       x+=lowbit(x);
    }
}
int sum(int x)
{
    int ans=;
    while(x>)
    {
        ans+=C[x];
        x-=lowbit(x);
    }
    return ans;
}
int ans[maxn];
int main()
{
    sieve();
    while(scanf("%d%d",&n,&m)==&&n)
    {
        int maW=;
        for(int i=; i<=n; i++)
            {
                scanf("%d",&w[i]);
               maW=max(maW,w[i]);
            }
        memset(Loc,,sizeof(Loc));
        for(int i=; i<=n; i++)
        {
            int ww=w[i];
            int L=;
            for(int j=; j<numofYZ[ ww ]; j++)
                L=max(L,Loc[ yinzi[ ww ][ j ] ]);
            wLR[i].L=L;
            for(int j=; j<numofYZ[ ww ]; j++)
                Loc[ yinzi[ ww ][ j ] ] = i;
        }
        for(int i=; i<=maW; i++)Loc[i]=n+;
        for(int i=n; i>; i--)
        {
            int ww=w[i];
            int R=n+;
            for(int j=; j<numofYZ[ ww ]; j++)
                R=min(R,Loc[ yinzi[ ww ][ j ] ]);
            wLR[i].R=R;
            G[R].push_back(i);
            for(int j=; j < numofYZ[ ww ]; j++)
                Loc[ yinzi[ww][ j ] ]=i;
        }
        for(int i=; i<m; i++)
        {
                Q[i].id=i;
                scanf("%d%d",&Q[i].L,&Q[i].R);
        }
        sort(Q,Q+m);
        int now=;
        memset(C,,sizeof(C));
        for(int i=; i<m; i++)
        {
            while(now<=Q[i].R)
            {
              for(int j=; j<G[now].size(); j++)
              {
                  int to=G[now][j];
                  add(to,-);
                  add(wLR[to].L,);
              }
              add(now,);
              add(wLR[now].L,-); now++;
            }
            ans[Q[i].id]=sum(Q[i].R)-sum(Q[i].L-);
        }
        for(int i=; i<m; i++)
            printf("%d\n",ans[i]);
        init(n);
    }
    return ;
}