bzoj3505

ans=C((n+1)*(m+1),3)-三点一线的情况
横线竖线我们可以先去掉
然后考虑斜线，由于对称性我们只要考虑斜率大于0的即可
有一个很显然的结论，但两点坐标差为x,y时，这条线段上的点数为gcd(x,y)
我们设左下角点为(0,0),则两端点坐标差为x,y的线段有(n-x+1)*(m-y+1)
要注意同在一条直线上不能重复计算，我们考虑线段更容易一点
所以，对于每条线段，三点一线的情况为除去两端点的线段上点数

 var i,j,n,m:longint;
     ans,tmp:int64;
 function calc(x:int64):int64;
   begin
     exit(x*(x-)*(x-) div );
   end;

 function gcd(a,b:longint):longint;
   begin
     if b= then exit(a)
     else exit(gcd(b,a mod b));
   end;

 begin
   readln(n,m);
   ans:=calc((n+)*(m+))-(n+)*calc(m+)-(m+)*calc(n+);
   for i:= to n do
     for j:= to m do
     begin
       tmp:=gcd(i,j)+;
       if tmp> then
         ans:=ans-*(n-i+)*(m-j+)*(tmp-);
     end;
   writeln(ans);
 end.