Прокомментируйте оставшиеся строчки, не понимаю как работает "Решето Эратосфена"

#include <iostream> //подключение стандартной библиотеки ввода-вывода
using namespace std; //пространство имён. Чтобы писать "cout << endl;" вместо "std::cout << endl;"
static const int N = 1000; //глобальная переменная. В данном случае (слово const) это константа, то есть изменить её значение невозможно. Сразу записываем её значение (1000)
int main()//главная функция. Здесь стартует программа.
{
int i, a[N]; //переменная i целого типа и целочисленный массив размером 1000
for (i = 2; i < N; i++) a[i] = 1; //Заполняем в цикле элементы массива a (a[i]), i=2,3...998,999  числом 1.
for (i = 2; i < N; i++){//Цикл от 2 до N-1(см.выше)
if (a[i])//Если элемент "i" массива "a" не является нулём, то:
for (int j = i; j*i < N; j++) a[i*j] = 0; //цикл по 'j'. Начальное значение j = i. Условие выполнения: пока j*i будут меньше N. После каждого шага цикла увеличиваем переменную j на один.
}
for (i = 2; i < N; i++) // от 2 до N-1
if (a[i]) cout << " " << i; //если элемент 'i' массива "а" отличен от нуля - выводим пробел и индекс этого элемента.
cout << endl;//перенос строки с очисткой буфера.
return 0; //выход из программы
}//конец главной функции