n, m, k = map(int, input().split())

MOD = 1000000007
def pow_mod(x, y):
    if y == 0: return 1
    ans = 1
    while y > 1:
        if y % 2 != 0: 
            ans *= x
            ans %= MOD
        x *= x
        x %= MOD
        y //= 2
    return ans * x % MOD
mod_inv = lambda x: pow_mod(x, MOD - 2)

a = 0
for i in range(1, n):
  a += (i * (n - i)) % MOD 
  a %= MOD
a *= (m * m) % MOD
a %= MOD

b = 0
for i in range(1, m):
  b += (i * (m - i)) % MOD 
  b %= MOD
b *= (n * n) % MOD
b %= MOD

w = (m * n - 2) % MOD
u = 1
v = 1
for i in range(1, k-1):
  u *= (w - i + 1)
  u %= MOD
  v *= i
  v %= MOD

ans = ((a + b) % MOD) 
ans *= u
ans %= MOD
ans *= mod_inv(v)
ans %= MOD


print(ans) n, m, k = map(int, input().split()) n map(int, input().split()) map int input().split() ().split () input split m k MOD = 1000000007 MOD 1000000007 def pow_mod(x, y):
    if y == 0: return 1
    ans = 1
    while y > 1:
        if y % 2 != 0: 
            ans *= x
            ans %= MOD
        x *= x
        x %= MOD
        y //= 2
    return ans * x % MOD pow_mod if y == 0: return 1
     y == 0 y 0 ans = 1 ans 1 while y > 1:
        if y % 2 != 0: 
            ans *= x
            ans %= MOD
        x *= x
        x %= MOD
        y //= 2
     y > 1 y 1 if y % 2 != 0: 
            ans *= x
            ans %= MOD
         y % 2 != 0 y % 2 y 2 0 ans *= x ans x ans %= MOD ans MOD x *= x x x x %= MOD x MOD y //= 2 y 2 return ans * x % MOD ans * x % MOD ans * x ans x MOD x x y y mod_inv = lambda x: pow_mod(x, MOD - 2) mod_inv lambda x: pow_mod(x, MOD - 2) pow_mod(x, MOD - 2) pow_mod x MOD - 2 MOD 2 x a = 0 a 0 for i in range(1, n):
  a += (i * (n - i)) % MOD 
  a %= MOD i range(1, n) range 1 n a += (i * (n - i)) % MOD a (i * (n - i)) % MOD i * (n - i) i n - i n i MOD a %= MOD a MOD a *= (m * m) % MOD a (m * m) % MOD m * m m m MOD a %= MOD a MOD b = 0 b 0 for i in range(1, m):
  b += (i * (m - i)) % MOD 
  b %= MOD i range(1, m) range 1 m b += (i * (m - i)) % MOD b (i * (m - i)) % MOD i * (m - i) i m - i m i MOD b %= MOD b MOD b *= (n * n) % MOD b (n * n) % MOD n * n n n MOD b %= MOD b MOD w = (m * n - 2) % MOD w (m * n - 2) % MOD m * n - 2 m * n m n 2 MOD u = 1 u 1 v = 1 v 1 for i in range(1, k-1):
  u *= (w - i + 1)
  u %= MOD
  v *= i
  v %= MOD i range(1, k-1) range 1 k-1 k 1 u *= (w - i + 1) u w - i + 1 w - i w i 1 u %= MOD u MOD v *= i v i v %= MOD v MOD ans = ((a + b) % MOD) ans (a + b) % MOD a + b a b MOD ans *= u ans u ans %= MOD ans MOD ans *= mod_inv(v) ans mod_inv(v) mod_inv v ans %= MOD ans MOD print(ans) print ans ans *= u u ans mod_inv = lambda x: pow_mod(x, MOD - 2) lambda x: pow_mod(x, MOD - 2) mod_inv a %= MOD MOD a b = 0 0 b b %= MOD MOD b b *= (n * n) % MOD (n * n) % MOD b n, m, k = map(int, input().split()) map(int, input().split()) n ans %= MOD MOD ans ans %= MOD MOD ans u %= MOD MOD u a *= (m * m) % MOD (m * m) % MOD a w = (m * n - 2) % MOD (m * n - 2) % MOD w a = 0 0 a v = 1 1 v ans *= mod_inv(v) mod_inv(v) ans MOD = 1000000007 1000000007 MOD m, k = map(int, input().split()) map(int, input().split()) m a %= MOD MOD a v *= i i v a += (i * (n - i)) % MOD (i * (n - i)) % MOD a k = map(int, input().split()) map(int, input().split()) k u = 1 1 u v %= MOD MOD v b += (i * (m - i)) % MOD (i * (m - i)) % MOD b b %= MOD MOD b u *= (w - i + 1) w - i + 1 u ans = ((a + b) % MOD) (a + b) % MOD ans def pow_mod(x, y):
    if y == 0: return 1
    ans = 1
    while y > 1:
        if y % 2 != 0: 
            ans *= x
            ans %= MOD
        x *= x
        x %= MOD
        y //= 2
    return ans * x % MOD def pow_mod(x, y):
    if y == 0: return 1
    ans = 1
    while y > 1:
        if y % 2 != 0: 
            ans *= x
            ans %= MOD
        x *= x
        x %= MOD
        y //= 2
    return ans * x % MOD pow_mod
n, m, k = map(int, input().split())

MOD = 1000000007
def pow_mod(x, y):
    if y == 0: return 1
    ans = 1
    while y > 1:
        if y % 2 != 0: 
            ans *= x
            ans %= MOD
        x *= x
        x %= MOD
        y //= 2
    return ans * x % MOD
mod_inv = lambda x: pow_mod(x, MOD - 2)

a = 0
for i in range(1, n):
  a += (i * (n - i)) % MOD 
  a %= MOD
a *= (m * m) % MOD
a %= MOD

b = 0
for i in range(1, m):
  b += (i * (m - i)) % MOD 
  b %= MOD
b *= (n * n) % MOD
b %= MOD

w = (m * n - 2) % MOD
u = 1
v = 1
for i in range(1, k-1):
  u *= (w - i + 1)
  u %= MOD
  v *= i
  v %= MOD

ans = ((a + b) % MOD) 
ans *= u
ans %= MOD
ans *= mod_inv(v)
ans %= MOD


print(ans)