for i=1to 500k=1:for j=1to k:k=k+j*!mod(i,j)next:if i