Cod sursa(job #1592883)

Utilizator AeroHHorea Stefan AeroH Data 8 februarie 2016 02:16:50
Problema Puteri Scor 0
Compilator cpp Status done
Runda Arhiva de probleme Marime 1.31 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],m[135];

long long rasp;

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

        dp[m[x[i]]][m[y[i]]][m[z[i]]]++;
    }
}

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)
            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;
}