Cod sursa(job #1592879)

Utilizator AeroHHorea Stefan AeroH Data 8 februarie 2016 02:01:01
Problema Puteri Scor 10
Compilator cpp Status done
Runda Arhiva de probleme Marime 1.16 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[64][64][64],i,j,v[65],pr[65];

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;
        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<=64;++i)
        if (isp(i))
            check(i);
    /*for(i=1;i<=64;++i)
        g<<v[i]<<" ";*/

    for(i=2;i<=64;++i)
        if (isp(i))
            rasp+=v[i];
    for(i=2;i<=64;++i)
        for(j=2;j<=64;++j)
            if(i!=j && isp(i) && isp(j))
                if (i*j<=64)
                    rasp-=v[i*j];
    rasp+=v[30];
    rasp+=v[42];

    g<<rasp;

    return 0;
}