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
//Jakub Rozek
//Od deski do deski
//czas: n^2
//pamiec: n

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

const long long N=3003;
const long long mod=1000000007;

long long n,m,odp,p;
long long dp[2][N][2];

int main()
{
    ios_base::sync_with_stdio(0);
    cin.tie(0);

    cin>>n>>m;

    dp[1][1][0]=m;
    p=1;
    for(int i=2; i<=n; ++i)
    {
        p=1-p;
        for(int j=1; j<=i; ++j)
        {
            dp[p][j][0]=dp[1-p][j][0]*(m-j)+dp[1-p][j-1][1]*(m-j+1);
            dp[p][j][0]%=mod;
            dp[p][j][1]=dp[1-p][j][0]*(j)+dp[1-p][j][1]*(j);
            dp[p][j][1]%=mod;
        }
    }

    for(int j=1; j<=n; ++j)
    {
        odp+=dp[p][j][1];
        odp%=mod;
    }

    cout<<odp<<"\n";
    return 0;
}