Pagini recente » Cod sursa (job #731206) | Cod sursa (job #1427931) | Cod sursa (job #688078) | Cod sursa (job #2846784) | Cod sursa (job #3146695)
#include <fstream>
#include <bitset>
using namespace std;
ifstream fin ("fractii.in");
ofstream fout ("fractii.out");
int k,j,n,i,d,x,v[500001];
long long sum,phi;
bitset <1000001> ciur;
int main()
{
fin>>n;
for (i=2; i<=n; i++)
{
if (ciur[i]==0)
{
v[++k]=i;
for (j=i+i; j<=n; j+=i)
ciur[j]=1;
}
}
for (i=2; i<=n; i++)
{
x=i;
phi=x;
for (d=1; d<=k&&v[d]*v[d]<=x; d++)
{
if (x%v[d]==0)
{
phi=phi/v[d]*(v[d]-1);
while (x%v[d]==0)
x/=v[d];
}
}
if (x!=1)
phi=phi/x*(x-1);
sum+=phi;
}
sum=sum*2+1;
fout<<sum;
return 0;
}