Cod sursa(job #1801071)

Utilizator iulianrotaruRotaru Gheorghe-Iulian iulianrotaru Data 8 noiembrie 2016 17:12:40
Problema Puteri Scor 100
Compilator cpp Status done
Runda Arhiva de probleme Marime 0.93 kb
#include <bits/stdc++.h>
using namespace std;
ifstream f("puteri.in");
ofstream g("puteri.out");
bool prim[129];
int n,i,j,a[1<<17],b[1<<17],c[1<<17],bit[129],l[129],s[65][65][65];
long long sol;
inline long long solve(int mod)
{
    int i,a_,b_,c_;
    long long sol=0;
    memset(s,0,sizeof s);
    for(i=0;i<=128;i++) l[i]=i%mod;
    for(i=0;i<n;++i)
    {
        a_=l[mod-l[a[i]]];
        b_=l[mod-l[b[i]]];
        c_=l[mod-l[c[i]]];
        if(a_<=64&&b_<=64&&c_<=64) sol+=s[a_][b_][c_];
        s[l[a[i]]][l[b[i]]][l[c[i]]]++;
    }
    return sol;
}
int main()
{
    f>>n;
    for(i=0;i<n;i++) f>>a[i]>>b[i]>>c[i];
    for(i=2;i<=128;i++) bit[i]=1;
    for(i=2;i<=128;i++)
    {
        if(!prim[i])
        {
            for(j=i;j<=128;j+=i) bit[j]*=-1,prim[j]=1;
            for(j=i*i;j<=128;j+=i*i) bit[j]=0;
        }
        if(bit[i]!=0) sol-=bit[i]*solve(i);
    }
    g<<sol;
    return 0;
}