#include <bits/stdc++.h>
using namespace std;
#define int long long
const int N=155,M=3e4+5,mods=1e9+7;
int pows(int a,int b){
	if(b==0)return 1;
	int res=pows(a,b>>1);
	res=res*res%mods;
	if(b&1)res=res*a%mods;
	return res;
}
int n,m,inv[M],res,f[2][N][M][2],rs[N][M][2],gs[M];
int C(int x){
	return x*(x+1)/2;
}
signed main(){
	cin>>n>>m;
	int K=n*(n+1)/2;
	for(int i=1;i<=n*n;i++)inv[i]=pows(i,mods-2);
	f[0][0][0][0]=-1;
	for(int i=1;i<=n+1;i++){
		for(int j=0;j<i;j++){
			for(int k=0;k+C(i-j-1)<=K;k++){
				f[0][j][k][0]%=mods;
				f[0][j][k][1]%=mods;
				if(i==m){
					f[1][j][k][0]+=f[0][j][k][0];
					f[1][j][k][1]+=f[0][j][k][1];
				}else{
					f[1][j][k][0]+=f[0][j][k][0];
					f[1][i][k+C(i-j-1)][0]-=f[0][j][k][0];
					f[1][j][k][1]+=f[0][j][k][1];
					f[1][i][k+C(i-j-1)][1]-=f[0][j][k][1];
					f[1][j][k][1]+=f[0][j][k][0];
				}
			}
		}
		swap(f[1],f[0]);
		memset(f[1],0,sizeof(f[1]));
	}
	for(int i=0;i<=K;i++)res-=inv[K-i]*(f[0][n+1][i][1]%mods)%mods;
	cout<<(res%mods*K%mods*inv[n-1]%mods+mods)%mods;
}