#include <bits/stdc++.h>
using namespace std;
typedef long long ll;
typedef long double ld;
typedef vector<int> VI;
typedef vector<ll> VL;
typedef vector<vector<int>> VII;
typedef vector<vector<ll>> VLL;
typedef vector<string> VS;

#define pii pair<int,int>
#define pll pair<ll,ll>
#define pis pair<int,string>
#define psi pair<string,int>
#define rep(i,num,n) for(int i=num;i<(int)(n);i++) //for_loop
#define REP(i,n) for(int i=0;i<(int)(n);i++)
#define rrep(i,num,n) for(int i=num-1;i>=(int)(n);i--) //reverse_for>
#define in(x,a,b) (a<=x && x<b)//範囲
#define reo return 0
#define all(v) v.begin(),v.end()
#define INF 2000000000
#define INFL 9000000000000000000
#define sz(c) ((int)(c).size())
#define ten(x) ((int)1e##x)
#define tenll(x) ((ll)1e##x)

const long double pi = 3.1415926535897932384626433832795028841971L;
const ll MOD=1000000007;
const int num_fac=200001;
ll modfact(ll x){//階乗
    static ll _fact[num_fac+1];
    if(_fact[0]==0){
        _fact[0]=1;
        for(int i=1;i<=num_fac;++i)_fact[i]=_fact[i-1]*i%MOD;
    }
    return _fact[x];
}
ll modpow(ll a,ll n){//累乗
    ll r=1;
    while(n)r=r*((n%2)?a:1)%MOD,a=a*a%MOD,n>>=1;
    return r;
}
ll moddiv(ll a,ll b){//除算
    ll ap_2=modpow(b,MOD-2);
    return (a*ap_2)%MOD;
}

ll aCb(ll a,ll b){
    return moddiv(modfact(a),modfact(a-b)*modfact((b))%MOD);
}

int main(){
    ll n,m;cin>>n>>m;
    if(n<m){
        cout<<0;reo;
    }
    ll ans=0;
    int m2=1;
    for(int i=m;i>=1;i--){
        ll t=(aCb(m,i)*modpow(i,n))%MOD;
        ans+=m2*t;
        ans=(ans+MOD)%MOD;
        m2*=-1;
    }
    cout<<ans<<endl;
}