Cod sursa(job #2102565)

Utilizator inquisitorAnders inquisitor Data 9 ianuarie 2018 00:14:28
Problema Ciurul lui Eratosthenes Scor 30
Compilator cpp Status done
Runda Arhiva educationala Marime 0.71 kb
#include <cstdio>

using namespace std;

unsigned int N, p = 1;

int v[(2000000 >> 6) + 10];

int main()
{
    freopen("ciur.in", "r", stdin);
    freopen("ciur.out", "w", stdout);

    scanf("%d", &N);

    unsigned int i, j;

    for(i = 1; ((i * i) << 1) + (i << 1) <= N; i += 1)
    {
        if((v[i >> 3] & (1 << (i & 7))) == 0)
        {
            for(j = ((i * i) << 1) + (i << 1); (j << 1) + 1 <= N; j += (i << 1) + 1)
            {
                v[j >> 3] |= (1 << (j & 7));
            }
        }
    }

    for(i = 1; 2 * i + 1 <= N; ++i)
    {
        if((v[i >> 3] & (1 << (i & 7))) == 0)
        {
            ++p;
        }
    }

    printf("%d", p);

    return 0;
}