#include <bits/stdc++.h>
using namespace std;
template <class T>
inline T checkmin(T &a, T b) {
  return (a < b) ? a : a = b;
}
template <class T>
inline T checkmax(T &a, T b) {
  return (a > b) ? a : a = b;
}
template <class T>
T GCD(T a, T b) {
  if (a < 0) return GCD(-a, b);
  if (b < 0) return GCD(a, -b);
  return (a == 0) ? b : GCD(b % a, a);
}
template <class T>
T LCM(T a, T b) {
  if (a < 0) return LCM(-a, b);
  if (b < 0) return LCM(a, -b);
  return (a == 0 || b == 0) ? 0 : a / GCD(a, b) * b;
}
template <class T>
T ExtGCD(T a, T b, T &x, T &y) {
  if (a < 0) {
    T c = ExtGCD(-a, b, x, y);
    x = -x;
    return c;
  }
  if (b < 0) {
    T c = ExtGCD(a, -b, x, y);
    y = -y;
    return c;
  }
  if (a == 0) {
    x = 0, y = 1;
    return b;
  } else {
    T c = ExtGCD(b % a, a, y, x);
    x -= b / a * y;
    return c;
  }
}
template <class T>
inline T sqr(T X) {
  return X * X;
}
namespace Poor {
const int MaxiN = 1005;
const int ModNum = 1000000007;
int N, M, K;
int Opt[MaxiN][MaxiN], S1[MaxiN][MaxiN], S2[MaxiN][MaxiN];
inline int Add(int a, int b) {
  if (a + b < ModNum) return a + b;
  return a + b - ModNum;
}
inline int Mul(int a, int b) { return ((long long)a * (long long)b) % ModNum; }
void Pretreat() {
  Opt[1][0] = 0;
  for (int i = 1; i < MaxiN; ++i) Opt[1][i] = 1;
  S1[1][0] = 0;
  for (int i = 1; i < MaxiN; ++i) S1[1][i] = Add(S1[1][i - 1], Opt[1][i]);
  S2[1][0] = 0;
  for (int i = 1; i < MaxiN; ++i) S2[1][i] = Add(S2[1][i - 1], S1[1][i]);
  for (int k = 2; k < MaxiN; ++k) {
    Opt[k][0] = 0;
    S1[k][0] = 0;
    S2[k][0] = 0;
    for (int i = 1; i < MaxiN; ++i) {
      Opt[k][i] = (i - 2 < 0) ? 0 : S2[k - 1][i - 2];
      S1[k][i] = Add(S1[k][i - 1], Opt[k][i]);
      S2[k][i] = Add(S2[k][i - 1], S1[k][i]);
    }
  }
}
void Run() {
  Pretreat();
  cin >> N >> M >> K;
  cout << Mul(Opt[K + 1][N], Opt[K + 1][M]) << endl;
}
}  // namespace Poor
int main() {
  Poor::Run();
  return 0;
}
