素数判断:
bool is_prime(int a)
{
for(int i=2;i*i<=a;i++)
if(a%i==0)
return false;
return true;
}
素数筛
两部分构成:埃氏筛+欧拉筛
- 埃氏筛:在给定范围内,找出n的倍数的所有的数并筛去,时间复杂度位O(nloglogn)。
void ai_prime(int m)
{
for(int i=2;i<=m;i++)
{
if(!visit[i]){
cnt++;
prime[cnt]=i;
}
for(int j=2;j*i<=n;j++)
{
visit[j*i]=1;
}
}
}
-
欧拉筛:让每个合数减少重复筛选的次数,时间复杂度O(n)。
eg:2与4的公倍数会被重复筛选
void ou_prime(int m)
{
for(int i=2;i<=m;i++)
{
if(!visit[i]){
cnt++;
prime[cnt]=i;
}
for(int j=1;j<=cnt;j++)//从已有的素数中筛选(判断)
{
if(i*prime[j]>m) break;//超出给定范围
visit[i*prime[j]]=1;//乘以相应的倍数
if(i%prime[j]==0) break;//过滤点重复的筛选
}
}
}
有空看看:
数论部分——Matrix67
米勒罗宾素数检测法——smiling~
https://blog.csdn.net/dy416524/article/details/86431057