积性函数与线性筛

这篇文章只是为了复习线性筛用的,所以非常水
主要是给莫反的文章准备个前置 给欧拉函数的文章填个坑

前置芝士:莫比乌斯函数、欧拉函数相关性质


积性函数

如果数论函数 f f f 保证,在 ( x , y ) = 1 (x,y)=1 (x,y)=1 时, f ( x y ) = f ( x ) × f ( y ) f(xy)=f(x)\times f(y) f(xy)=f(x)×f(y) ,那么 f f f 是积性函数。
特别地,如果对于任意 x , y x,y x,y 都有这个性质,则 f f f 为完全积性函数。

所有的积性函数都可以线性筛。下面会有几个线性筛的板子,先放上对应的函数:

  1. p ( x ) = [ i s p r i m e ( x ) ] p(x)=[isprime(x)] p(x)=[isprime(x)] x x x 是质数时值为1,否则为0,我也不知道叫啥名字)
  2. μ ( x ) \mu(x) μ(x) 莫比乌斯函数,不多介绍了
  3. φ ( x ) \varphi(x) φ(x) 欧拉函数,也不多说了

p ( x ) p(x) p(x) 其实没啥用,它对应的线性筛其实就是最原本的线性筛筛素数。就是用来保持队形的


线性筛

一种在线性复杂度内筛出 1 ∼ n 1\sim n 1n 中所有质数的算法。
在筛素数的基础上也可以筛各种积性函数。学习数论应该经常会用到这个技能

没啥技术含量 直接上板子就能看懂了 所以直接上板子吧

#include<cstdio>
const int n=1000000; 
int cnt,t,p[100000],v[n+10],mu[n+10],phi[n+10];
void wpri() { //prime
	for (int i=2; i<=n; ++i) {
		v[i]||(p[++cnt]=i);
		for (int j=1; j<=cnt&&1ll*p[j]*i<=n; ++j) {
			v[t=p[j]*i]=1;
			if (i%p[j]==0) break;
		}
	}
}
void wmu() { //mu
	for (int i=2; i<=n; ++i) {
		v[i]||(p[++cnt]=i,mu[i]=-1);
		for (int j=1; j<=cnt&&1ll*p[j]*i<=n; ++j) {
			v[t=p[j]*i]=1;
			if (i%p[j]==0) break;
			mu[t]=-mu[i];
		}
	}
}
void wphi() { //phi
	for (int i=2; i<=n; ++i) {
		v[i]||(p[++cnt]=i);
		for (int j=1; j<=cnt&&1ll*p[j]*i<=n; ++j) {
			v[t=p[j]*i]=1;
			if (i%p[j]) phi[t]=phi[i]*(p[j]-1);
			else { phi[t]=phi[i]*p[j]; break; }
		}
	}
}

int main() {
	//v[1]=1; wpri();
	//v[1]=mu[1]=1; wmu();
	//v[1]=phi[1]=1; wphi();
	return 0;
}

线性筛能干的事情中,筛这三类东西应该就是最常用的用途了
可以背板子 反正不难
另外下标为1的位置不要忘记赋个初值

好的这样一篇文章就写完了


评论
添加红包

请填写红包祝福语或标题

红包个数最小为10个

红包金额最低5元

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

抵扣说明:

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

余额充值