# PBinfo #3466 – num_sum_div – Rezolvare C++

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

## Enunț

Cerința
Pentru un număr natural x notăm cu S suma divizorilor săi diferiți de x. Dacă S este strict mai mică decât x, atunci x se numește număr deficient, dacă S este egală cu x, atunci x se numește număr perfect, iar dacă S este strict mai mare decât x, atunci x se numește număr abundent.

Se dă un șir de n numere naturale. Să se calculeze câte numere sunt deficiente, perfecte, respectiv abundente.

Date de intrare
Fișierul de intrare num_sum_div.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 num_sum_div.out va conține pe prima linie 3 numere separate prin câte un spațiu, reprezentând câți termeni ai șirului sunt numere deficiente, perfecte, respectiv abundente.

Restricții și precizări

1 ≤ n ≤ 10.000
numerele de pe a doua linie a fișierului de intrare vor fi mai mici decât 5.000.000.000
punctajul pe test va fi acordat doar dacă toate cele 3 numere afișate sunt corecte

Exemplu:
num_sum_div.in

10
2 6 8 20 21 13 14 15 12 24
num_sum_div.out

6 1 3
Explicație
În fișierul de intrare sunt 6 numere deficiente, un singur număr perfect și 3 numere abundente.

## Cod sursă C++

```cpp
#include <bits/stdc++.h>
using namespace std;
class InParser
{
private:
    FILE *fin;
    char *buff;
    int sp;
    char read_ch()
    {
        ++sp;
        if(sp==4096)
        {
            sp=0;
            fread(buff,1,4096,fin);
        }
        return buff[sp];
    }
public:
    InParser(const char* nume)
    {
        fin=fopen(nume,"r");
        buff=new char[4096]();
        sp=4095;
    }
    InParser& operator >> (int &n)
    {
        char c;
        while(!isdigit(c=read_ch())&&c!='-');
        int sgn=1;
        if (c=='-')
        {
            n=0;
            sgn=-1;
        }
        else
        {
            n=c-'0';
        }
        while(isdigit(c=read_ch()))
        {
            n=10*n+c-'0';
        }
        n*=sgn;
        return *this;
    }
    InParser& operator >> (long long &n)
    {
        char c;
        n=0;
        while(!isdigit(c=read_ch())&&c!='-');
        long long sgn=1;
        if(c=='-')
        {
            n=0;
            sgn=-1;
        }
        else
        {
            n=c-'0';
        }
        while(isdigit(c=read_ch()))
        {
            n=10*n+c-'0';
        }
        n*=sgn;
        return *this;
    }
};
class OutParser
{
private:
    FILE *fout;
    char *buff;
    int sp;
    void write_ch(char ch)
    {
        if(sp==50000)
        {
            fwrite(buff,1,50000,fout);
            sp=0;
            buff[sp++]=ch;
        }
        else
        {
            buff[sp++]=ch;
        }
    }
public:
    OutParser(const char* name)
    {
        fout=fopen(name,"w");
        buff=new char[50000]();
        sp=0;
    }
    ~OutParser()
    {
        fwrite(buff,1,sp,fout);
        fclose(fout);
    }
    OutParser& operator <<(int vu32)
    {
        if(vu32<=9)
        {
            write_ch(vu32+'0');
        }
        else
        {
            (*this) <<(vu32/10);
            write_ch(vu32%10+'0');
        }
        return *this;
    }
    OutParser& operator <<(long long vu64)
    {
        if(vu64<=9)
        {
            write_ch(vu64+'0');
        }
        else
        {
            (*this) <<(vu64/10);
            write_ch(vu64%10+'0');
        }
        return *this;
    }
    OutParser& operator <<(char ch)
    {
        write_ch(ch);
        return *this;
    }
    OutParser& operator <<(const char *ch)
    {
        while(*ch)
        {
            write_ch(*ch);
            ++ch;
        }
        return *this;
    }
};
InParser f ("num_sum_div.in");
OutParser g ("num_sum_div.out");
const int NN = 70715;
bool ciur[NN + 1];
int prime[32768], np, k1, k2, k3;
int main ()
{
    for (int i = 2; i * i <= NN; i++)
    if (ciur[i] == 0)
    for (int j = 2; i * j <= NN; j++)
    ciur[i * j] = 1;
    prime[++np] = 2;
    for (int i = 3; i <= NN; i += 2)
    if (ciur[i] == 0)
    prime[++np] = i;
    long long n, x, m, S;
    f >> n;
    while (n--)
    {
        f >> x;
        m = x, S = 1;
        int d = 1;
        while (m > 1)
        {
            long long f = 1, t = 1;
            while (m % prime[d] == 0)
                t += (f *= prime[d]), m /= prime[d];
            S *= t;
            d++;
            if (1LL * prime[d] * prime[d] > m && m > 1)
                    S *= (m + 1), m = 1;
        }
        if (S < 2 * x) k1++;
        else if (S == 2 * x) k2++;
        else if (S > 2 * x) k3++;
    }
    g << k1 << " " << k2 << " " << k3;
    return 0;
}
```

## Resurse

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