#include <bits/stdc++.h>
using namespace std;
long long mul(long long x, long long y) {
  return ((x % 1000000007) * (y % 1000000007)) % 1000000007;
}
long long add(long long x, long long y) {
  return ((x % 1000000007) + (y % 1000000007)) % 1000000007;
}
long long power(long long a, long long b) {
  if (b == 0) return 1;
  if (b == 1) return a;
  long long temp = power(a, b / 2);
  if (b % 2 == 0) {
    return mul(temp, temp);
  } else {
    return mul(a, mul(temp, temp));
  }
}
long long div1(long long x, long long y) {
  y = power(y, 1000000007 - 2);
  return ((x % 1000000007) * (y % 1000000007)) % 1000000007;
}
int main() {
  ios_base::sync_with_stdio(false);
  cin.tie(NULL);
  long long a, b, i, x, y;
  cin >> a >> b;
  long long ans = 0, two = 2;
  for (i = 1; i <= 1; i++) {
    x = div1(mul(i, mul(b, mul(a, add(a, 1)))), two);
    y = mul(i, a);
    ans = add(ans, add(x, y));
  }
  ans = div1(mul(b - 1, mul(b, ans)), two);
  cout << ans << endl;
  return 0;
}
