#include <bits/stdc++.h>
using namespace std;

const int64_t MOD = 1e9+7;;
void add(int64_t& a, int64_t b){
    a = (a+b) % MOD;
}
void mul(int64_t& a, int64_t b){
    a = a*b % MOD;
}

int64_t power_mod(int64_t num, int64_t power){
    int64_t prod = 1;
    num %= MOD;
    while(power > 0){
        if(power&1) prod = prod * num % MOD;
        num = num * num % MOD;
        power >>= 1;
    }
    return prod;
}

int main(){
    int N, K;
    cin >> N >> K;
    int64_t ans = N*K*(K+1)/2 * power_mod(K+1, N-1) % MOD;

    static int64_t dp[101][101*101];
    dp[N][0] = 1;
    int mx = 0;
    for(int i=N; i>0; i--){
        int mx2 = mx;
        for(int j=0; j<=mx; j++) for(int k=0; k<=K; k++){
            int op = 0;
            if(k <= i) op = (j+k) / i;
            int64_t sub = op * dp[i][j] % MOD * power_mod(K+1, i-1) % MOD;
            add(ans, MOD - sub);
            add(dp[i-1][j+op], dp[i][j]);
            mx2 = max(mx2, j+op);
        }
        mx = mx2;
    }
    cout << ans << endl;
    return 0;
}
