莫比乌斯函数线性筛

int mo[maxn], vis[maxn], prime[maxn];

void init () {
	mo[1] = 1;
	mem(vis, 0);
	int len = 0;
	for (int i = 2; i < maxn; ++i) {
		if (!vis[i]) {
			prime[len++] = i;
			mo[i] = -1;
		}
		for (int j = 0; j < len && (LL)prime[j] * i < maxn; ++j) {
			vis[ i * prime[j] ] = 1;
			if (i % prime[j] == 0) {
				mo[ i * prime[j] ] = 0;
				break;
			}
			mo[ i * prime[j] ] = -mo[i];
		}
	}
	// for (int i = 2; i < maxn; ++i) mo[i] += mo[i-1];
}
评论 2
添加红包

请填写红包祝福语或标题

红包个数最小为10个

红包金额最低5元

当前余额3.43前往充值 >
需支付:10.00
成就一亿技术人!
领取后你会自动成为博主和红包主的粉丝 规则
hope_wisdom
发出的红包
实付
使用余额支付
点击重新获取
扫码支付
钱包余额 0

抵扣说明:

1.余额是钱包充值的虚拟货币,按照1:1的比例进行支付金额的抵扣。
2.余额无法直接购买下载,可以购买VIP、付费专栏及课程。

余额充值