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

mod = 10**9 + 7

# (n * m - 2)C(k - 2)
# 距離xの組 -  n - x + 1
# 横 n - x + 1
# 縦 m * m
ans = 0
for y in range(1, n):
    ans += y * (n - y) * m * m % mod
for x in range(1, m):
    ans += x * (m - x) * n * n % mod


def inv(x):
    return pow(x, mod - 2, mod)


# comb
for i in range(k - 2):
    ans = ans * (n * m - 2 - i) * inv(i + 1) % 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 = 10**9 + 7 mod 10**9 + 7 10**9 10 9 7 ans = 0 ans 0 for y in range(1, n):
    ans += y * (n - y) * m * m % mod y range(1, n) range 1 n ans += y * (n - y) * m * m % mod ans y * (n - y) * m * m % mod y * (n - y) * m * m y * (n - y) * m y * (n - y) y n - y n y m m mod for x in range(1, m):
    ans += x * (m - x) * n * n % mod x range(1, m) range 1 m ans += x * (m - x) * n * n % mod ans x * (m - x) * n * n % mod x * (m - x) * n * n x * (m - x) * n x * (m - x) x m - x m x n n mod def inv(x):
    return pow(x, mod - 2, mod)


# comb inv return pow(x, mod - 2, mod) pow(x, mod - 2, mod) pow x mod - 2 mod 2 mod x x for i in range(k - 2):
    ans = ans * (n * m - 2 - i) * inv(i + 1) % mod i range(k - 2) range k - 2 k 2 ans = ans * (n * m - 2 - i) * inv(i + 1) % mod ans ans * (n * m - 2 - i) * inv(i + 1) % mod ans * (n * m - 2 - i) * inv(i + 1) ans * (n * m - 2 - i) ans n * m - 2 - i n * m - 2 n * m n m 2 i inv(i + 1) inv i + 1 i 1 mod print(ans) print ans k = map(int, input().split()) map(int, input().split()) k ans = 0 0 ans ans += y * (n - y) * m * m % mod y * (n - y) * m * m % mod ans n, m, k = map(int, input().split()) map(int, input().split()) n mod = 10**9 + 7 10**9 + 7 mod def inv(x):
    return pow(x, mod - 2, mod)


# comb def inv(x):
    return pow(x, mod - 2, mod)


# comb inv ans += x * (m - x) * n * n % mod x * (m - x) * n * n % mod ans m, k = map(int, input().split()) map(int, input().split()) m ans = ans * (n * m - 2 - i) * inv(i + 1) % mod ans * (n * m - 2 - i) * inv(i + 1) % mod ans
n, m, k = map(int, input().split())

mod = 10**9 + 7

# (n * m - 2)C(k - 2)
# 距離xの組 -  n - x + 1
# 横 n - x + 1
# 縦 m * m
ans = 0
for y in range(1, n):
    ans += y * (n - y) * m * m % mod
for x in range(1, m):
    ans += x * (m - x) * n * n % mod


def inv(x):
    return pow(x, mod - 2, mod)


# comb
for i in range(k - 2):
    ans = ans * (n * m - 2 - i) * inv(i + 1) % mod

print(ans)
