#include <iostream> #include <cmath> using namespace std; const int MOD = 1000000009; /* int mnoz_iter (int val, int czynnik, int cnt) { if (cnt == 0) return val; return mnoz_iter(val * czynnik, czynnik + 1, cnt - 1); } int newton (int n, int k) // n po k { return mnoz_iter(1, n - k + 1, k)/mnoz_iter(1, 1, k); }*/ int pom(int n, int m) { int sum = 0; for(int i = 2; i < n - 1; i++) { sum = (sum + (int(m * pow(m, i - 2)) % MOD * (pom(n - i, m) % MOD)) % MOD) % MOD; } sum = (sum + int(m * pow(m, n - 2)) % MOD) % MOD; return sum; } int drwal(int n, int m) //zwraca odpowiedz do problemu, n miejsc, m rodzajow drzew { int sum = 0; for(int i = 2; i < n - 1; i++) { sum = (sum + (int(m * pow(m, i - 2)) % MOD * (drwal(n - i, m) % MOD)) % MOD) % MOD; } sum = (sum + int(m * pow(m, n - 2)) % MOD) % MOD; int temp = pom(n, m); temp -= m * pow(m, n - 2); temp /= m; sum = (sum - temp) % MOD; return sum; } int main() { int n, m; cin >> n >> m; if(n == 1) { cout << 0; return 0; } cout << drwal(n, m); return 0; }
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 51 52 53 54 55 56 57 58 59 60 61 62 63 64 | #include <iostream> #include <cmath> using namespace std; const int MOD = 1000000009; /* int mnoz_iter (int val, int czynnik, int cnt) { if (cnt == 0) return val; return mnoz_iter(val * czynnik, czynnik + 1, cnt - 1); } int newton (int n, int k) // n po k { return mnoz_iter(1, n - k + 1, k)/mnoz_iter(1, 1, k); }*/ int pom(int n, int m) { int sum = 0; for(int i = 2; i < n - 1; i++) { sum = (sum + (int(m * pow(m, i - 2)) % MOD * (pom(n - i, m) % MOD)) % MOD) % MOD; } sum = (sum + int(m * pow(m, n - 2)) % MOD) % MOD; return sum; } int drwal(int n, int m) //zwraca odpowiedz do problemu, n miejsc, m rodzajow drzew { int sum = 0; for(int i = 2; i < n - 1; i++) { sum = (sum + (int(m * pow(m, i - 2)) % MOD * (drwal(n - i, m) % MOD)) % MOD) % MOD; } sum = (sum + int(m * pow(m, n - 2)) % MOD) % MOD; int temp = pom(n, m); temp -= m * pow(m, n - 2); temp /= m; sum = (sum - temp) % MOD; return sum; } int main() { int n, m; cin >> n >> m; if(n == 1) { cout << 0; return 0; } cout << drwal(n, m); return 0; } |