#include <bits/stdc++.h>
long long pi = acos(-1);
long long z = 1000000007;
long long inf = 100000000000000000;
long long p1 = 37;
long long p2 = 53;
long long mod1 = 202976689;
long long mod2 = 203034253;
long long gcd(long long a, long long b) {
  if (b > a) return gcd(b, a);
  if (b == 0) return a;
  return gcd(b, a % b);
}
long long power(long long a, long long b, long long p) {
  if (b == 0) return 1;
  long long c = power(a, b / 2, p);
  if (b % 2 == 0)
    return ((c * c) % p);
  else
    return ((((c * c) % p) * a) % p);
}
long long inverse(long long a, long long n) { return power(a, n - 2, n); }
long long max(long long a, long long b) {
  if (a > b) return a;
  return b;
}
long long left(long long i) { return ((2 * i) + 1); }
long long right(long long i) { return (2 * (i + 1)); }
using namespace std;
void solve() {
  long long n, m;
  cin >> n >> m;
  long long ans = (n % z) * (m % z);
  ans %= z;
  long long to_sub = 0;
  m = min(n, m);
  long long lim = m + 1;
  long long inv = inverse(2, z);
  for (long long k = 1; k * k <= n; k++) {
    long long l = (n) / (k + 1);
    l++;
    long long u = (n / k);
    u = min(m, u);
    if (l > m) continue;
    lim = min(lim, l);
    long long curr = (u % z) * ((u + 1) % z) - ((l % z) * ((l - 1) % z));
    curr %= z;
    curr *= inv;
    curr %= z;
    curr *= k;
    curr %= z;
    curr += z;
    curr %= z;
    to_sub += curr;
    to_sub %= z;
  }
  for (long long i = 1; i < lim; i++) {
    long long curr = (n / i);
    curr %= z;
    curr *= i;
    curr %= z;
    to_sub += curr;
    to_sub %= z;
  }
  ans -= to_sub;
  ans %= z;
  ans += z;
  ans %= z;
  cout << ans << endl;
  return;
}
int main() {
  ios_base::sync_with_stdio(0);
  cin.tie(0);
  cout.tie(0);
  long long t;
  t = 1;
  while (t--) {
    solve();
  }
  return 0;
}
