hdu4777 树状数组

 

题目链接：http://acm.hdu.edu.cn/showproblem.php?pid=4777

 

题目大意：

给n个数，求LR区间内有多少个数和当前区间所有其他数互质。

 

解题思路：

训练赛完全想不到，首先我们需要预处理出 对于一个数 i 它最左和最右到哪不互质，记为 pl pr，

那么当我们查询区间 L R时，实际上就是查询区间中有多少 pl<L pr>R的数。

首先基本没有涉及修改操作的题目，几乎都是需要的离线的，这题也不例外，做法也比较多，我们先按照R排序，

当我们遇到i的时候，在pl处+1，在pr处-1，遇到pr的时候，在i处+1,pl处-1。

这其中个人感觉很难理解，最好自己推一下样例，应该会容易理解很多。

 

Ac代码：

#include<bits/stdc++.h>
using namespace std;
const int maxn=2e5+5;
bool vis[maxn];
vector<int> pri[maxn],vec[maxn];    //储存每个数的素因子
int n,m,a[maxn],pos[maxn],ans[maxn];
struct node
{
    int l,r;
    int idx;
    bool operator<(const node &p) const
    {
        return r<p.r;
    }
}pa[maxn],q[maxn];
void init(int si)   //得到每个数的素因子
{
    for(int i=0;i<maxn;i++) pri[i].clear();
    for(int i=2;i<maxn;i++)
        if(pri[i].size()==0)
            for(int j=i;j<maxn;j+=i)
                pri[j].push_back(i);
}
void vs()   //找出每个数左边界和右边界
{
    for(int i=1;i<=n;i++) pa[i].l=0,pa[i].r=n+1;
    for(int i=0;i<maxn;i++) pos[i]=0;
    for(int i=1;i<=n;i++)
    {
        int now=a[i];
        for(int j=0;j<pri[now].size();j++)
        {
            int t=pri[now][j];
            pa[i].l=max(pa[i].l,pos[t]);
            pa[pos[t]].r=min(pa[pos[t]].r,i);
            pos[t]=i;
        }
    }
}

int bit[maxn];
int lowbit(int x) { return x&(-x); }
void add(int x,int val)
{
    if(x==0) return ;
    while(x<=n)
    {
        bit[x]+=val;
        x+=lowbit(x);
    }
}
int sum(int x)
{
    int s=0;
    while(x)
    {
        s+=bit[x];
        x-=lowbit(x);
    }
    return s;
}

int main()
{
    init(maxn);
    while(scanf("%d%d",&n,&m)!=EOF)
    {
        if(n==0&&m==0) break;
        memset(bit,0,sizeof bit);
        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].idx=i;
        vs();
        sort(q+1,q+1+m);
        for(int i=0;i<n+1;i++) vec[i].clear();
        for(int i=1;i<=n;i++) vec[pa[i].r].push_back(i);
        int id=1;
        for(int i=1;i<=m;i++)   //区间处理
        {
            while(id<=n&&id<=q[i].r)
            {
                add(pa[id].l,1);
                int sz=vec[id].size();
                for(int j=0;j<sz;j++)
                {
                    int v=vec[id][j];
                    add(pa[v].l,-1);
                    add(v,1);
                }
                id++;
            }
            ans[q[i].idx]=sum(q[i].r)-sum(q[i].l-1);
            ans[q[i].idx]=q[i].r-q[i].l+1-ans[q[i].idx];
        }
        for(int i=1;i<=m;i++) printf("%d\n",ans[i]);
    }
}