#include<bits/stdc++.h>
const int N = 105;
const int mod = 1e9 + 7;
using namespace std;

int n, k, f[N][N*N], ans;

void add(int&a, int b){
    a += b; if (a >= mod) a -= mod; if (a < 0) a += mod;
}

int main(){
    ios_base::sync_with_stdio(false); cin.tie(0); cout.tie(0);
    cin >> n >> k;
    f[n+1][0] = 1;
    int cnt = 0;
    ans = 1LL * k * (k+1) / 2 % mod;
    for (int i = 2; i <= n; i++) ans = 1LL * ans * (k+1) % mod;
    ans = 1LL * ans * n % mod;
    for (int i = n; i >= 1; i--) {
        int cur = cnt;
        cnt += (cnt + min(i, k)) / i;
        for (int j = 0; j <= cnt; j++){
            for (int l = min(cur, j); l >= 0; l--){
                int Min = max((j - l) * i - l, 0), Max = min(min(k, i), (j-l+1)*i-1-l);
                if (l == j){
                    if (k > i) add(f[i][j], 1LL * f[i+1][j] * (k-i) % mod);
                }
                if (Max < Min) continue;
                add(f[i][j], 1LL * f[i+1][l] * (Max - Min + 1) % mod);
            }
        }
    }
    for (int i = 0; i <= cnt; i++){
        add(ans, -1LL * f[1][i] * i % mod);
    }
    cout << ans << endl;
}
