Cod sursa(job #761633)

Utilizator warchildmdMihail Burduja warchildmd Data 26 iunie 2012 21:08:16
Problema Fractii Scor 80
Compilator cpp Status done
Runda Arhiva de probleme Marime 2.25 kb
#include <cstdio>
#include <cmath>

long long SUM;
int N;
int primes[1000001];
int p[100000];
int fis[1000001];
int nr_primes;

int fi(int n)
{
    if (primes[n] == 0)
    {
        return n - 1;
    }
    int result = n;
    /*if (n % 2 == 0) result -= result / 2;
    while (n % 2 == 0) n /= 2;
    for(int i=3; i*i <= n; i+=2)
    {
        if (n % i == 0) result -= result / i;
        while (n % i == 0) n /= i;
    }*/
    for(int i = 0; p[i] * p[i] <= n; i++)
    {
        if (n % p[i] == 0) result -= result / p[i];
        while (n % p[i] == 0) n /= p[i];
    }
    if (n > 1) result -= result / n;
    return result;
}

void check_primes()
{
    p[nr_primes] = 2;
    nr_primes++;

    for(int i = 4; i <= N; i+=2)
    {
        primes[i] = 1;
    }

    for(int i = 3; i <= N; i+=2)
    {
        if(primes[i] == 0)
        {
            p[nr_primes] = i;
            nr_primes++;

            for(int j = 2*i; j <= N; j+=i)
            {
                primes[j] = 1;
            }
        }
    }
}

int pow(int x, int y)
{
    if (y == 0)
    {
        return 1;
    }
    if (y == 1)
    {
        return x;
    }
    if (y % 2 == 0)
    {
        return pow(x, y / 2) * pow(x, y / 2);
    }
    return pow(x, y / 2) * pow(x, y / 2) * x;
}

int euler(int x)
{
    if (x == 1)
    {
        return 1;
    }
    if (primes[x] == 0)
    {
        return x - 1;
    }
    int PRODUCT = 1;
    int current = 0;
    while(x > 1)
    {
        if (x % p[current] == 0)
        {
            int count = 0;
            int total = 1;
            while (x % p[current] == 0)
            {
                total *= p[current];
                x /= p[current];
                count++;
            }
            //return fis[x / total] * fis[total];
            PRODUCT = PRODUCT * (p[current] - 1) * (total / p[current]);
        }
        current++;
    }
    return PRODUCT;
}

int main()
{
    freopen("fractii.in", "r", stdin);
    freopen("fractii.out", "w", stdout);
    scanf("%d", &N);
    check_primes();

    for (int i = 2; i <= N; i++)
    {
        fis[i] = fi(i);
        SUM += fis[i];
    }

    SUM = 2 * SUM + 1;

    printf("%lld\n", SUM);
    return 0;
}