#include <bits/stdc++.h>
using namespace std;

#define rep(i,x,y) for(ll i=(x);i<(y);i++)
#define rrep(i,x,y) for(ll i=(ll)(y)-1;i>=(x);i--)
#define all(x) (x).begin(),(x).end()
#ifdef LOCAL
#define debug(x) cerr << #x << " = " << (x) << " (L" << __LINE__ << ")" << " " << __FILE__ << endl
#define itrdebug(x) cerr << #x << " "; for (auto & el : (x)) {cerr << (el) << " ";} cerr << endl
#define dassert(...) assert(__VA_ARGS__)
#else
#define debug(x)
#define itrdebug(x)
#define dassert(...)
#endif
//#define int long long

typedef long long ll;
const ll MOD = 1e9 + 7;
const long double EPS = 1e-8;

namespace mod {
    long long pow(long long a, long long x, long long p) {
        long long res = 1;
        while(x > 0) {
            if( x % 2 != 0) {
                res = (res * a) % p;
            }
            a = (a * a) % p;
            x /= 2;
        }
        return res;
    }

    long long inverse(long long a, long long p) {
        //calculates the modular multiplicative of a mod m.
        //(assuming p is prime).
        return pow(a, p-2, p);
    }

    long long factorial(long long a, long long p) {
        long long fact = 1;
        for (long long i=2; i<=a; i++) {
            fact *= i;
            fact %= p;
        }
        return fact;
    }

    template <long long MOD>
    struct Int {
        long long x;
        Int(long long x = 0) : x((x % MOD + MOD) % MOD){}
        Int operator-() const { return Int<MOD>(-x);}
        Int& operator+=(const Int<MOD> a) {
            if ((x += a.x) >= MOD) x -= MOD;
            return *this;
        }
        Int& operator-=(const Int<MOD> a) {
            if ((x += MOD-a.x) >= MOD) x -= MOD;
            return *this;
        }
        Int& operator*=(const Int<MOD> a) {
            (x *= a.x) %= MOD; return *this;
        }
        Int operator+(const Int<MOD> a) const { return Int<MOD>(*this) += a;}
        Int operator-(const Int<MOD> a) const { return Int<MOD>(*this) -= a;}
        Int operator*(const Int<MOD> a) const { return Int<MOD>(*this) *= a;}
        Int pow(long long t) const {
            if (!t) return 1;
            Int<MOD> a = pow(t >> 1);
            a *= a;
            if (t&1) a *= *this;
            return a;
        }

        Int inv() const { return pow(MOD - 2);}
        Int& operator/=(const Int<MOD> a) { return *this *= a.inv();}
        Int operator/(const Int<MOD> a) const { return Int<MOD>(*this) /= a;}
    };
}
std::istream& operator>>(std::istream& is, const mod::Int<MOD>& a) { return is >> a.x;}
std::ostream& operator<<(std::ostream& os, const mod::Int<MOD>& a) { return os << a.x;}

// #include "mod.hpp"

namespace mod {
    struct Combination {
        Combination(long long max_n, long long p) : _mod(p) {
            factorials_ = std::vector<long long>(max_n+1, 1);
            inverse_factorials_ = std::vector<long long>(max_n+1, 1);
            for (long long i=1; i<max_n+1; i++) {
                factorials_[i] = factorials_[i-1] * i;
                factorials_[i] %= p;
            }
            inverse_factorials_[max_n] = mod::pow(factorials_[max_n],  p-2, p);
            for (long long i=max_n-1; i>=1; i--) {
                inverse_factorials_[i] = inverse_factorials_[i+1] * (i+1);
                inverse_factorials_[i] %= p;
            }
        }

        long long operator()(long long n, long long k) {
            long long result = factorials_[n];
            result *= inverse_factorials_[k];
            result %= _mod;
            result *= inverse_factorials_[n-k];
            result %= _mod;
            return result;
        }

    private:
        std::vector<long long> factorials_, inverse_factorials_;
        long long _mod;
    };
}

long long combination(long long n, long long k) {
    long long numerator = 1; // n * (n-1) * ... * (n-k+1)
    for (int i=0; i<k; i++) {
        numerator = numerator * (n-i);
    }

    long long denominator = 1; // k!
    for (int i=1; i<=k; i++) {
        denominator = denominator * i;
    }

    return numerator / denominator;
}

void solve(long long N, long long M){
    if (N==1 && M==1) {
        cout << 0 << endl;
        return;
    }
    mod::Combination comb(5e5, MOD);
    mod::Int<MOD> answer = 0;
    vector<mod::Int<MOD>> factorials(N+1);
    factorials[0] = 1;
    rep(i,1,N+1) factorials[i] = factorials[i-1] * i;
    rep(i,0,N+1) {
        mod::Int<MOD> cnt = factorials[N] * comb(M,N) * comb(N,i) * comb(M-i,N-i) * factorials[N-i];
        if (i%2==0) {
            answer += cnt;
        } else {
            answer -= cnt;
        }
    }
    cout << answer << endl;
}

signed main(){
    // ios_base::sync_with_stdio(false);
    // cin.tie(NULL);
    long long N;
    scanf("%lld",&N);
    long long M;
    scanf("%lld",&M);
    solve(N, M);
    return 0;
}


Measured memory_bytes: 14712