常见素数筛法

标签：io   os   sp   for   on   数据   ef   amp   as   

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

#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;
}

常见素数筛法

标签：io   os   sp   for   on   数据   ef   amp   as   

原文地址：http://blog.csdn.net/u013351484/article/details/41803017