线性筛法求素数(埃氏筛法、欧拉筛法)

筛法都是初始化把所有数都先设为素数,然后筛除合数。

1)埃氏筛法

理解起来比较简单,就是从小到大的枚举每一个数,标记它的所有倍数都是合数(非素数)放到u[ i ]为false,u中的下标对应的就是这个数的值。

bool u[MAXN];
int num,su[MAXN];
// 1. 埃氏筛法
void Prime () { 
	for (int i=2; i<=MAXN; i++) 
		u[i]=true;
	for (int i=2; i<=MAXN; i++) {
		if(u[i])
			for(int j=2;j*i<=MAXN;j++)
			u[j*i]=false;
		
		for (int i=2; i<=MAXN; i++) 
			if(u[i]) su[num++]=i;
		
	}
}

2)欧拉筛法

从小到大的筛选素数,事实上还是会有很多重复筛选的情况,所以欧拉筛选就出现了,每个合数只会被筛一遍。
其思想就是:每个合数只被它的最小质因子筛选一次。

void findSu() {
	//memset(u, true, sizeof(u));
	for(int i=0;i<=MAXN;i++) u[i]=true;
	for (int i=2; i<=MAXN; i++) {
		if(u[i]) su[num++]=i;   //素数就是true
		
		for (int j=0; j<num && i*su[j]<=MAXN; j++) {
			u[i*su[j]]=false;
			if(i%su[j]==0) break; 
		}
	}
}

解释一下if(i%su[j]==0) break; 的意思:
因为任何合数都能表示成多个素数的积.所以,每个合数都有一个最小质因子.我们通过这个最小质因子就可以判断什么时候不用筛选下去了.

i % su[j]==0, i是su[j]的倍数时,如果继续循环,对于i * su[j+1]这个数一定在后面还会被筛一次,为什么呢?

因为这个时候 i=su[j]* k(k是某个正整数),所以 i * su[j+1] = su[j]* k* su[j+1].而su[j]必定小于su[j+1],所以i* su[j+1]必定被后面当i=su[j+1]* k的时候再一次筛掉,所以我们break;

例如 :i = 8 ,j = 1,su[j] = 2,如果不跳出循环,su[j+1] = 3,8 * 3 = 2 * 4 * 3 = 2 * 12,在i = 12时会重复筛选一次。欧拉筛法的原理便是通过最小素因子来消除。

总结

埃氏筛法是从小到大枚举每个数的倍数,而欧拉筛法从小到大是枚举每个素数的倍数。

  • 0
    点赞
  • 2
    收藏
    觉得还不错? 一键收藏
  • 0
    评论

“相关推荐”对你有帮助么?

  • 非常没帮助
  • 没帮助
  • 一般
  • 有帮助
  • 非常有帮助
提交
评论
添加红包

请填写红包祝福语或标题

红包个数最小为10个

红包金额最低5元

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

抵扣说明:

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

余额充值