#include <cstdio>
#include <cstdlib>
#include <cstring>
#include <cmath>
#include <cassert>
#include <iostream>
#include <sstream>
#include <vector>
#include <queue>
#include <set>
#include <map>
#include <utility>
#include <numeric>
#include <algorithm>
#include <bitset>
#include <complex>
#include <array>
#include <list>
#include <stack>
#include <valarray>

using namespace std;

typedef unsigned uint;
typedef long long Int;
typedef unsigned long long UInt;

const int INF = 1001001001;
const Int INFLL = 1001001001001001001LL;

template<typename T> void pv(T a, T b) { for (T i = a; i != b; ++i) cout << *i << " "; cout << endl; }
template<typename T> void chmin(T& a, T b) { if (a > b) a = b; }
template<typename T> void chmax(T& a, T b) { if (a < b) a = b; }
int in() { int x; scanf("%d", &x); return x; }
double fin() { double x; scanf("%lf", &x); return x; }
Int lin() { Int x; scanf("%lld", &x); return x; }

int N, K;
Int memo[128][11111];

const Int MO = 1000000007;

Int powmod(Int x, int e) {
  Int v = 1;
  for (; e; e=e>>1, x=(x*x)%MO) {
    if (e&1) v=(v*x)%MO;
  }
  return v;
}

Int solve(int pos, int add) {
  if (pos < 0) return 0;
  if (memo[pos][add] >= 0) {
    return memo[pos][add];
  }
  Int& res = memo[pos][add];
  res = 0;

  for (int u = 0; u <= K; ++u) {
    if (u <= pos + 1) {
      int last = (u + add) % (pos + 1);
      int aa = (u + add) / (pos + 1);
      (res += solve(pos - 1, aa + add)) %= MO;
      (res += powmod(K + 1, pos) * last % MO) %= MO;
    } else {
      int last = (u + add);
      (res += solve(pos - 1, add)) %= MO;
      (res += powmod(K + 1, pos) * last % MO) %= MO;
    }
  }
  return res;
}

int main() {
  N = in();
  K = in();
  for (int i = 0; i < 128; ++i) {
    for (int j = 0; j < 11111; ++j) {
      memo[i][j] = -1;
    }
  }

  printf("%lld\n", solve(N - 1, 0));

  return 0;
}
