发布时间:2023-07-11 16:36:55
int find_primes (int max, int **arr)
{
#define MAX_SIZE max
bool is_prime[max];
int primes[max];
#undef MAX_SIZE
memset (is_prime, true, sizeof (is_prime));
is_prime[0] = is_prime[1] = false;
int cnt = 0;
for (int i = 2; i < max ; ++i) {
if (is_prime[i]) primes[cnt++] = i;
for (int j = 0; j < cnt && i * primes[j] < max; ++j) {
is_prime[i * primes[j]] = false;//每个质数的倍数都不是质数
if (i % primes[j] == 0) break;
//说明i不是质数,存在primes[j]可以整除i,
//此时无需再继续判断,因为i*primes[j]是合数,
//更大的质数*i也是合数,只需要随着i的增加
//就能用小的质数(primes[j])筛掉i对应的合数。
}
}
int *cloned = malloc (sizeof (int) * cnt);
memcpy (cloned, primes, sizeof (int) *cnt);
*arr = cloned;
return cnt;
}