#include <ext/pb_ds/assoc_container.hpp> 
#include <ext/pb_ds/tree_policy.hpp>  
#include <functional>
#include <iostream> 
#include<bits/stdc++.h>
using namespace __gnu_pbds;  
using namespace std;  
// Declaring ordered_set 
typedef tree<long long int, null_type, less<long long int>, rb_tree_tag,  
            tree_order_statistics_node_update>  
    ordered_set;  
// https://www.geeksforgeeks.org/order_of_key-in-c/
 
typedef pair<long long int,long long int> pa;
 
#define resize(A) A.resize(unique(A.begin(), A.end()) - A.begin())
#define pb push_back
#define  MAX 500005

#define endl "\n"                   //Take care in case of interactive problems


long long int mod=998244353;
long long int MOD=1e9+7;
long long int INF=1e18;
 
 
 
 
////////////////////////////////////////////// sieve  //////////////////////////////////////////////////////////////////

/*bool C[MAX];
vector<long long  int>V;
 void sieve()
 {
    long long int i,j;
    for(i=2;i<MAX;i++)
    {
      if(!C[i])
     { 
         
          V.pb(i);
        for(j=2*i;j<MAX;j=j+i)
        {
            C[j]=true;
        }
      }
        
    }
     
 }
 
 */
 
 
//////////////////////////////////////////////////////////D-S-U/////////////////////////////////////////////////////////
/*long long int parent[MAX];
long long int size[MAX];
void initialise()
{ //cout<<MAX<<" ";
    for(long  long int i=0;i<MAX;i++)
    {
        parent[i]=i;
        //cout<<i<<" ";
        size[i]=1;
    }
}
long long int  root(long long int x)
{
    while(x!=parent[x])
    {
        parent[x]=parent[parent[x]];
        x=parent[x];
    }
    return x;
}
 
    
void union1(long long int x,long long int y)
{
    if(size[x]>=size[y])
    {
        size[x]+=size[y];
        parent[y]=x;
        
    }
    else
    {
        size[y]+=size[x];
        parent[x]=y;
    }
}*/
//////////////////////////////////////////////////////MODPOWER//////////////////////////////////////////////////////////
 
 
long long int modulo(long long int a,long long  int b, long long int n){
    long long int  x=1, y=a%n; 
    while (b > 0) {
        if (b%2 == 1) {
            x = (x*y) % n; // multiplying with base
        }
        y = (y*y) % n; // squaring the base
        b /= 2;
    }
    return x % n;
}
 
 
////////////////////////////////////////////////BIT/////////////////////////////////////////////////////////////////////
/* long long int bit[MAX+1]={0};
long long int query(long long int in,long long int n)
{
    //in=in+1;
    long long int sum=0;
    if(in==0)
        return 0;
    while(in>0)
    {
      sum=sum+bit[in];
        
        in=in-(in&(-in));
    }
    return sum;
}
void update(long long int in,long long int n,long long int val)
{
    // in=in+1;
    
    while(in<=n)
    {
       // cout<<in<<" ";
        bit[in]+=val;
        in=in+ (in &(-in));
    }
}
 
 */
 

//////////////////////////////////////////////// DFS //////////////////////////////////////////////////////////////////
/*
  bool visited[MAX];
 long long int parent[MAX];
 long long int dist[MAX];


void dfs(vector<long long int>adj[],long long int i,long long int k)
{
    visited[i]=true;
    for(auto j:adj[i])
    {
        if(!visited[j])
        {
            dist[j]=dist[i]+1;
            parent[j]=i;
            dfs(adj,j,k);
            
        }
        
    }
}
 */
 

////////////////////////////////////////////////////////////////////////////////////////////////////////////////////////

 
long long int R[4]={1,0,-1,0};
long long int C[4]={0,1,0,-1};
 
 


/*
           Come Rain, Come Snow.
           Wherever your dreams take you, you gotta go.
*/



int main()
{
 ios_base::sync_with_stdio(false); 
 cin.tie(NULL);
 cout.tie(NULL);
long long int  sum,d,i,w,j,k=0,l=0,r=0,m,n,p=0,cnt=0,cnt1=0,x,y=0,z=0,a,b,c,low,mid,high,ans=0,ans1=0,maxi=0,mini=INF,t1,x1,y1,z1,t,o,i1,j1;
string s,s1,s2,s3,s4;
    
   long long int fact[MAX];
    fact[0]=1;
    for(i=1;i<MAX;i++)
    {
        fact[i]=fact[i-1]*i;
        fact[i]%=MOD;
    }
    cin>>n>>m;
    p=fact[m];
    r=modulo(fact[m-n],MOD-2,MOD);
    p=(p*r)%MOD;
    for(j=0;j<=n;j++)
    {
        if(j%2==0)
        {
            a=fact[n]*fact[m-j];
            a=a%MOD;
            b=fact[n-j]*fact[j];
            b=b%MOD;
            b=b*fact[m-n];
            b=b%MOD;
            b=modulo(b,MOD-2,MOD);
            a=a*b;
            a=a%MOD;
            cnt+=a;
            cnt%=MOD;
        }
        else
        {
            a=fact[n]*fact[m-j];
            a=a%MOD;
            b=fact[n-j]*fact[j];
            b=b%MOD;
            b=b*fact[m-n];
            b=b%MOD;
            b=modulo(b,MOD-2,MOD);
            a=a*b;
            a=a%MOD;
            cnt1+=a;
            cnt1%=MOD;
        }
    }
    cnt=(cnt-cnt1+MOD)%MOD;
    p=p*cnt;
    cout<<p%MOD;
    
}


Measured memory_bytes: 7564