算法基礎-求1000以內的質數

//篩選1000以內的所有質數
//說明：除了自身以外，無法被其它整數整除的數稱之為質數。
//
#include <stdio.h>
#include <stdlib.h>
#define N 1000

int main(void) {
    int i, j;
    int prime[N + 1];
    //初始化數組都等於1
    for(    i = 2; i <= N; i++){
        prime[i]=1;
    }
    //嵌套for循環篩選出能整除的相對於的數組至於0
    for (i = 2; i * i <= N; i++) {
        if (prime[i] == 1) {
            for (j = 2 * i; j <= N; j++) {
                if (j % i == 0)
                    prime[j] = 0;
            }
        }
    }
    //遍歷數組、數組還等於1的是質數打印出來
    for (i = 2; i < N; i++) {
        if (prime[i] == 1) {
            printf("%4d", i);
            //隔幾個換行
            if ((i % 16) == 0){//不知道為什么沒有起作用
                printf("\n");
            }
        }
    }
    printf("\n");
    return 0;
}

找一些簡單的數學算法敲一敲。