Pagini recente » Cod sursa (job #1870215) | Cod sursa (job #566888) | Cod sursa (job #1448919) | Cod sursa (job #2401305) | Cod sursa (job #1693168)
#include <fstream>
#include <vector>
#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;
P.push_back(2);
for (i=3; i<=n; i+=2)
{
if (Prim[i]==0)
{
P.push_back(i);
for (j=i; j<=n; j+=2*i)
Prim[j]=1;
}
}
}
long long Euler(int x)
{
long long a=x, i;
for (i=0; P[i]<=x && i<P.size(); i++)
{
if (x%P[i]==0)
{
a/=P[i];
a*=P[i]-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;
}