#include<bits/stdc++.h>
#define ll long long
#define fi first
#define se second
#define pb push_back
using namespace std;
const int N=155;
const int mod=1e9+7;
int n,m;
ll f[N][N*N][2],ifac[N*N],res;
int get(int x){
    return x*(x+1)/2;
}
void add(ll &x,ll y){
    x=(x+y)%mod;
}
ll fpow(ll x,ll y=mod-2){
    ll z(1);
    for(;y;y>>=1){
        if(y&1)z=z*x%mod;
        x=x*x%mod;
    }return z;
}
void init(int n){
    ifac[1]=1;for(int i=2;i<=n;i++)ifac[i]=mod-ifac[mod%i]*(mod/i)%mod;
}
int main(){
    //freopen("data.in","r",stdin);
    ios::sync_with_stdio(false);
    cin.tie(0),cout.tie(0);
    cin>>n>>m,init(n*(n+1));
    if(n==2){
        cout<<0;
        return 0;
    }
    f[0][0][0]=1;
    for(int i=0;i<=n-1;i++){
        for(int j=0;j<=i*(i+1)/2;j++){
            for(int k=0;k<2;k++){
                if(!f[i][j][k])continue;
                for(int l=i+1;l<=n;l++){
                    int op=(i<=n-m&&l>n-m);
                    add(f[l][j+get(l-i-1+op)][k],-f[i][j][k]);
                    if(!k)add(f[l][j+get(l-i-1+op)][1],-f[i][j][k]*(l-i-1));
                }
            }
        }
    }
    for(int i=0;i<n*(n+1)/2;i++){
        res=(res+f[n][i][1]*ifac[n*(n+1)/2-i])%mod;
    }
    res=res*fpow(n-1)%mod*n*(n+1)%mod*(mod+1>>1)%mod;
    cout<<(res+mod)%mod;
}