素数筛法

ll pri[maxn],pri_num;
ll mu[maxn];//莫比乌斯函数值
bool vis[maxn];

void mobius()  {//筛法求莫比乌斯函数
    pri_num = 0;//素数个数
    memset(vis, false, sizeof(vis));
    vis[1] = true;
    mu[1] = 1;
    for(int i = 2; i <maxn; i++){
        if(!vis[i]){
            pri[pri_num++] = i;
            mu[i] = -1;
        }
        for(int j=0; j<pri_num && i*pri[j]<maxn ; j++){
            vis[i*pri[j]]=true;//标记非素数
            //eg:i=3,i%2,mu[3*2]=-mu[3]=1;----;i=6,i%5,mu[6*5]=-mu[6]=-1;
            if(i%pri[j])mu[i*pri[j]] = -mu[i];
            else {
                mu[i*pri[j]] = 0;
                break;
            }

        }
    }
}

 

相关文章:

  • 2022-03-11
  • 2021-10-24
  • 2022-02-26
猜你喜欢
  • 2021-06-13
  • 2022-12-23
  • 2021-08-18
  • 2022-12-23
  • 2021-06-21
  • 2021-12-08
相关资源
相似解决方案