【数论】关于质数

质数判定

试除法

bool is_prime(int n)
{
	if(n<2)return false;
	
	
	for(int i=2;i<n/i;i++)
	{
		if(n%i==0)return false;
	}
	
	return true;
}

注：

for循环终止条件要注意，若写i<sqrt(n)，sqrt函数太慢了，若写i*i<n，i*i可能超过int的范围；

数的因子成双出现，所以若存在大于sqrt(n)的因子，也一定存在小于sqrt(n)的因子。所以循环枚举到i<=sqrt(n)即可

分解质因数

void divide(int x)
{
	for(int i=2;i<=n/i;i++)
		if(n%i==0)
		{
			int s=0;
			while(n%i==0)
			{
				n/=i;
				s++;
			}
			printf("%d %d\n",i,s);//输出质因子和质因子幂次
		}
	if(n>1)printf("%d %d\n",n,1); //单独处理大于sqrt(n)的质因子
}

埃式筛法

#include<iostream>
#include<cstdio>
#include<cstring>
#include<algorithm>
using namespace std;

const int N=100010;

int primes[N],cnt;
bool st[N];

void get_primes(int n)
{
	for(int i=2;i<=n;i++)
	{
		if(!st[i])
		{
			primes[cnt++]=n;
			for(int j=i+i;j<=n;j+=i)st[j]=true;
		}
	}
}


int main()
{
	int n;
	cin>>n;
	
	get_primes(n);
	cout<<cnt<<endl;
	
	return 0;
}

注：埃式筛法只用质数去筛选质数，例如数i只会遍历一遍小于i的质数，如果都不是i的因子，说明数i是质数。算法时间复杂度为O（nloglogn）；

线性筛

#include<iostream>
#include<cstdio>
#include<cstring>
#include<algorithm>
using namespace std;

const int N=100010;

int primes[N],cnt;
bool st[N];

void get_primes(int n)
{
	for(int i=2;i<=n;i++)
	{
		if(!st[i])primes[cnt++]=i;
		for(int j=0;primes[j]<=n/i;j++)
		{
			st[primes[j]*i]=true;
			if(i%primes[j]==0)break;
		}
	}
}


int main()
{
	int n;
	cin>>n;
	
	get_primes(n);
	cout<<cnt<<endl;
	
	return 0;
}

注：在线性筛法中，每个和数只会被最小质因子筛去，当数据范围为1e6时，与埃式筛法时速差不多，当数据范围为1e7时，比埃式筛法快一倍