Cod sursa(job #514479)

Utilizator ShootMeBistriceanu Andrei ShootMe Data 18 decembrie 2010 19:55:17
Problema Fractii Scor 30
Compilator cpp Status done
Runda Arhiva de probleme Marime 1.19 kb
#include <stdio.h>

#define MAXN 1000001

int n;
int ciur[MAXN];

int max(int a, int b) {
	return a < b ? b : a;
}

void init_ciur(int val) {
	for (int i = 0; i <= n; i++) {
		ciur[i] = val;
	}
}

void make_fi_of_n_ciur() {
	for (int i = 2; i <= n; i++) {
		if (ciur[i] & 1) {
			int max_steps = (n - i) / i;
			int pos = i << 1;
			ciur[i] = i - 1;
			for (int steps = 0; steps < max_steps; steps++, pos += i) {
				int pow = 1;
				int aux = pos;
				while (aux % i == 0) {
					pow *= i;
					aux /= i;
				}
				ciur[pos] *= (i-1)*(pow / i);
			}
		}
	}
}

//int cmmdc(int a, int b) {
//	if (a % b != 0) {
//		return cmmdc(b, a % b);
//	} else {
//		return b == 1 ? 0 : b;
//	}
//}
//
//void brute() {
//	int sum = 1;
//	for (int i = 1; i <= n; i++) {
//		for (int j = i + 1; j <= n; j++) {
//			if (cmmdc(i, j) == 0) {
//				sum += 2;
//			}
//		}
//	}
//	printf("\n\n\n%d\n\n\n", sum);
//}

int main() {
	freopen("fractii.in", "r", stdin);
	freopen("fractii.out", "w", stdout);

	scanf("%d", &n);

	init_ciur(1);

	make_fi_of_n_ciur();

	int sum = 1;
	for (int i = 2; i <= n; i++) {
		sum += 2*ciur[i];
	}
	printf("%d\n", sum);

	return 0;
}