#include <bits/stdc++.h>
using namespace std;
const long long mod = 1000000007;
int gcd(int a, int b) { return b == 0 ? a : gcd(b, a % b); }
int func1(int a, int b) {
  return ((a * b) << 2) - ((a + b) << 1) - (gcd(a, b) << 1);
}
inline int func2(int a, int b) { return ((a * b) << 1) - (gcd(a, b) << 1); }
inline int func3(int a, int b) {
  return 5 * a * b - ((a + b) << 1) - (gcd(a, b) << 1);
}
int call(int a, int b) {
  if ((a & 1) && !(b & 1)) return func1(a, b);
  if (!(a & 1) && (b & 1)) return func1(b, a);
  if ((a & 1) && (b & 1))
    return func2(a, b);
  else
    return func3(a, b);
}
int main() {
  long long a, b, ans = 0;
  scanf("%lld%lld", &a, &b);
  for (int i = 1; i <= a; ++i)
    for (int j = 1; j <= b; ++j) {
      long long temp, sz = (a - i + 1) * (b - j + 1);
      if ((i & 1) && (j & 1))
        temp = 0;
      else
        temp = ((i - 1) << 1) + ((j - 1) << 1) + 4;
      temp = (temp + call(i, j)) % mod;
      temp = temp * sz % mod;
      ans = (ans + temp) % mod;
    }
  ans = ans * 6 % mod;
  printf("%lld\n", ans);
}
