3 3 1
 300 2 300
 4000 4000 1
 4000 4000 100
 4 2 3
 3 300 300
 4000 100 4000
 3 300 1
 300 300 1
 4000 1000 3000
 3 2 4000
 4 3 2
 3 2 300
 4000 3998 2
 100 200 300
 300 300 300
 3 3 3
 10 10 10
 10 4 9
 239 300 231
 4000 2 3998
 200 100 300
 4000 2000 2000
 4000 4000 4000
 10 7 5
 4000 100 3900
 3 4000 4000

12
 775907030
 63263244
 994443885
 24
 496527918
 908339579
 828107078
 775907030
 876839920
 938379934
 48
 196174631
 296557186
 316471646
 375912430
 72
 318389383
 135283173
 774612666
 686088712
 949581532
 310481606
 997463324
 130636800
 221262673
 680114446

#include <bits/stdc++.h>
using namespace std;
const long long mod = 1000000009;
const int N = 4005;
long long fac[N], ifac[N], inv[N], pw2[N];
void pre() {
  fac[0] = ifac[0] = fac[1] = ifac[1] = inv[1] = 1;
  for (int i = 2; i < N; i++) {
    fac[i] = i * fac[i - 1] % mod;
    inv[i] = inv[mod % i] * (mod - mod / i) % mod;
    ifac[i] = ifac[i - 1] * inv[i] % mod;
  }
  pw2[0] = 1;
  for (int i = 1; i < N; i++) {
    pw2[i] = pw2[i - 1] << 1;
    if (pw2[i] >= mod) pw2[i] -= mod;
  }
}
int main() {
  ios::sync_with_stdio(false);
  cin.tie(NULL);
  cout.tie(NULL);
  cout << setprecision(32);
  pre();
  int n, w, b;
  cin >> n >> w >> b;
  long long ans = 0;
  for (int r = 1; r <= w; r++) {
    if (1 <= n - r && n - r <= b) {
      long long tmp = (fac[w - 1] * ifac[r - 1]) % mod;
      (tmp *= ifac[w - r]) %= mod;
      (tmp *= fac[b - 1]) %= mod;
      (tmp *= ifac[n - r - 1]) %= mod;
      (tmp *= ifac[b - n + r]) %= mod;
      (tmp *= (r - 1)) %= mod;
      ans += tmp;
      if (ans >= mod) ans -= mod;
    }
  }
  (ans *= fac[w]) %= mod;
  (ans *= fac[b]) %= mod;
  cout << ans << endl;
  return 0;
}
