Pagini recente » Cod sursa (job #846841) | Cod sursa (job #2592353) | Cod sursa (job #1665124) | Cod sursa (job #391335) | Cod sursa (job #761633)
Cod sursa(job #761633)
#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;
}