#include<bits/stdc++.h>
using namespace std;
const int mod=1e9+7;
int qpow(int a,int b){
    int c=1;for(;b;b>>=1){
        if(b&1)c=1ll*a*c%mod;
        a=1ll*a*a%mod;
    }
    return c;
}
int n,m,dp[160][12010][2],w1[160],w2[160],t;
void ad(int &x,int y){x+=y;if(x>=mod)x-=mod;}
int main(){
    scanf("%d%d",&n,&m);t=n*(n+1)/2;
    for(int i=1;i<=n;i++)w1[i]=i*(i-1)/2,w2[i]=i*(i+1)/2;
    dp[0][0][0]=1;
    for(int i=0;i<n;i++)for(int j=0;j<=t;j++)for(int o=0;o<2;o++)if(dp[i][j][o])
    for(int k=i+1;k<=n;k++){
        int wo=(i<m&&k>=m)?w2[k-i]:w1[k-i];
        if(!o)ad(dp[k][j+wo][1],mod-dp[i][j][o]);
        ad(dp[k][j+wo][o],mod-dp[i][j][o]);
    }
    int ans=0;
    for(int j=0;j<t;j++){
        int e=1ll*t*qpow(t-j,mod-2)%mod;
        ad(ans,1ll*e*dp[n][j][0]%mod*n%mod*qpow(n-1,mod-2)%mod);
        ad(ans,mod-1ll*e*dp[n][j][1]%mod*qpow(n-1,mod-2)%mod);
    }
    printf("%d",ans);
    return 0;
}