#include <bits/stdc++.h>
using namespace std;
long long int inv[3001], fact[3001], fact_inv[3001];
void offline() {
  inv[1] = 1;
  for (int i = 2; i < 3001; i++)
    inv[i] = 1000000007LL -
             ((1000000007LL / i) * inv[1000000007LL % i] % 1000000007LL);
  fact[0] = fact_inv[0] = 1;
  for (int i = 1; i < 3001; i++) fact[i] = (fact[i - 1] * i) % 1000000007LL;
  for (int i = 1; i < 3001; i++)
    fact_inv[i] = (fact_inv[i - 1] * inv[i]) % 1000000007LL;
}
long long int comb(int n, int r) {
  if (r < 0) return 0;
  if (n < r) return 0;
  return (((fact[n] * fact_inv[r]) % 1000000007LL) * fact_inv[n - r]) %
         1000000007LL;
}
map<pair<pair<int, int>, int>, long long int> myMap;
long long int func(int n, int m, int k) {
  if (k < 0) return 0;
  if ((n <= 0) || (m <= 0)) return k == 0;
  auto key = make_pair(make_pair(n, m), k);
  if (myMap.find(key) != myMap.end()) return myMap[key];
  long long int result = 0;
  for (int i = 1; i <= n; i++) {
    for (int j = 1; j <= m; j++) {
      long long int cnt = ((n + 1 - i) * (m + 1 - j)) % 1000000007LL;
      result = (result + (cnt * func(i - 2, j - 2, k - 1))) % 1000000007LL;
    }
  }
  return myMap[key] = result;
}
int main() {
  offline();
  int n, m, k;
  scanf("%d%d%d", &n, &m, &k);
  long long int result =
      (comb(n - 1, 2 * k) * comb(m - 1, 2 * k)) % 1000000007LL;
  printf("%lld\n", result);
  return 0;
}
