#include <bits/stdc++.h>
using namespace std;
const long long N = 102;
const long long mod = 1e9 + 7;
void add(long long &a, long long b) {
  a += b;
  if (a >= mod) a -= mod;
  if (a < 0) a += mod;
}
long long qpow(long long a, long long exp) {
  long long ans = 1;
  while (exp) {
    if (exp & 1) ans = ans * a % mod;
    a = a * a % mod;
    exp >>= 1;
  }
  return ans;
}
long long f[N][N * N];
signed main() {
  ios_base::sync_with_stdio(0);
  cin.tie(0);
  cout.tie(0);
  long long n, k;
  cin >> n >> k;
  f[n + 1][0] = 1;
  for (long long i = n; i >= 1; --i)
    for (long long p = 0; p < n * k; ++p) {
      long long d1 = p + min(i, k);
      long long d2 = p - (p % i) + i;
      add(f[i][p + p / i], f[i + 1][p] * (min(d2 - 1, d1) + 1 - p) % mod);
      add(f[i][p + p / i + 1], f[i + 1][p] * max(d1 - d2 + 1, 0ll) % mod);
      add(f[i][p], f[i + 1][p] * max(k - i, 0ll) % mod);
    }
  long long ans = qpow(k + 1, n - 1) * k * (k + 1) * n / 2;
  ans %= mod;
  for (long long p = 0; p < n * k; ++p) add(ans, -p * f[1][p] % mod);
  cout << ans << endl;
}
