Cod sursa(job #1657413)

Utilizator BourucLiviuBouruc Petru Liviu BourucLiviu Data 20 martie 2016 14:38:12
Problema Ciurul lui Eratosthenes Scor 40
Compilator cpp Status done
Runda Arhiva educationala Marime 0.59 kb
#include <fstream>
using namespace std;

ifstream cin ("ciur.in");
ofstream cout("ciur.out");

int prim[1000001 / 31 + 1];

int main()
{
    int n, l = 0;
    cin >> n;
    for(int i = 2; i <= n; ++i) prim[i >> 5] |= (1 << (i & 31));
    for(int i = 4; i <= n; i += 2) prim[i >> 5] &=~ (1 << (i & 31));
    for(int i = 3; i * i <= n; i += 2)
        if(prim[i >> 5] & (1 << (i & 31)))
            for(int j = i * i; j <= n; j += (i << 1)) prim[j >> 5] &=~ (1 << (j & 31));
    for(int i = 2; i <= n; ++i)
        if(prim[i >> 5] & (1 << (i & 31))) l++;
    cout << l;
    return 0;
}