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
#include <iostream>
#include <cstdint>
#include <vector>
#include <algorithm>
#include <iterator>


using namespace std;

constexpr uint64_t p = 1'000'000'007;

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

    int64_t n,m;
    cin>>n>>m;

    vector <vector<uint64_t>> otwartych( n+1, vector<uint64_t> (n+1,0) );
    vector <vector<uint64_t>> domknietych( n+1, vector<uint64_t> (n+1,0) );
    ///tab[ pozycja, liczba koncowek]


    otwartych[1][1]=m;

    for (int pozycja=2; pozycja <=n; pozycja ++ ){
        for (int koncowek=1; koncowek<=pozycja; koncowek++ ){
            otwartych[pozycja][koncowek] =  ((m-koncowek)*otwartych[pozycja-1][koncowek] + ( m - (koncowek-1))*domknietych[pozycja-1][koncowek-1])%p;
            domknietych[pozycja][koncowek] =  ( (koncowek)*otwartych[pozycja-1][koncowek] + ((koncowek))*domknietych[pozycja-1][koncowek]     )%p;
        }

    }
    uint64_t aku=0;
    for (int i=1; i<=n; i++){
        aku = (aku + domknietych[n][i])%p;
    }
    cout<<aku<<endl;

    return 0;
}