#include<cstdio>
#include<cstring>
#include<algorithm>
#include<cmath>
#include<iostream>
using namespace std;
int getint()
{
    int res;
    scanf("%d",&res);
    return res;
}
int n,m;
const int N = 65;
const int M = 6000;
const int Mod = 1000000007;
int dp[N][N][N];
int zuhe[N][N];
int ec[N];
int erci[M];
void ad(int &x,int y)
{
    x+=y;
    x%=Mod;
}
int ja(int x,int y)
{
    x+=y;
    x%=Mod;
    return x;
}
int cg(int x,int y)
{
    long long z=x;
    z*=y;
    z%=Mod;
    x=z;
    return x;
}
void Init()
{
    int i,j,k,s,t;
    for(i=0;i<N;i++)
    {
        zuhe[i][0]=1;
        for(j=1;j<=i;j++)
            zuhe[i][j]=ja(zuhe[i-1][j-1],zuhe[i-1][j]);
    }
    ec[0]=1;
    for(i=1;i<N;i++)
        ec[i]=cg(ec[i-1],2);
    for(i=0;i<N;i++)
        ad(ec[i],Mod-1);
    erci[0]=1;
    for(i=1;i<M;i++)
        erci[i]=cg(erci[i-1],2);
    dp[1][0][1]=1;
    for(i=1;i<N;i++)
    {
        for(j=0;j<i;j++)
        {
            for(k=1;k<=i;k++)
            {
                int now=dp[i][j][k];
                for(s=1;s+i<N;s++)
                {
                    now=cg(now,ec[k]);
                    int val=now;
                    val=cg(val,erci[(s*(s-1))/2]);
                    val=cg(val,zuhe[s+i-1][s]);
                    ad(dp[s+i][j+1][s],val);
                }
            }
        }
    }
}
void f()
{
    int i,j,k;
    n=getint();
    m=getint();
    int ans=0;
    for(i=1;i<=n;i++)
    {
        for(j=0;j<m;j++)
        {
            for(k=1;k<=i;k++)
            {
                int now=dp[i][j][k];
                now=cg(now,zuhe[n-1][n-i]);
                now=cg(now,erci[((n-i)*(n-i-1))/2]);
                ad(ans,now);
            }
        }
    }
    printf("%d\n",ans);
}
int main()
{
    int T=getint();
    Init();
    while(T--)
        f();
    return 0;
}