素数(质数)指的是指大于一的自然数 的因子只有它本身和1的数,1,0不是素数也不是合数。
方法一:一个数 n如果是合数,那么它的所有的因子不超过sqrt(n)--n的开方,所以,判定单个数是否是素数通常是通过这个原理,求a是否为素数,做一个小于等于sqrt(a)的循环即可判定。复杂度是o(n*sqrt(n))
int num=0;
for(int i=2;i<=n;i++)
{
for(int j=2;j*j<=i;j++)
{
if(i%j==0)
break;
if(j*j>i)
{
prime[num++]=i;
}
}
}
方法二: 埃拉托色尼(Eratosthenes)筛法
这种筛法在N为10^6甚至更大时可以很快速的求出1到N之间的所有素数,把之前找到的素数的倍数标记为合数,没被标记的自然为素数,时间复杂度o(nloglog(n)),这种方法在解决一般的素数问题中足矣。
bool prime[N];
int p[N];
int k=0;
void isprime(int n)
{
long long i,j;
memset(prime,true,sizeof(prime));
for(int i=2;i<=n;i++)
{
if(prime[i])
{
p[k++]=i;
{
for(int j=i*i;j<n;j+=i)
{
prime[j]=false;
}
}
}
}
}
方法三:线性筛法 它与埃拉托色尼筛法类似,但它保证使得任何一个合数,只被它最小的质因数标记过一次。所以整个算法是线性的。 时间复杂度为o(n),但就实际问题而言,埃拉托色尼筛法和线性筛法在运算速度上并没有太大的区别
bool prime[N];
int p[N];
void Init()
{
int cnt = 0;
memset(prime,true,sizeof(prime));
for (int i = 2; i < N; i++)
{
if (prime[i])
p[cnt++] = i;
for (int j = 0; j < cnt && p[j] * i < N; j++)
{
prime[i*p[j]] = false;
if (i % p[j] == 0)
break;
}
}
}