常见素数筛法

列出几种常用的素数筛选法,附上计时器。。。

#include<cstdio>
#include<cstdlib>
#include<cmath>
#include<map>
#include<queue>
#include<stack>
#include<vector>
#include<algorithm>
#include<cstring>
#include<string>
#include<iostream>
#include<windows.h>
#define ms(x,y) memset(x,y,sizeof(x))
const int MAXN=90000000;
const int INF=1<<30;
const int sq=(int)sqrt(double(MAXN));  
using namespace std;
int u;
int len;
bool vis[MAXN+5];
int prime[1000000];

/*
 * 奇数筛法,(偶数肯定不是素数,故只需要判断奇数即可)即是vis数组里面虚拟的值都是奇数
 */
void print_prime()  
{  
	if(MAXN&1==0) len=MAXN/2-1;
	else len=(MAXN+1)/2-1;
	bool *p;
    for(int i=0; i*i<len; i++){
        if(vis[i]){
			int b=2*i+3;//第i个奇数的值
            for(p=vis+i+b; p<vis+len; p+=b) *p=false;
        }  
    }  
}  

/*
 * 普通筛法
 */
void print_prime1()
{
	for(int i=2; i*i<=MAXN; i++){
		if(vis[i]==true){
			for(int j=i*i; j<=MAXN; j+=i) vis[j]=false;
		}
	}
}

/*
 * 减少重复筛选
 */
void print_prime2()
{
	for(int i=2; i<=MAXN; i+=2)
		vis[i] = false;
	vis[2] = true;
	for(int i=3; i*i<=MAXN; i++){
		if(vis[i]==true){
			for(int j=i*i; j<=MAXN; j+=2*i) vis[j]=false;
		}
	}
}

/*
 * 小数据常用方法
 */
void print_prime3()
{
	for(int i=2; i<=MAXN; i++){
		int ok=1;
		for(int j=2; j*j<=i; j++){
			if(i%j==0){
				ok=0;
				break;
			}
		}
	}
}

int main()
{
	ms(vis,true);

	FILETIME beg,end;//<windows.h>里的计时器
	GetSystemTimeAsFileTime(&beg);
	print_prime2(); 
	GetSystemTimeAsFileTime(&end);
	long time = 100*(end.dwLowDateTime-beg.dwLowDateTime);
	cout<<time<<endl;
#if 0
	for(int i=2; i<=MAXN; i++) if(vis[i]==true) u++; 
	cout<<u<<endl;
#endif
	u=1;
	for(int i=0; i<len; i++) if(vis[i]==true) u++; 
	cout<<u<<endl;
	return 0;
}


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

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

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

请填写红包祝福语或标题

红包个数最小为10个

红包金额最低5元

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

抵扣说明:

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

余额充值