ll n,s,t,z,a[2d5],b[1d3],c[3d5];

{
	rd(n,a(n));
	rep[a](x,n){
		rep(k,2,1d3){
			b[k]+=x%k;
		}
		c[x]+=1;
		s+=x;
	}
	c[1..2d5]+=c[0..];
	rep[a](x,n){
		if(x<1d3){
			z+=b[x];
		}else{
			z+=s;
			t=2d5/x*x;
			while(t){
				z-=x*(c[2d5]-c[t-1]);
				t-=x;
			}
		}
	}
	wt(z);
}