Pagini recente » Cod sursa (job #3158332) | Monitorul de evaluare | Cod sursa (job #1150433) | Cod sursa (job #1766040) | Cod sursa (job #3304391)
#include <fstream>
using namespace std;
ifstream in("puteri.in");
ofstream out("puteri.out");
int n, ans;
short 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;
int 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 += 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;
}