#include <iostream>
#include <cstdio>
#include <cstdlib>
#include <algorithm>
#include <cmath>
#include <vector>
#include <set>
#include <map>
#include <unordered_set>
#include <unordered_map>
#include <queue>
#include <ctime>
#include <cassert>
#include <complex>
#include <string>
#include <cstring>
using namespace std;

#ifdef LOCAL
	#define eprintf(...) fprintf(stderr, __VA_ARGS__)
#else
	#define eprintf(...) 42
#endif

typedef long long ll;
typedef pair<int, int> pii;
#define mp make_pair

const int MOD = (int)1e9 + 7;
int add(int x, int y) {
	x += y;
	if (x >= MOD) return x - MOD;
	return x;
}
int sub(int x, int y) {
	x -= y;
	if (x < 0) return x + MOD;
	return x;
}
int mult(int x, int y) {
	return ((ll)x * y) % MOD;
}

const int N = 102;
const int S = N * N;
int dp[N][S];
int ANS;
int n, k;
int prod[N];

int main()
{
//	freopen("input.txt", "r", stdin);
//	freopen("output.txt", "w", stdout);

	scanf("%d%d", &n, &k);

	prod[0] = 1;
	for (int i = 1; i <= n; i++)
		prod[i] = mult(prod[i - 1], k + 1);

	dp[n][0] = 1;
	for (int p = n; p > 0; p--)
		for (int x = 0; x <= k * (n - p); x++)
			for (int y = 0; y <= k; y++) {
				if (y > p) {
					dp[p - 1][x] = add(dp[p - 1][x], dp[p][x]);
					continue;
				}
				int z = x + (x + y) / p;
				dp[p - 1][z] = add(dp[p - 1][z], dp[p][x]);
			}
	ANS = k * (k + 1) / 2;
	ANS = mult(ANS, prod[n - 1]);
	ANS = mult(ANS, n);

	for (int p = 1; p <= n; p++)
		for (int x = 0; x <= k * (n - p); x++)
			for (int y = 0; y <= k; y++) {
				if (y > p) continue;
				int z = (x + y) / p;
				ANS = sub(ANS, mult(z, mult(dp[p][x], prod[p - 1])));
			}
	printf("%d\n", ANS);

	return 0;
}
