Pagini recente » Cod sursa (job #1765880) | Cod sursa (job #2483244) | Cod sursa (job #505873) | Cod sursa (job #2277955) | Cod sursa (job #1645933)
#include <bits/stdc++.h>
using namespace std;
typedef long long ll;
const int Pmax = 1e3;
int n , i , j;
int phi[Pmax+10];
unordered_map < int , ll > ans;
void bagaphi(int n)
{
for (i = 1; i <= n; ++i)
phi[i] = i;
for (i = 2; i <= n; ++i)
if (phi[i] == i)
for (j = i; j <= n; j += i) phi[j] = phi[j] / i * (i - 1);
for (i = 1; i <= n; ++i)
phi[i] += phi[i-1];
}
ll solve(int n)
{
if (n <= Pmax) return phi[n];
ll ret = ans[n];
if (ret) return ret;
ret = n * (n + 1) / 2;
int l , r , d;
for (l = 2; l <= n; l = r + 1)
{
d = n / l;
r = n / d;
ret -= 1LL * (r - l + 1) * solve(d);
}
ans[n] = ret;
return ret;
}
int main()
{
freopen("fractii.in","r",stdin);
freopen("fractii.out","w",stdout);
scanf("%d", &n);
bagaphi(min(Pmax , n));
printf("%lld\n", 2LL * solve(n) - 1LL);
return 0;
}