数组与素数筛法

数组与预处理命令

有关数组知识的补充

• 不论是几维数组，在计算机中的存储都是连续的。
• 函数内部在栈区开辟的数组里面的值不一定是“干净的”即不都为0，因此最好是初始化一下：int a[2] = {0},，而int b[2] = {1},则是int b[0] = 1,而已，要想全赋值为相同的数的话就得用memset，初始化为0的那样一种方式是一个特例。
• 函数外部在堆区开辟的数组里面一定是“干净的”。

素数筛

1. 标记(利用数组)一个范围内的数字是否是合数，没有被标记的则为素数
2. 算法的空间复杂度为O(n)，时间复杂度为O(n * loglogn)(不是真正的O(n)，因为标记过程中会重复标记)
3. 总体思想是是用素数去标记掉不是素数的数字，例如我已经知道了i是素数，那么2 * i、3 * i、4 * i…就都不是素数
4. 约数个数定理和约束和定理：任何一个正整数都可以写成素因子连乘的形式。

#include <stdio.h>
#define MAX_N 100

int prime[MAX_N + 5];
void init_prime() {
    for (int i = 2; i <= MAX_N; i++) {
        if (prime[i] == 0) {
            for (int j = 2 * i; j <= MAX_N; j += i) {
                prime[j] = 1;
            }
        }
    }
}

int main() {
    init_prime();
    for (int i = 2; i <= MAX_N; i++) {
        if (prime[i] == 0) {
            printf("%d\n", i);
        }
    }
    return 0;
}

根据对偶逻辑优化代码(减少缩进)：

#include <stdio.h>
#define MAX_N 100

int prime[MAX_N + 5];
void init_prime() {
    for (int i = 2; i <= MAX_N; i++) {
        if (prime[i]) continue;
        for (int j = 2 * i; j <= MAX_N; j += i) {
            prime[j] = 1;
        }
    }
}

int main() {
    init_prime();
    for (int i = 2; i <= MAX_N; i++) {
        if (prime[i]) continue;
        printf("%d\n", i);
    }
    return 0;
}

将找到的素数顺序存储到一个数组中而且不用再去开辟空间：

#include <stdio.h>
#define MAX_N 100

int prime[MAX_N + 5];
void init_prime() {
    for (int i = 2; i <= MAX_N; i++) {
        if (prime[i]) continue;
        prime[++prime[0]] = i;//改动的部分
        for (int j = 2 * i; j <= MAX_N; j += i) {
            prime[j] = 1;
        }
    }
}

int main() {
    init_prime();
    for (int i = 1; i <= prime[0]; i++) {
        printf("%d\n", prime[i]);
    }
    return 0;
}

• prime[0]可以当作计数量，统计素数的个数。在全局区下开辟的空间数组内都是0，因此prime[0] = 0.
• 标记比存的速度快
• 减少重新标记：for (int j = i * i; j <= MAX_N; j += i)。但这种方法有bug：令MAX_N = 100000 → segment fault，i * i超出int范围了
  解决方案：从数据的表示着手达到优化的目的。

#include <stdio.h>
#define MAX_N 100000

int prime[MAX_N + 5];
void init_prime() {
    for (int i = 2; i <= MAX_N; i++) {
        if (prime[i]) continue;
        prime[++prime[0]] = i;
        for (int j = i; j <= MAX_N / i; j++) {
            prime[j * i] = 1;
        }
    }
}

int main() {
    init_prime();
    printf("%d\n", prime[0]);
    return 0;
}

• 第9行也可以改成for (int j = i, int I = MAX_N / i; j <= I; j++)这样改以后MAX_N / i在该层循环中执行一次，减少执行的次数