Cod sursa(job #1561066)

Utilizator jurjstyleJurj Andrei jurjstyle Data 3 ianuarie 2016 17:09:11
Problema Ciurul lui Eratosthenes Scor 100
Compilator cpp Status done
Runda Arhiva educationala Marime 0.8 kb
/*
 * Nu ne intereseaza numerele pare
 * pentru ca 2 este singurul numar prim par
 * si il tratam separat
 * Pentru a tine minte ca pe pozitia i nu este un numar prim,
 * marcam in vector bitul [i / 2]
 * Daca in vector este marcat bitul i, atunci inseamna ca
 * numarul 2*i+1 nu este prim
 */

#include <cstdio>

const int NMAX = 1 << 17;

unsigned char V[NMAX];

int main(void) {
    freopen("ciur.in", "rt", stdin);
    freopen("ciur.out", "wt", stdout);
    int N, i, j, i2, nrprime = 0;

    scanf(" %d", &N);

    for (i = 3; i <= N; i += 2) {
        if (V[i >> 4] & (1 << ((i >> 1) & 7))) continue;
        ++nrprime;

        for (j = i + (i2 = i + i); j <= N; j += i2)
            V[j >> 4] |= 1 << ((j >> 1) & 7);
    }

    printf("%d\n", nrprime + 1);

    return 0;
}