素数的筛法,n*long(n)和O(n)的写法


需求:

产生区间[0, N]之间的素数。


最简单的方法是,一个一个套,但是每试探一个数n, 需要枚举sqrt(n)次,所以这种算法的时间复杂度就是O(N*sqrt(N)), 也就是N*N。

那么效率更加高一点的方法是,把每个素数的2,3,4,... x, (x <= sqrt(N)) 倍数筛选出来,留下的就是素数。 该算法效率是 N*log(N)。

void sieveOne(){
	int m = sqrt(N+0.5);
	for(int i=2; i<= m; i++){
		if(!arr[i]){
			for(int j=i*i; j<= N; j+=i){
				arr[j] = 1;
			}
		}
	}
}

其实还有更加高效的算法。 根据因式分解定理,每一个合数都可以分解为 N= p^k*q的形式,又因为每个合数的分解方式是唯一的,所以该算法的效率是O(N)。


void sieveTwo(){
	int p,q;
	int end = sqrt(N+0.5);
	for(p = 2; p!=end; ++p ){
		if(!arr[p]){
			for(q = p; p*q <= N; ++q){
				for(int k= p*q; k<= N; k*=p) arr[k] = 1;
			}
		}
	}
}




评论
添加红包

请填写红包祝福语或标题

红包个数最小为10个

红包金额最低5元

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

抵扣说明:

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

余额充值