这道题初看确实没什么思路,感觉之前的数论知识都用不上,只好自己找规律
首先当n>=k 这部分是很容易直接算出的
下面我们先来尝试这穷举i,
不难发现当穷举i时,总存在一段连续的除数,k div i=p定值
设这段是i~j,则这部分的的余数和signma(k-p*q) (i<=q<=j) 即为k*(j-i+1)-p*(i+j)*(j-i+1)/2
由于随着i的增大,k div i逐渐变小,是接近单调的
因此这样一段连续的除数我们可以通过二分确定其范围
这样就可以AC了

 var ans:int64;
n,k,i,j,p:longint; function find(l,r:longint):longint;
var m,w:longint;
begin
w:=l;
while l<=r do
begin
m:=(l+r) shr ;
if (k div m=p) then
begin
w:=m;
l:=m+;
end
else r:=m-;
end;
exit(w);
end; begin
readln(n,k);
if n>k then
begin
ans:=ans+int64(n-k)*int64(k);
n:=k-;
end;
if n=k then dec(n);
i:=;
while i<=n do
begin
p:=k div i;
j:=find(i,n);
ans:=ans-int64(j-i+)*int64(i+j)*int64(p) div ;
i:=j+;
end;
ans:=ans+int64(n)*int64(k);
writeln(ans);
end.
05-11 19:23