1
 2
 3
 4
 5
 6
 7
 8
 9
10
11
12
13
14
15
16
17
18
19
20
21
22
23
24
25
26
27
28
29
30
31
32
33
34
35
36
37
38
39
40
41
42
43
44
45
46
47
48
49
50
#include <bits/stdc++.h>

using namespace std;
long long MOD = 1000000007ll;

long long DP[3001][1501][2];

long long solve() {
    long long n,m;
    cin>>n>>m;

    DP[1][1][0]=m;

    long long x = 1;
    for(long long i=2;i<=n;i++) {
        for(long long j=1;j<=x;j++) {
            // koniec = 1
            // konczymy cykl
            DP[i][j][1] += (DP[i-1][j][1]*j)%MOD;
            DP[i][j][1]%=MOD;
            // rozpoczynamy cykl
            if(m-j>0 && DP[i-1][j][1]>0) {
                DP[i][j+1][0] += (DP[i-1][j][1]*(m-j))%MOD;
                DP[i][j+1][0]%=MOD;
                x = max(x,j+1);
            }
            //koniec = 0
            DP[i][j][1] += (DP[i-1][j][0]*j)%MOD;
            DP[i][j][1]%=MOD;
            // rozpoczynamy cykl
            if(m-j>0) {
                DP[i][j][0] += (DP[i-1][j][0]*(m-j))%MOD;
                DP[i][j][0]%=MOD;
            }
        }
    }
    // jesli n to kończymy
    long long wynik = 0;
    for(long long j=1;j<=x;j++) {
        wynik+=DP[n][j][1];
        wynik%=MOD;
    }
    return wynik;
}

int main() {
    ios_base::sync_with_stdio(false);
    cout<<solve()<<"\n";
    return 0;
}