Cod sursa(job #1592882)

Utilizator AeroHHorea Stefan AeroH Data 8 februarie 2016 02:10:26
Problema Puteri Scor 10
Compilator cpp Status done
Runda Arhiva de probleme Marime 1.3 kb
#include <fstream>
#include <cstring>
#define ll long long int
#define punct pair<char,char>
using namespace std;

ifstream f("puteri.in");
ofstream g("puteri.out");

int n,x[100001],y[100001],z[100001],t1,t2,t3,dp[65][65][65],i,j,v[135],pr[135];

long long rasp;

void check(int mod)
{
    int i,a,b,c;
    memset(dp,0,sizeof(dp));
    for(i=1;i<=n;++i)
    {
        a = mod - x[i]%mod;
        b = mod - y[i]%mod;
        c = mod - z[i]%mod;
        if(a%mod <= 64 && b%mod <= 64 && c%mod <= 64)
        v[mod]+=dp[a%mod][b%mod][c%mod];

        dp[x[i]%mod][y[i]%mod][z[i]%mod]++;
    }
}

bool isp(int t)
{
    for(int i=2;i*i<=t;++i)
        if (t%i == 0)
            return 0;
    return 1;
}

int main()
{
    f>>n;
    for(i=1;i<=n;++i)
        f>>x[i]>>y[i]>>z[i];

    for(i=2;i<=128;++i)
        if (isp(i))
            check(i);

    for(i=2;i<=128;++i)
        if (isp(i))
            rasp+=v[i];

    for(i=2;i<=128;++i)
        for(j=2;j<=128;++j)
            if(i!=j && isp(i) && isp(j))
                if (i*j<=128)
                    rasp-=v[i*j];
    rasp+=v[30];
    rasp+=v[42];
    rasp+=v[66];
    rasp+=v[78];
    rasp+=v[102];
    rasp+=v[114];
    rasp+=v[70];
    rasp+=v[110];
    rasp+=v[105];
    g<<rasp;

    return 0;
}