#pragma GCC optimize("Ofast","inline")
#include<bits/stdc++.h>
#define int long long
using namespace std;
const int N=152,M=N*(N+1)/2,mod=1e9+7;

int n,m,f[N][M][2],t,inv[M],ans;

int ksm(int a,int b)
{
    int ans=1;
    while(b)
    {
        if(b&1)ans=ans*a%mod;
        b>>=1;a=a*a%mod;
    }
    return ans;
}

int g(int len)
{
    return len*(len+1)/2;
}

void ad(int&x,int y)
{
    x+=y;x-=x>=mod?mod:0;
}

void add(int s,int k)
{
    ad(ans,inv[k]*s%mod);
}

main()
{
    cin>>n>>m;t=g(n);m=min(m,n-m+1);
    for(int i=0;i<=t;i++)inv[i]=ksm(1+mod-i*ksm(t,mod-2)%mod,mod-2);
    for(int i=1;i<=m;i++)
    {
        memset(f,0,sizeof f);
        f[i][0][0]=mod-1;
        f[i][0][1]=(i-1)*(mod-1)%mod;
        for(int j=i;j<=n;j++)
        for(int k=0;k<=g(j-1);k++)
        for(int o=0;o<2;o++)if(f[j][k][o])
        {
            int s=f[j][k][o],u=j<m?m:n;
            if(j>=m)
            {
                if(o)add(s,k+g(i+n-j));
                else add(s*(n-j)%mod,k+g(i+n-j));
            }
            for(int p=j+1;p<=u;p++)
            {
                ad(f[p][k+g(p-j-1)][o],mod-s);
                if(!o)ad(f[p][k+g(p-j-1)][1],(mod-s)*(p-j-1)%mod);
            }
        }
    }
    cout<<ans*ksm(n-1,mod-2)%mod<<'\n';
}