Cod sursa(job #2731813)

Utilizator AACthAirinei Andrei Cristian AACth Data 28 martie 2021 12:59:03
Problema Pairs Scor 0
Compilator cpp-64 Status done
Runda Arhiva de probleme Marime 1.29 kb
#include <bits/stdc++.h>
using namespace std;
ifstream f("pairs.in");
ofstream g("pairs.out");
#define int long long
const int Max = 1e6 + 1;
void nos()
{
    ios_base::sync_with_stdio(false);
    cin.tie(NULL);
}
int n;
vector < int > multime;
int maxx;
void read()
{
    f>>n;
    int i;
    for(i=1;i<=n;i++)
    {
        int x;
        f>>x;
        maxx = max(maxx,x);
        multime.push_back(x);
    }
}

int divisible[Max];
int dp[Max];
void baga(int n)
{
    divisible[1] ++;
    int d = 2;
    while(n)
    {
        if(n%d == 0)
        {
            divisible[d]++;
            while(n%d == 0)
                n /= d;
        }
        d++;
        if(d*d > n)
            break;
    }
    if(n!=1)
        divisible[n] ++;
}
void solve()
{
    for(auto it : multime)
        baga(it);
    //dp[i] nr de perechi care au cmmdc i
    //un fel de ciur invers?

    int i,j;
    for(i = maxx; i>= 1;i--)
    {
        int val = divisible[i];
        dp[i] = (val*(val-1)) / 2;
        for(j=2;j<=maxx / i; j++)
            dp[i] -= dp[i*j];
        //scadem multipli, principiul includerii excluderii
    }
    g<<dp[1];
}
void restart()
{

}

int32_t main()
{

        read();
        solve();
        restart();
    return 0;
}