from collections import defaultdict

mod = 10 ** 9 + 7

memo = defaultdict(int)


def cmb(n, r):
    if (n, r) in memo:
        return memo[(n, r)]

    if n < r or r < 0:
        return 0

    if r > n - r:
        return cmb(n, n - r)

    res = f[n] * pow(f[n - r], mod - 2, mod) * pow(f[r], mod - 2, mod) % mod
    memo[(n, r)] = res
    return res


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

f = [1]
for i in range(1, 2 * 10 ** 5 + 1):
    f.append(f[-1] * i % mod)

ans = 0
for d in range(1, m):
    ans += d * (m - d) * n ** 2
    ans %= mod

for d in range(1, n):
    ans += d * (n - d) * m ** 2
    ans %= mod

ans *= cmb(n * m - 2, k - 2)
ans %= mod

print(ans) from collections import defaultdict mod = 10 ** 9 + 7 mod 10 ** 9 + 7 10 ** 9 10 9 7 memo = defaultdict(int) memo defaultdict(int) defaultdict int def cmb(n, r):
    if (n, r) in memo:
        return memo[(n, r)]

    if n < r or r < 0:
        return 0

    if r > n - r:
        return cmb(n, n - r)

    res = f[n] * pow(f[n - r], mod - 2, mod) * pow(f[r], mod - 2, mod) % mod
    memo[(n, r)] = res
    return res cmb if (n, r) in memo:
        return memo[(n, r)]

     (n, r) in memo (n, r) n r memo return memo[(n, r)] memo[(n, r)] memo (n, r) n r if n < r or r < 0:
        return 0

     n < r or r < 0 n < r n r r < 0 r 0 return 0 0 if r > n - r:
        return cmb(n, n - r)

     r > n - r r n - r n r return cmb(n, n - r) cmb(n, n - r) cmb n n - r n r res = f[n] * pow(f[n - r], mod - 2, mod) * pow(f[r], mod - 2, mod) % mod res f[n] * pow(f[n - r], mod - 2, mod) * pow(f[r], mod - 2, mod) % mod f[n] * pow(f[n - r], mod - 2, mod) * pow(f[r], mod - 2, mod) f[n] * pow(f[n - r], mod - 2, mod) f[n] f n pow(f[n - r], mod - 2, mod) pow f[n - r] f n - r n r mod - 2 mod 2 mod pow(f[r], mod - 2, mod) pow f[r] f r mod - 2 mod 2 mod mod memo[(n, r)] = res memo[(n, r)] memo (n, r) n r res return res res n n r r n, m, k = map(int, input().split()) n map(int, input().split()) map int input().split() ().split () input split m k f = [1] f [1] 1 for i in range(1, 2 * 10 ** 5 + 1):
    f.append(f[-1] * i % mod) i range(1, 2 * 10 ** 5 + 1) range 1 2 * 10 ** 5 + 1 2 * 10 ** 5 2 10 ** 5 10 5 1 f.append(f[-1] * i % mod) f.append f append f[-1] * i % mod f[-1] * i f[-1] f -1 i mod ans = 0 ans 0 for d in range(1, m):
    ans += d * (m - d) * n ** 2
    ans %= mod d range(1, m) range 1 m ans += d * (m - d) * n ** 2 ans d * (m - d) * n ** 2 d * (m - d) d m - d m d n ** 2 n 2 ans %= mod ans mod for d in range(1, n):
    ans += d * (n - d) * m ** 2
    ans %= mod d range(1, n) range 1 n ans += d * (n - d) * m ** 2 ans d * (n - d) * m ** 2 d * (n - d) d n - d n d m ** 2 m 2 ans %= mod ans mod ans *= cmb(n * m - 2, k - 2) ans cmb(n * m - 2, k - 2) cmb n * m - 2 n * m n m 2 k - 2 k 2 ans %= mod ans mod print(ans) print ans ans %= mod mod ans ans *= cmb(n * m - 2, k - 2) cmb(n * m - 2, k - 2) ans k = map(int, input().split()) map(int, input().split()) k ans = 0 0 ans ans %= mod mod ans def cmb(n, r):
    if (n, r) in memo:
        return memo[(n, r)]

    if n < r or r < 0:
        return 0

    if r > n - r:
        return cmb(n, n - r)

    res = f[n] * pow(f[n - r], mod - 2, mod) * pow(f[r], mod - 2, mod) % mod
    memo[(n, r)] = res
    return res def cmb(n, r):
    if (n, r) in memo:
        return memo[(n, r)]

    if n < r or r < 0:
        return 0

    if r > n - r:
        return cmb(n, n - r)

    res = f[n] * pow(f[n - r], mod - 2, mod) * pow(f[r], mod - 2, mod) % mod
    memo[(n, r)] = res
    return res cmb mod = 10 ** 9 + 7 10 ** 9 + 7 mod memo = defaultdict(int) defaultdict(int) memo ans += d * (m - d) * n ** 2 d * (m - d) * n ** 2 ans n, m, k = map(int, input().split()) map(int, input().split()) n f = [1] [1] f ans += d * (n - d) * m ** 2 d * (n - d) * m ** 2 ans m, k = map(int, input().split()) map(int, input().split()) m ans %= mod mod ans
from collections import defaultdict

mod = 10 ** 9 + 7

memo = defaultdict(int)


def cmb(n, r):
    if (n, r) in memo:
        return memo[(n, r)]

    if n < r or r < 0:
        return 0

    if r > n - r:
        return cmb(n, n - r)

    res = f[n] * pow(f[n - r], mod - 2, mod) * pow(f[r], mod - 2, mod) % mod
    memo[(n, r)] = res
    return res


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

f = [1]
for i in range(1, 2 * 10 ** 5 + 1):
    f.append(f[-1] * i % mod)

ans = 0
for d in range(1, m):
    ans += d * (m - d) * n ** 2
    ans %= mod

for d in range(1, n):
    ans += d * (n - d) * m ** 2
    ans %= mod

ans *= cmb(n * m - 2, k - 2)
ans %= mod

print(ans)
