推荐视频:518 筛法求莫比乌斯函数
前提知识:莫比乌斯函数
点击查看代码
#include<bits/stdc++.h>
using namespace std;
#define LL long long
const int N = 1e8 + 10;
int p[N], cnt;
int mu[N];//d[i]记录i的约数的和
bool vis[N];
void get_mu(int n) {//筛法求约数的个数
mu[1] = 1;
for (int i = 2; i <= n; i++) {
if (!vis[i]) {
p[++cnt] = i;
mu[i]=-1;
}
for (int j = 1; 1LL * i * p[j]<=n; j++) {
int m = i * p[j];
vis[m]=1;
if (i % p[j] == 0) {//含有相同质因子
mu[m] = 0;
break;
}
else {//若i不能被p[j]整除
mu[m] = -mu[i];
}
}
}
}
int main() {
int n;
cin >> n;
get_mu(n);
return 0;
}