线性筛素数法

原创 2015年11月21日 19:42:41

线性筛素数法代码:

#include<iostream>
using namespace std;    
const long N = 200000;   
long prime[N] = {0},num_prime = 0;    
int isNotPrime[N] = {1, 1};   
int main()    
{     
     	for(long i = 2 ; i < N ; i ++)       
       	{            
		if(! isNotPrime[i])               
	 		prime[num_prime ++]=i;  
		//关键处1        
		for(long j = 0 ; j < num_prime && i * prime[j] <  N ; j ++)
    		{               
		      	isNotPrime[i * prime[j]] = 1;  
	  		if( !(i % prime[j] ) )  //关键处2                  
				break;           
		}        
	}        
	return 0;   
}  

思想就是不断i乘上一个质数将其筛掉。


稍作修改,可以舍弃求模运算,并且可以只考虑奇数,并且可以同时求出质因子中的最小质数与最大质数。当然下面的maxprime数组对于求质数的过程没有用,可以删掉。

也可以在maxprime数组处稍作变动,同时求得一个数是不是唯一的质因子的幂次方。时间复杂度依然还是原先的O(n),并且没有求模运算。


const int MaxPrimeNum = 78499, MaxSize = 1001, MaxSize1 = 31251;
int MinPrime[MaxSize], MaxPrime[MaxSize];		//分别记录可以整除MaxSize的最小与最大prime的index      //MaxPrime可以改成bit类型数组来记录是否是prime或者是否是Prime幂次方
int prime[MaxPrimeNum];


//自己修改,欧拉筛法求n以内的所有素数,存在prime数组里,index从1到num_prime-1 。
//并同时求出n以内的所有数的最小质因子与最大质因子,分别存在MaxPrime数组与MinPrime数组里
//复杂度O(n)
void GetPrimeSolution::euler_sieve(int n)
{
	int i, j, len, num_prime = 1, tmpMaxIndex,    tmp;
	i = 2;
	prime[len = tmpMaxIndex = MaxPrime[i] = MinPrime[i] = num_prime++] = i;

	for (i = 3; i < n; i += 2)
	{
		if (!MinPrime[i])
		{
			prime[len = tmpMaxIndex = MaxPrime[i] = MinPrime[i] = num_prime++] = i;
		}
		else
		{
			tmpMaxIndex = MaxPrime[i];
			len = MinPrime[i];
		}     

		for (j = 2; j <= len && (tmp = i * prime[j]) < n; ++j)
		{
			MaxPrime[tmp] = tmpMaxIndex;
			MinPrime[tmp] = j;
		}
	}

/*	for (i = 0; i < num_prime; ++i)
	{
		printf("%d\t", prime[i]);
	}

	for (i = 0; i < n; ++i)
	{
		//printf("%d\t%d\t%d\n", i, MinPrime[i], MaxPrime[i]);
		printf("%d\t%d\t%d\n", i, prime[MinPrime[i]], prime[MaxPrime[i]]);
	}*/
}



const int MaxPrimeNum = 78499, MaxSize = 1001;
int MinPrime[MaxSize], MaxPrime[MaxSize];		//分别记录可以整除MaxSize的最小与最大prime的index      //MaxPrime可以改成bit类型数组来记录是否是prime或者是否是Prime幂次方
int prime[MaxPrimeNum];


//自己修改,欧拉筛法求n以内的所有素数,存在prime数组里,index从1到num_prime-1 。
//并同时求出n以内的所有数的最小质因子与最大质因子,分别存在MaxPrime数组与MinPrime数组里
//复杂度O(n)
void euler_sieve(int n)
{
	int i, j, len, num_prime = 1, tmp;
	i = 2;
	prime[len = MinPrime[i] = num_prime++] = i;

	for (i = 3; i < n; i += 2)
	{
		if (!MinPrime[i])
		{
			prime[len = MinPrime[i] = num_prime++] = i;
		}
		else
		{
			len = MinPrime[i];
		}

		for (j = 2; j <= len && (tmp = i * prime[j]) < n; ++j)
		{
			MinPrime[tmp] = j;
		}
	}

	for (i = 0; i < num_prime; ++i)
	{
	printf("%d\t", prime[i]);
	}
/*
	for (i = 0; i < n; ++i)
	{
	//printf("%d\t%d\t%d\n", i, MinPrime[i], MaxPrime[i]);
	printf("%d\t%d\t%d\n", i, prime[MinPrime[i]], prime[MaxPrime[i]]);
	}*/
}

int main()
{
	euler_sieve(100);
	return 0;
}
















版权声明:本文为博主原创文章,未经博主允许不得转载。

当我真正理解素数线性筛法

参考自:点击链接 主要代码:const int MAXN = 10000010; bool com[MAXN]; int primes, prime[MAXN/10];//数组不必开的太大 ...

普通筛素数及线性筛素数

一般筛法:一般筛法: const int MAXN = 70000; void Prime() // 筛法求素数表p[i]=1为素数,类似于hash映射! { int i, j; ...
  • leolin_
  • leolin_
  • 2011年07月28日 20:47
  • 7725

素数线性筛模板与证明

#include const int TOP=1e7+10000;//+10000是为了多筛一个素数 bool e[TOP]; int p[TOP/5]; int pnum; void prime()...

线性筛

#include #include #include #include #include using namespace std; const int maxn = 1e7 + ...
  • SM_545
  • SM_545
  • 2017年10月11日 21:50
  • 128

线性筛素数

给定一个范围N,你需要处理M个某数字是否为质数的询问(每个数字均在范围1-N内)i%prime[j]==0是O(n)筛素数的关键 为什么这里要break?先可以得出一个结论,此时的prime[j]为...

线性筛素数

线性筛素数的总结。

欧拉线性筛法求素数(顺便实现欧拉函数的求值)

我们先来看一下最经典的埃拉特斯特尼筛法。时间复杂度为O(n loglog n) int ans[MAXN]; void Prime(int n) { int cnt=0; memset(prime...
  • NK_test
  • NK_test
  • 2015年05月29日 23:02
  • 7061

素数筛法 埃氏筛 线性筛

素数筛法? 呃~,太暴力的试1~sqrt(n)时候是因数就不说了。 这里介绍埃氏筛和线性筛。埃氏筛埃氏筛主要原理是: 每次搜索到一个质数将他的倍数标记为合数; 本质原理就是:一个合数一定表示为...
  • sdfzsyh
  • sdfzsyh
  • 2017年07月06日 22:10
  • 211

素数线性筛

ps:证明转自神牛博客。代码:#include #includeusing namespace std; const int N=100100; int v[N],p[N],n,m,tot;int m...

3813: 奇数国 线段树+线性筛素数+欧拉函数+乘法逆元+压位

感觉就是几个算法的叠加。。没什么思维的难度。。无脑码代码呗。 首先维护区间乘法线段树就好,然后那个存在ax+by=1的限制条件,根据裴蜀定理其实据说gcd(x,y)=1即x,y互质,那么这个可以用欧...
内容举报
返回顶部
收藏助手
不良信息举报
您举报文章:线性筛素数法
举报原因:
原因补充:

(最多只允许输入30个字)