#include<bits/stdc++.h>

using namespace std;

const int mod = 1e9 + 7, N = 105;
int n, K, suma, dp[N][N * N], sumt;

int binpw (int aa, int mm) {
    int ret = 1;
    for (; mm; mm >>= 1, aa = (1LL * aa * aa) % mod) if(mm & 1) ret = (1LL * ret * aa) % mod;
    return ret;
}

int main () {
    ios_base::sync_with_stdio(0); cin.tie(0); cout.tie(0);
    cin >> n >> K;
    suma = (1LL * ( (K * (K + 1) ) / 2 * n) % mod * binpw(K + 1, n - 1) ) % mod;

    dp[n + 1][0] = 1;
    for (int i = n; i > 0; --i) {
        for (int j = 0; j < N * N; ++j) {
            if(!dp[i + 1][j]) continue;

            for (int k = 0; k <= K; ++k) {
                if(k <= i && j + (j + k) / i < N * N) {
                    dp[i][j + (j + k) / i] += dp[i + 1][j];
                    dp[i][j + (j + k) / i] %= mod;
                }
                if(k > i) {
                    dp[i][j] = (dp[i][j] + dp[i + 1][j]) % mod;
                }
            }
        }
    }

    for (int i = 0; i < N * N; ++i) sumt = (sumt + 1LL * i * dp[1][i]) % mod;

    cout << (suma - sumt + mod) % mod;

    return 0;
}
