#include <bits/stdc++.h>
using namespace std;
template <typename M, typename N>
istream& operator>>(istream& is, pair<M, N>& p) {
  is >> p.first >> p.second;
  return is;
}
template <typename M>
istream& operator>>(istream& is, vector<M>& v) {
  for (auto& it : v) is >> it;
  return is;
}
template <typename M>
istream& operator>>(istream& is, vector<vector<M>>& v) {
  for (auto& it : v) is >> it;
  return is;
}
template <typename M, typename N>
ostream& operator<<(ostream& os, const pair<M, N>& p) {
  os << p.first << ' ' << p.second << '\n';
  return os;
}
template <typename M>
ostream& operator<<(ostream& os, const vector<M>& v) {
  for (auto it : v) os << it << ' ';
  cout << '\n';
  return os;
}
template <typename M>
ostream& operator<<(ostream& os, const vector<vector<M>>& v) {
  for (auto it : v) os << it;
  return os;
}
const long double pi = acos(-1);
const long long inf = (long long)2e9 + 10;
const int mod = 1000000007;
long long n, m;
const int N = 2e5 + 10;
long long binpow(long long a, long long b, long long m) {
  a %= m;
  long long res = 1;
  while (b > 0) {
    if (b & 1) res = res * a % m;
    a = a * a % m;
    b >>= 1;
  }
  return res;
}
vector<long long int> fact, fact_inverse;
template <typename type>
void tabulate(type n) {
  fact.resize(n + 1);
  fact_inverse.resize(n + 1);
  fact[0] = 1;
  fact_inverse[0] = binpow(1, mod - 2, mod);
  for (int i = 1; i < n + 1; i++) {
    fact[i] = fact[i - 1] * i;
    fact[i] %= mod;
    fact_inverse[i] = binpow(fact[i], mod - 2, mod);
  }
}
template <typename type>
long long int nCk(type n, type k) {
  return n < k ? 0
               : fact[n] * fact_inverse[n - k] % mod * fact_inverse[k] % mod;
}
template <typename type>
long long int nPk(type n, type k) {
  return n < k ? 0 : fact[n] * fact_inverse[n - k] % mod;
}
int main() {
  ios_base::sync_with_stdio(false);
  cin.tie(NULL);
  tabulate(4000);
  int n, m, k;
  cin >> n >> m >> k;
  n -= 2;
  m -= 2;
  cout << nCk(n + 1, 2 * k) * nCk(m + 1, 2 * k) % mod << '\n';
  return 0;
}
