CIUR

void ciur(int N, int s[])
{

    vector<bool> prime(N + 1, false);


    // pentru toate numerele pare
    for (int i = 2; i <= N; i += 2)
        s[i] = 2;

    // pentru numere impare <=N
    for (int i = 3; i <= N; i += 2) {
        if (prime[i] == false) {

            // s(i) pentru un număr prim este el însuși cel mai mic factor prim
            s[i] = i;
  		/// aici s-ar putea memora toate numerele prime într-un alt vector
            // pentru toți multipli lui
            for (int j = i; j * i <= N; j += 2) {
                if (prime[i * j] == false) {
                    prime[i * j] = true;

                    // i este cel mai mic factor prim
                    // number "i*j".
                    s[i * j] = i;
                }
            }
        }
    }
}