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-1)*(x-2) div 6);
  end;

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

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