#include<bits/stdc++.h>
#define int long long
using namespace std;
const int mn=2e6+5,mod=1e9+7;
int d[mn],n,m;
int fpow(int x,int y)
{
    int res=1;
    while(y)
    {
        if(y&1)
        {
            res=res*x%mod;
        }
        x=x*x%mod;
        y>>=1;
    }
    return res;
}
int jc[mn],inv[mn];
void init()
{
    jc[0]=1;
    for(int i=1;i<=2e6;i++)
    {
        jc[i]=jc[i-1]*i%mod;
    }
    inv[2000000]=fpow(jc[2000000],mod-2);
    // cerr<<inv[2000000]<<'\n';
    for(int i=2e6;i>=1;i--)
    {
        inv[i-1]=inv[i]*i%mod;
    }
    d[0]=1;
    for(int i=2;i<=2e6;i++)
    {
        d[i]=(i-1)*((d[i-1]+d[i-2])%mod)%mod;
    }
}
int a(int x,int y)
{
    if(x<y || y<0)return 0;
    return jc[x]*inv[x-y]%mod;
}
int c(int x,int y)
{
    if(x<y || y<0)return 0;
    // cerr<<jc[x]<<" "<<inv[x-y]<<'\n';
    return jc[x]*inv[y]%mod*inv[x-y]%mod;
}
void solve()
{
    scanf("%lld%lld",&n,&m);
    int ans=0;
    int f=1,bt=1;
    for(int i=0;i<=n;i++)
    {
        // bt*=i;
        ans+=c(n,i)*a(m,i)%mod*a(m-i,n-i)%mod*a(m-i,n-i)%mod*f;
        // cerr<<n<<" "<<i<<" "<<c(n,i)<<" "<<a(m,i)<<" "<<a(m-i,n-i)<<" "<<a(m-i,n-i)<<'\n';
        ans%=mod;
        ans+=mod;
        ans%=mod;
        f*=-1;
    }
    printf("%lld\n",ans);
}
signed main()
{
    init();
    int T;
    scanf("%lld",&T);
    while(T--)
    {
        solve();
    }
}