素数:只能被1和其本身整除
普通求一个数是不是素数,可能就是:方法1
bool JuNum(int t_v)
{
for (int i = 2; i <= sqrt(t_v + 0.5); (i==2)?(i++):(i+=2))
{
if (!(t_v%i)) return false;
}
return true;
}
下面介绍一种更高效的方法,一种素数打表法,它判断并统计小于一定值内的素数时可以说是线性的:方法2
#define N 100000
#define M 10000
bool isPrime[N+5];
int prime[M];
int cnt;
void JuPr()
{
int i, j;
memset(isPrime, 0, sizeof(isPrime));
cnt = 0;
for (i = 2; i <= N; i++)
{
if (!isPrime[i]) prime[cnt++] = i;
for (j = 0; j < cnt && (i*prime[j] <= N); j++)
{
isPrime[i*prime[j]] = 1;
//非素数都可以分解成两个数的乘积的形式或者若干个素数的乘积
//假如非素数n,那么n=A*B(A是素数,B是任意整数)成立
//那么有1,B=A^K(K是某一整数);2,A和B互质
//如果i%prime[j] == 0,那么i = prime[j] * J,且J>prime[j]
//所以一定存在i*prime[j+1] = I(>i) * prime[k]
if (!(i%prime[j])) break; //避免对一个数重复多次判定
}
}
}