#include <bits/stdc++.h>
using namespace std;
const long long mod = 1000000007ll;
void ad(long long &a, long long b) {
  a += b;
  a %= mod;
}
void mn(long long &a, long long b) {
  a += mod - b;
  a %= mod;
}
void ml(long long &a, long long b) {
  a *= b;
  a %= mod;
}
int main() {
  int N, K;
  scanf("%d%d", &N, &K);
  static long long P[200];
  P[0] = 1ll;
  for (int i = 0; i <= N; i++) {
    P[i + 1] = P[i] * (K + 1) % mod;
  }
  static long long dp[10010][200] = {0};
  long long ans = 0ll;
  dp[0][N] = 1;
  for (int i = N - 1; i >= 0; i--) {
    for (int y = 0; y <= N * N; y++) {
      if (!dp[y][i + 1]) continue;
      for (int x = 0; x <= K; x++) {
        if (x <= i + 1 && i + 1 - x <= y) {
          int p = (x + y) / (i + 1);
          ad(ans, P[i] * p % mod * dp[y][i + 1] % mod);
          ad(dp[y + p][i], dp[y][i + 1]);
        } else {
          ad(dp[y][i], dp[y][i + 1]);
        }
      }
    }
  }
  printf("%lld\n",
         (K * (K + 1) / 2 * P[N - 1] % mod * N % mod + mod - ans) % mod);
  return 0;
}
