Pagini recente » Cod sursa (job #1378850) | Cod sursa (job #235919) | Cod sursa (job #767397) | Cod sursa (job #1129626) | Cod sursa (job #11019)
Cod sursa(job #11019)
Program cel;
Type vector=array[2..1000000] of boolean;
vector1=array[2..1000000] of longint;
var f,g:text;
a:vector;
phi:vector1;
n,i,j,p,nr,fm,d,x:longint;
Begin
Assign(f,'fractii.in');Reset(f);
Assign(g,'fractii.out');Rewrite(g);
Read(f,n);
fillchar(a,sizeof(a),true);
For i:=2 to trunc(sqrt(n)) do
If a[i]=true then Begin
j:=2;
While i*j<=n do
Begin
a[i*j]:=false;
inc(j);
end;
end;
For i:=2 to n do
If a[i]=true then phi[i]:=i-1;
i:=2;
While i<=trunc(sqrt(n)) do
Begin
p:=i;
While i*p<=n do
Begin
phi[i*p]:=((i*p) div p)*phi[p];
a[i*p]:=true;
i:=i*p;
end;
i:=p+1;
end;
For i:=2 to n do
If a[i]=false then Begin
d:=2;
phi[i]:=1;
x:=i;
repeat
fm:=0;
While x mod d=0 do
Begin
x:=x div d;
inc(fm);
end;
If fm<>0 then phi[i]:=phi[i]*phi[round(exp(fm*ln(d)))];
inc(d);
until x=1;
end;
nr:=0;
For i:=2 to n do nr:=nr+phi[i];
nr:=1+2*nr;
Write(g,nr);
Close(f);
Close(g);
end.