Pagini recente » Cod sursa (job #2482132) | Cod sursa (job #3216329) | Cod sursa (job #3193782) | Cod sursa (job #841428) | Cod sursa (job #1693307)
#include <fstream>
#include <vector>
#include <cmath>
#define NM 1000005
using namespace std;
ifstream fin("fractii.in");
ofstream fout("fractii.out");
int n;
bool Prim[NM];
vector <int> P;
void GenCiur()
{
int i, j, lim;
P.push_back(2);
lim=sqrt(n);
for (i=3; i<=lim; i+=2)
{
if (Prim[i]==0)
{
P.push_back(i);
for (j=i; j<=lim; j+=2)
Prim[j]=1;
}
}
}
long long Euler(int x)
{
long long a=x, i;
int sz=P.size();
for (i=0; i<sz; i++)
{
if (x%P[i]==0)
{
a/=P[i];
a*=P[i]-1;
while (x%P[i]==0)
{
x /= P[i];
}
}
if (x==1)
break;
}
if (x!=1)
{
a/=x;
a*=x-1;
}
return a;
}
int main()
{
long long i, sol=0;
fin>>n;
GenCiur();
for (i=2; i<=n; i++)
{
sol+=Euler(i);
}
fout<<2*sol+1;
return 0;
}