Cod sursa(job #3304392)

Utilizator unomMirel Costel unom Data 23 iulie 2025 10:56:27
Problema Puteri Scor 10
Compilator cpp-64 Status done
Runda Arhiva de probleme Marime 1.47 kb
#include <fstream>

using namespace std;

ifstream in("puteri.in");
ofstream out("puteri.out");
int n;
long long ans;
int frec[65][64][64][64];
int mobius[105];

void build_mobius()
{
    mobius[1] = -1;

    for(int i = 1; i<=64; i++)
    {
        if(mobius[i])
        {
            mobius[i] = -mobius[i];

            for(int j = 2 * i; j<=64; j += i)
            {
                mobius[j] += mobius[i];
            }
        }
    }

    for(int i = 1; i<=64; i++)
    {
        mobius[i] *= -1;
    }
}

int main()
{
    in>>n;

    build_mobius();

//    for(int i = 1; i<=64; i++)
//    {
//        out<<mobius[i]<<" ";
//    }
//    out<<'\n';

    int a, b, c;
    for(int i = 1; i<=n; i++)
    {
        in>>a>>b>>c;

        long long rez = 0;
        for(int j = 2; j<=64; j++)
        {
            int ra = j - (a % j);
            int rb = j - (b % j);
            int rc = j - (c % j);

            if(ra == j)
            {
                ra = 0;
            }
            if(rb == j)
            {
                rb = 0;
            }
            if(rc == j)
            {
                rc = 0;
            }

            rez += 1LL * frec[j][ra][rb][rc] * mobius[j];

            //out<<rez<<" ";
        }
        //out<<'\n';

        ans += rez;

        for(int j = 2; j<=64; j++)
        {
            frec[j][a % j][b % j][c % j]++;
        }
    }

    out<<ans;

    return 0;
}