JZOJ 1845——约数

---

先将1~n所有非素数除1外的最小公因数求出来
最后循环求出加上每一个i的除1外的最小公约数，如果为0，则加本身

---

代码如下:

var  n,i,j:longint;  ans:int64;  w:array[2..10000000] of longint;begin  readln(n);  for i:=2 to trunc(sqrt(n)) do    if w[i]=0 then for j:=1 to n div i do if w[i*j]=0 then w[i*j]:=i;  for i:=2 to n do if w[i]>0 then ans:=ans+w[i] else ans:=ans+i;  writeln(ans);end.