# PBinfo #3459 – CountPrime – Rezolvare C++

Sursa: Codul lui Andrei
Limbaj: C++
Autor rezolvare: deivid

## Enunț

Cerința
Să se calculeze cate numere prime sunt în intervalul [st, dr].

Date de intrare
Fișierul de intrare countprime.in conține pe prima linie două numere st și dr.

Date de ieșire
Fișierul de ieșire countprime.out va conține pe prima linie numărul cnt, reprezentând numărul de numere prime din intervalul dat.

Restricții și precizări

1 ≤ st ≤ dr < 232
dr - st ≤ 1.000.000

Exemplu:
countprime.in

2 10
countprime.out

4
Explicație
Sunt 4 numere prime în intervalul [2, 10].

## Cod sursă C++

```cpp
#include <fstream>
using namespace std;
ifstream fin("countprime.in");
ofstream fout("countprime.out");
typedef unsigned long long ull;
ull check[] = {2, 3, 5, 7, 11, 13, 17, 19, 23, 29, 31, 37, 41};
ull power(ull val, ull put, ull mod){
    ull rez = 1;
    ull pw = val;
    while(put){
        if(put & 1){
            rez *= pw;
            rez %= mod;
        }
        pw *= pw;
        pw %= mod;
        put >>= 1;
    }
    return rez;
}
bool prim(ull num, ull baze){
    ull cn = num - 1;
    ull exp = 0;
    while(cn % 2 == 0)
    {
        exp++;
        cn /= 2;
    }
    ull x = (num - 1) / ((ull)1 << exp);
    ull pw = power(baze, x, num);
    if(pw == 1 || pw == num - 1)
        return true;
    for(ull i = 1; i < exp; ++i){
        pw = power(pw, 2, num);
        if(pw == num - 1)
            return true;
    }
    return false;
}
int main()
{
    ull st, dr;
    fin >> st >> dr;
    ull cnt = 0;
    for(ull i = st; i <= dr; ++i){
        bool ok = 0;
        for(int k = 0; k < 13 && !ok; ++k)
            if(i == check[k])
                ++cnt, ok = 1;
        if(!ok){
            ok = 1;
            for(int k = 0; k < 13 && ok; ++k)
                ok = prim(i, check[k]);
            cnt += ok;
        }
    }
    fout << cnt << '\n';
    return 0;
}
```

## Resurse

- [Pagina originală](https://codulluiandrei.ro/rezolvari/pbinfo/3459)
- [Sursa C++](https://codulluiandrei.ro/rezolvari/pbinfo/3459/cod.cpp)
- [Sursa pe GitHub](https://raw.githubusercontent.com/codulluiandrei/pbinfo/main/pbinfo-3459/main.cpp)
- [Cursuri de informatică pentru BAC](https://codulluiandrei.ro/academia/info-curs/pachet-cpp/)