bzoj1257

这道题初看确实没什么思路，感觉之前的数论知识都用不上，只好自己找规律
首先当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 1;
      if (k div m=p) then
      begin
        w:=m;
        l:=m+1;
      end
      else r:=m-1;
    end;
    exit(w);
  end;

begin
  readln(n,k);
  if n>k then
  begin
    ans:=ans+int64(n-k)*int64(k);
    n:=k-1;
  end;
  if n=k then dec(n);
  i:=1;
  while i<=n do
  begin
    p:=k div i;
    j:=find(i,n);
    ans:=ans-int64(j-i+1)*int64(i+j)*int64(p) div 2;
    i:=j+1;
  end;
  ans:=ans+int64(n)*int64(k);
  writeln(ans);
end.