Cod sursa(job #2229071)

Utilizator inquisitorAnders inquisitor Data 5 august 2018 19:57:04
Problema Ciurul lui Eratosthenes Scor 100
Compilator cpp Status done
Runda Arhiva educationala Marime 1.58 kb
#include <stdio.h>

int p[2000000];

int sieve(int n)
{
    int i, j, nr = (n > 1) + (n > 2);

    for(i = 5; i * i <= n; i += 4)
    {
        if((p[i >> 5] & (1 << (i & 31))) == 0)
        {
            j = i * i;

            if((i - 5) % 3)
            {
                p[j >> 5] |= (1 << (j & 31));

                j += 4 * i;
            }

            for(; j <= n; j += 4 * i)
            {
                p[j >> 5] |= (1 << (j & 31));

                if((j += 2 * i) <= n) p[j >> 5] |= (1 << (j & 31));
            }
        }

        i += 2;

        if(i * i <= n)
        {
            if((p[i >> 5] & (1 << (i & 31))) == 0)
            {
                j = i * i;

                if((i - 5) % 3)
                {
                    p[j >> 5] |= (1 << (j & 31));

                    j += 4 * i;
                }

                for(; j <= n; j += 4 * i)
                {
                    p[j >> 5] |= (1 << (j & 31));

                    if((j += 2 * i) <= n) p[j >> 5] |= (1 << (j & 31));
                }
            }
        }
    }

    for(i = 5; i <= n; i += 4)
    {
        if((p[i >> 5] & (1 << (i & 31))) == 0)
        {
            ++nr;
        }

        i += 2;

        if(i <= n)
        {
            if((p[i >> 5] & (1 << (i & 31))) == 0)
            {
                ++nr;
            }
        }
    }

    return nr;
}

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

    int N; scanf("%d", &N);

    printf("%d ", sieve(N));

    return 0;
}