[ZOJ3886]Nico Number（线段树）

题目：

我是超链接

题意：

niconico数：所有比x小的数且与x互质的数，从小到大排列是一个等差数列，则x为niconico数。
长度为10^5的数组。三种操作：

第一：询问区间[l,r] niconico数的个数

第二：将区间 [ l , r ] 的数对v取模

第三：将某一个数更新为x

题解：

首先我们打个表看看这个数，然后我们可以发现，一个数是niconico数有以下3种情况

• 是质数
• 是2^x
• 是6

真是坑人，0&1也属于niconico。。

代码：

#include <cstdio>
#include <cstring>
#include <iostream>
using namespace std;
const int N=100005;
const int C=1e7;
bool vis[C+5],ss[C+5];int pri[C+5],sum[N*4],maxx[N*4],tot,a[N];
void init()
{
    for (int i=2;i<=C;i++)
    {
        if (!ss[i]) vis[i]=true,pri[++tot]=i;
        for (int j=1;j<=tot && pri[j]*i<=C;j++)
        {
            ss[pri[j]*i]=1;
            if (i%pri[j]==0) break;
        }
    }
    for (int i=2;i<=C;i*=2) vis[i]=1;
    vis[6]=vis[1]=vis[0]=1;
}
void updata(int now)
{
    sum[now]=sum[now<<1]+sum[now<<1|1];
    maxx[now]=max(maxx[now<<1],maxx[now<<1|1]);
}
void build(int now,int l,int r)
{
    if (l==r)
    {
        maxx[now]=a[l];
        sum[now]=(vis[a[l]]?1:0);
        return;
    }
    int mid=(l+r)>>1;
    build(now<<1,l,mid);
    build(now<<1|1,mid+1,r);
    updata(now);
}
int qurry(int now,int l,int r,int lrange,int rrange)
{
    if (lrange<=l && rrange>=r) return sum[now];
    int mid=(l+r)>>1,ans=0;
    if (lrange<=mid) ans+=qurry(now<<1,l,mid,lrange,rrange);
    if (rrange>mid) ans+=qurry(now<<1|1,mid+1,r,lrange,rrange);
    return ans;
}
void change(int now,int l,int r,int x,int v)
{
    if (l==r)
    {
        maxx[now]=v;
        sum[now]=(vis[maxx[now]]?1:0);
        return;
    }
    int mid=(l+r)>>1;
    if (x<=mid) change(now<<1,l,mid,x,v);
    else change(now<<1|1,mid+1,r,x,v);
    updata(now);
}
void qm(int now,int l,int r,int lrange,int rrange,int x)
{
    if (maxx[now]<x) return;    
    if (l==r)
    {
        maxx[now]%=x;
        sum[now]=(vis[maxx[now]]?1:0);
        return;
    }
    int mid=(l+r)>>1;
    if (lrange<=mid) qm(now<<1,l,mid,lrange,rrange,x);
    if (rrange>mid) qm(now<<1|1,mid+1,r,lrange,rrange,x);
    updata(now);
}
int main()
{
    init();int n,m,id;
    while (scanf("%d",&n)!=EOF)
    {
        for (int i=1;i<=n;i++) scanf("%d",&a[i]);
        build(1,1,n);
        scanf("%d",&m);
        for (int i=1;i<=m;i++)
        {
            int l,r,x;scanf("%d",&id);
            if (id==1)
            {
                scanf("%d%d",&l,&r);printf("%d\n",qurry(1,1,n,l,r));
            }else if (id==2)
            {
                scanf("%d%d%d",&l,&r,&x);
                qm(1,1,n,l,r,x);
            }else
            {
                scanf("%d%d",&l,&x);
                change(1,1,n,l,x);
            }
        }
    }
}