#ifdef _DEBUG
#include ".vscode/debug_print.hpp"
#else
#define debug(...)
#endif
#include <bits/stdc++.h>

using ll = long long;
using ld = long double;
using namespace std;
constexpr int INF32       = 1050000000;
constexpr long long INF64 = 4000000000000000000;
constexpr int MOD7        = 1000000007;
constexpr int MOD53       = 998244353;
constexpr long double PI  = 3.14159265358979323846;
template <class T> inline bool chmin(T &a, T b) {
    if (a > b) {
        a = b;
        return true;
    }
    return false;
}
template <class T> inline bool chmax(T &a, T b) {
    if (a < b) {
        a = b;
        return true;
    }
    return false;
}
void print() { std::cout << '\n'; }
template <class H, class... T> void print(H &&head, T &&... args) {
    std::cout << head;
    sizeof...(args) == 0 ? std::cout << "" : std::cout << ' ';
    print(std::forward<T>(args)...);
}
template <class T> void print(std::vector<T> &v) {
    for (int i = 0; i < v.size(); i++) {
        std::cout << v[i];
        i == v.size() - 1 ? std::cout << '\n' : std::cout << ' ';
    }
}
template <class T> void print(std::vector<std::vector<T>> &v) {
    for (int i = 0; i < v.size(); i++) {
        for (int j = 0; j < v[i].size(); j++) {
            std::cout << v[i][j];
            j == v[i].size() - 1 ? std::cout << '\n' : std::cout << ' ';
        }
    }
}
void scan() {}
template <class H, class... T> void scan(H &&head, T &&... args) {
    std::cin >> head;
    scan(std::forward<T>(args)...);
}
template <class T> void scan(std::vector<T> &v) {
    for (auto &&i : v) {
        scan(i);
    }
}

void vprint() { std::cout << '\n'; }
template <class H, class... T> void vprint(H &&head, T &&... args) {
    std::cout << head;
    sizeof...(args) == 0 ? std::cout << "" : std::cout << '\n';
    vprint(std::forward<T>(args)...);
}
template <class T> void vprint(std::vector<T> &v) {
    for (auto &&i : v) {
        std::cout << i << '\n';
    }
}

template <int mod> class modint {
  private:
    long long val;

  public:
    constexpr modint(long long x = 0) noexcept : val(x % mod) {
        if (val < 0)
            val += mod;
    };
    constexpr modint operator-() const noexcept { return val ? mod - val : 0; }
    constexpr modint operator+(const modint &x) const noexcept {
        return modint(*this) += x;
    }
    constexpr modint &operator++() noexcept {
        val += 1;
        if (mod <= val)
            val -= mod;
        return *this;
    }
    constexpr modint operator++(int) noexcept {
        modint<mod> ret = *this;
        val += 1;
        if (mod <= val)
            val -= mod;
        return ret;
    }
    constexpr modint operator-(const modint &x) const noexcept {
        return modint(*this) -= x;
    }
    constexpr modint &operator--() noexcept {
        val -= 1;
        if (val < 0)
            val += mod;
        return *this;
    }
    constexpr modint operator--(int) noexcept {
        modint<mod> ret = *this;
        val -= 1;
        if (val < 0)
            val += mod;
        return ret;
    }
    constexpr modint operator*(const modint &x) const noexcept {
        return modint(*this) *= x;
    }
    constexpr modint operator/(const modint &x) const noexcept {
        return modint(*this) /= x;
    }
    constexpr modint &operator+=(const modint &x) noexcept {
        val += x.val;
        if (mod <= val)
            val -= mod;
        return *this;
    }
    constexpr modint &operator-=(const modint &x) noexcept {
        val -= x.val;
        if (val < 0)
            val += mod;
        return *this;
    }
    constexpr modint &operator*=(const modint &x) noexcept {
        val = val * x.val % mod;
        return *this;
    }
    constexpr modint &operator/=(const modint &x) noexcept {
        val *= x.inv();
        return *this;
    }
    constexpr modint inv() const noexcept {
        long long a = val, b = mod, u = 1, v = 0, t;
        while (0 < b) {
            t = a / b;
            a -= t * b;
            swap(a, b);
            u -= t * v;
            swap(u, v);
        }
        return modint(u);
    }
    constexpr bool operator==(const modint &x) const noexcept {
        return this->val == x.val;
    }
    constexpr bool operator!=(const modint &x) const noexcept {
        return this->val != x.val;
    }
    constexpr bool operator<(const modint &x) const noexcept {
        return val < x.val;
    }
    constexpr bool operator<=(const modint &x) const noexcept {
        return val <= x.val;
    }
    constexpr bool operator>(const modint &x) const noexcept {
        return x.val < val;
    }
    constexpr bool operator>=(const modint &x) const noexcept {
        return x.val <= val;
    }
    friend constexpr ostream &operator<<(ostream &os, modint<mod> &x) noexcept {
        return os << x.val;
    }
    friend constexpr istream &operator>>(istream &is, modint<mod> &x) noexcept {
        return is >> x.val;
    }
    friend constexpr modint<mod> power(const modint<mod> &radix,
                                       long long x) noexcept {
        if (x == 0)
            return 1;
        auto t = power(radix, x / 2);
        t      = t * t;
        if (x & 1)
            t = t * radix;
        return t;
    }
};

using mint = modint<MOD7>;

int main() {
    std::cin.tie(0);
    std::ios::sync_with_stdio(false);

    int N, K;
    scan(N, K);

    mint ans = 0;
    auto sum = [](ll l, ll r) { return (l + r) * (r - l + 1) / 2; };
    for (size_t i = K; i <= N + 1; i++) {
        ll l = sum(0, i - 1);
        ll r = sum(N - i + 1, N);
        ans += r - l + 1;
    }
    print(ans);

    return 0;
}

Measured memory_bytes: 3660