Python :

# coding=utf-8
#!/python2
n,m=sorted(map(int,raw_input().split()))
mu=[i*i for i in xrange(0,n+1)]
for i in xrange(1,n+1):
	for j in xrange(i+i,n+1,i): mu[j]-=mu[i]
print (n*(n-1)*m*(m-1)*(11*m*n*(n+1)*(m+1)+6*(n*(n+1)+m*(m+1)))/144-sum(mu[d]*(n-d+n%d)*(n/d)/2*(m-d+m%d)*(m/d)/2 for d in xrange(1,n+1)))/3%1004535809