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 ;
       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.