#pragma GCC optimize("Ofast", "inline", "-ffast-math")
#include <bits/stdc++.h>
using namespace std;

using LL = long long;
constexpr int mod = 1E9 + 9;

struct mint {
    LL x;
    mint(int _x = 0) {
        x = _x;
    }

    mint operator+(mint p) {
        LL res = x + p.x;
        while (res >= mod) res -= mod;
        return res;
    }

    mint operator-(mint p) {
        LL res = x - p.x;
        while (res < 0) res += mod;
        return res;
    }

    mint operator*(mint p) {
        return x * p.x % mod;
    }

    mint pow(LL b) {
        mint res = 1, a = *this;
        while (b) {
            if (b & 1) res = res * a;
            a = a * a;
            b >>= 1;
        }
        return res;
    }

    mint operator/(mint p) {
        return (*this) * p.pow(mod - 2);
    }
};

// void solve() {
//     int m, k;
//     cin >> m >> k;
//
//     vector<mint> pre(k + 1);
//     pre[0] = 1;
//     for (int i = 1; i <= k; i++) {
//         mint sum = mint(i) * (m - 1) + 1;
//         pre[i] = pre[i - 1] * (sum - 1) / sum;
//     }
//
//     vector<mint> suf(k + 1);
//     suf[k] = 1;
//     for (int i = k - 1; i >= 1; i--) {
//         mint sum = mint(i) * (m - 1) + 1;
//         suf[i] = suf[i + 1] * (sum - 1) / sum;
//     }
//
//     vector<mint> p(k + 1);
//     mint pp = 1;
//     p[0] = p[1] = 0;
//     for (int i = 2; i <= k; i++) {
//         p[i] = p[i - 1];
//         p[i] = (mint(1) - p[i - 1]) * suf[i - 1] + p[i];
//     }
//
//     for (int i = 1; i <= k; i++) {
//         cout << p[i].x << " \n"[i == k];
//     }
//
//     mint sum = 1, ans = 0;
//     for (int i = 1; i <= k; i++) {
//         sum = sum * m;
//         ans = ans + sum * i * (mint(1) - p[i]);
//     }
//
//     cout << ans.x << "\n";
// }

void solve() {
    int m, k;
    cin >> m >> k;

    vector<mint> dep(k + 1);
    dep[0] = 0;
    mint ans = 0;

    vector<mint> sum(k + 1);
    sum[0] = 1;
    for (int i = 1; i <= k; i++) {
        sum[i] = sum[i - 1] + m - 1;
    }

    for (int i = 1; i <= k; i++) {
        //mint t = dep[i - 1] / sum[i - 1];
        dep[i] = //(t + 1) * m - dep[i - 1];
            (dep[i-1]*sum[i-1]+(dep[i-1]+1)*m-dep[i-1])/sum[i];
        ans = ans + (dep[i-1]+1)*m;
    }
    cout << ans.x << "\n";
}

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

    int t = 1;
    // cin >> t;
    while (t--) {
        solve();
    }
}
/*
 *2 3
333333337
500000005 666666673 750000007
0 750000007 375000004
17
*/