Cod sursa(job #232869)

Utilizator crawlerPuni Andrei Paul crawler Data 16 decembrie 2008 10:06:23
Problema Ciurul lui Eratosthenes Scor 100
Compilator cpp Status done
Runda Arhiva educationala Marime 0.58 kb
#include <stdio.h>
#include <math.h>
#include <bitset>

using namespace std; 

bitset<1000100> a;

int main()
{
    freopen("ciur.in","r",stdin);
    freopen("ciur.out","w",stdout);
    
    int n,ret=0,lim;
    
    scanf("%d",&n);
    if (n>=2) ++ret;
    lim = (int)sqrt(n);
        
    for (int i=3;i<=lim;i+=2) if(a[i/2] == 0)
    {
        ++ret;
        for (int j=i*i;j<=n;j+=2*i) a[j/2] = 1;        
    }
    ++lim;
    if (lim%2 == 0) ++lim;
    for (int i=lim;i<=n;i+=2) if(a[i/2] == 0) ++ret;
    
    printf("%d\n",ret);
        
    return 0;    
}