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]);
    }
}

 

發表評論
所有評論
還沒有人評論,想成為第一個評論的人麼? 請在上方評論欄輸入並且點擊發布.
相關文章