Pagini recente » Cod sursa (job #2482884) | Cod sursa (job #1615064) | Cod sursa (job #1396359) | Cod sursa (job #760741) | Cod sursa (job #2488053)
#include <fstream>
#include <string.h>
#define DIM 1000000
using namespace std;
ifstream fin ("fractii.in");
ofstream fout("fractii.out");
long long sol;
int x,i,j,c,n,k,y,v[1000010],p[80000];
int main () {
fin>>n;
for(i=2;i<=n;i++)
if(v[i]==0){
p[++k]=i;
v[i]=i-1;
sol+=v[i];
for(j=i+i;j<=n;j+=i)
v[j]=1;
}else{
x=i;
for(j=1;j<=k;j++){
if(x%p[j]==0 && x/p[j]>1){
y=x;
c=0;
while(y%p[j]==0)
y/=p[j],c++;
if(y==1){
v[i]=x/p[j]*(p[j]-1);
sol+=v[i];
break;
}
v[i]=(p[j]-1)*v[y]*(x/y/p[j]);
sol+=v[i];
break;
}
}
}
fout<<sol*2+1;
}