#include <iostream>
#include <fstream>
#include <cassert>
#include <typeinfo>
#include <vector>
#include <stack>
#include <cmath>
#include <set>
#include <map>
#include <string>
#include <algorithm>
#include <cstdio>
#include <queue>
#include <iomanip>
#include <cctype>
#include <random>
#include <complex>
#define syosu(x) fixed<<setprecision(x)
using namespace std;
typedef long long ll;
typedef unsigned long long ull;
typedef pair<int,int> P;
typedef pair<double,double> pdd;
typedef pair<ll,ll> pll;
typedef vector<int> vi;
typedef vector<vi> vvi;
typedef vector<double> vd;
typedef vector<vd> vvd;
typedef vector<ll> vl;
typedef vector<vl> vvl;
typedef vector<char> vc;
typedef vector<vc> vvc;
typedef vector<string> vs;
typedef vector<bool> vb;
typedef vector<vb> vvb;
typedef vector<P> vp;
typedef vector<vp> vvp;
typedef vector<pll> vpll;
typedef pair<int,P> pip;
typedef vector<pip> vip;
const int inf=1<<30;
const ll INF=1ll<<59;
const double pi=acos(-1);
const double eps=1e-7;
const ll mod=1e9+7;
const int dx[4]={0,1,0,-1},dy[4]={1,0,-1,0};

const int M=4000;

int n,k;

int main(){
	cin>>n>>k;
	vl dp,DP,dp1,DP1;
	dp=DP=dp1=DP1=vl(M);
	dp1[0]=1;
	int L=0;
	for(int i=n;i>0;i--){
		for(int j=0;j<M;j++) DP[j]=DP1[j]=0;
		for(int j=0;j<=k;j++) for(int l=0;l<=L;l++){
			if(j<=i){
				int I=l+(l+j)/i;
				(DP1[I]+=dp1[l])%=mod;
				(DP[I]+=dp[l]+(l+j)%i*dp1[l])%=mod;
			}
			else{
				(DP1[l]+=dp1[l])%=mod;
				(DP[l]+=dp[l]+(l+j)*dp1[l])%=mod;
			}
		}
		dp=DP;
		dp1=DP1;
		L+=(i+L)/i;
	}
	int res=0;
	for(int i=0;i<M;i++) (res+=dp[i])%=mod;
	cout<<res<<endl;
}