X mod Y =X-[X/Y]
=>sum{k mod i,{i,1,k}}=sum{i,{i,1,k}}-sum{[i/n],{i,1,k}}
sum{i,{i,1,k}}=k*(k+1)/2
sum{[i/n],{i,1,k}} 可以求出1~trunc(sqrt(n))中的n的因子并求出对称的因子,通过n div (n div i)得到上限
X mod Y =X-[X/Y]
=>sum{k mod i,{i,1,k}}=sum{i,{i,1,k}}-sum{[i/n],{i,1,k}}
sum{i,{i,1,k}}=k*(k+1)/2
sum{[i/n],{i,1,k}} 可以求出1~trunc(sqrt(n))中的n的因子并求出对称的因子,通过n div (n div i)得到上限