Cod sursa(job #232870)

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

using namespace std; 

bitset<1000002> a;

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