# PBinfo #3315 – Eratostene4 – Rezolvare C++

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

## Enunț

Cerința
Se dau n numere naturale. Pentru fiecare număr aflaţi câţi divizori liberi de pătrate are acesta.

Date de intrare
Fișierul de intrare eratostene4.in conține pe prima linie numărul n, iar pe a doua linie n numere naturale separate prin spații.

Date de ieșire
Fișierul de ieșire eratostene4.out va conține pe prima linie, pentru fiecare număr din fişierul de intrare, numărul divizorilor liberi de pătrate ai acestuia. 

Restricții și precizări

1 ≤ n ≤ 100.000
numerele de pe a doua linie a fișierului de intrare vor fi mai mici decât 10.000.000
un număr natural se numeşte liber de pătrate dacă nu se divide cu pătratul unui număr prim

Exemplu:
eratostene4.in

3
20 8 5
eratostene4.out

4 2 2
Explicație
Divizorii lui 20, liberi de pătrate, sunt: 1, 2, 5, 10.
Divizorii lui 8, liberi de pătrate, sunt: 1, 2. 
Divizorii lui 5, liberi de pătrate, sunt: 1, 5.

## Cod sursă C++

```cpp
#include <fstream>
#define N 10000000
using namespace std;
ifstream fi("eratostene4.in");
ofstream fo("eratostene4.out");
int i, j, n, b[21];
char p[N+1], a[N+1];
int main()
{
    //ciurul lui eratostene + numararea factorilor primi din descompunere
    p[1] = 1;
    for(i = 2; i <= N; i++)
        if(p[i]==0)
    {
        a[i] = 1;
        j = i+i;
        while(j <= N)
        {
            p[j] = 1;
            a[j]++;
            j = j + i;
        }
    }
    //puterile lui 2
    b[0] = 1;
    for(i = 1; i <= 20; i++)
        b[i] = b[i-1] * 2;
    //prelucrare sir
    fi >> n;
    for(i = 1; i <= n; i++)
    {
        fi >> j;
        fo << b[a[j]] << " ";
    }
    return 0;
}
```

## Resurse

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