莫比乌斯函数

推荐视频: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;
}
posted @ 2023-09-29 20:32  不o凡  阅读(18)  评论(0)    收藏  举报