#include <bits/stdc++.h>
using namespace std;
const double eps = 1e-7;
const string nl = "\n";
template <typename T>
string convtos(T a) {
  stringstream ss;
  ss << a;
  return ss.str();
}
template <typename T>
T convtot(string s) {
  T x;
  stringstream ss;
  ss << s;
  ss >> x;
  return x;
}
void printtime() { cerr << double(clock()) / CLOCKS_PER_SEC << nl; }
void zfunc(const string &s, vector<int> &z) {
  int n = s.size();
  z.resize(n);
  z[0] = 0;
  for (int l = 0, r = 0, i = 1; i < n; ++i) {
    z[i] = (i <= r ? min(r - i + 1, z[i - l]) : 0);
    while (i + z[i] < n && s[z[i]] == s[z[i] + i]) z[i]++;
    if (i + z[i] > r) r = z[i] + i - 1, l = i;
  }
}
void pfunc(const string &s, vector<int> &p) {
  int n = s.size();
  p.resize(n);
  p[0] = 0;
  for (int i = 1; i < n; ++i) {
    int l = p[i - 1];
    while (l > 0 && s[l] != s[i]) l = p[l - 1];
    if (l == 0)
      p[i] = (s[i] == s[0]);
    else
      p[i] = l + 1;
  }
}
template <typename T>
T xgcd(T a, T b, T &x, T &y) {
  if (a == 0) {
    x = 0, y = 1;
    return b;
  }
  T x1, y1, ans;
  ans = xgcd(b % a, a, x1, y1);
  x = y1 - (b / a) * x1;
  y = x1;
  return ans;
}
void prime(int N, vector<int> &pr) {
  pr.resize(0);
  vector<int> lp(N, 0);
  for (size_t i = 2; i < N; ++i) {
    if (lp[i] == 0) {
      lp[i] = i;
      pr.push_back(i);
    } else
      for (size_t j = 0; j < pr.size() && pr[j] <= lp[i] && i * pr[j] < N; ++j)
        lp[i * pr[j]] = pr[j];
  }
}
void solve();
int main() {
  ios_base::sync_with_stdio(0);
  solve();
  return 0;
}
long long d[1100][1100], sum[1100];
long long p = 1e9 + 7;
int f(int n, int k) {
  if (k == 0) return 1;
  if (n <= 2) return 0;
  if (d[n][k] != -1) return d[n][k];
  long long res = 0, c;
  for (int l = n - 2; l >= 1; --l) {
    c = f(l, k - 1);
    c = (c * (n - 2 - l + 1)) % p;
    res = (res + c) % p;
  }
  d[n][k] = int(res);
  return d[n][k];
}
void solve() {
  int n, m, k;
  cin >> n >> m >> k;
  int i, j, z;
  for (i = 0; i <= k; ++i) d[0][i] = d[1][k] = 0;
  for (i = 0; i <= n; ++i) d[i][0] = 1;
  int N = max(n, m);
  sum[0] = 0;
  for (i = 1; i <= N; ++i) sum[i] = (sum[i - 1] + 1) % p;
  for (j = 1; j <= k; ++j) {
    for (i = 2; i <= N; ++i) {
      d[i][j] = (d[i - 1][j] + sum[i - 2]) % p;
    }
    sum[0] = 0;
    for (z = 1; z <= N; ++z) sum[z] = (sum[z - 1] + d[z][j]) % p;
  }
  cout << (d[n][k] * d[m][k]) % p;
}
