MOD = 10**9 + 7

N, M, K = map(int, input().split())

def getComb(n, k, MOD):
    if n < k:
        return 0
    if n-k < k:
        k = n-k
    comb = 1
    for x in range(n-k+1, n+1):
        comb = (comb * x) % MOD
    d = 1
    for x in range(1, k+1):
        d = (d * x) % MOD
    comb *= pow(d, MOD-2, MOD)
    return comb % MOD

comb = getComb(N*M-2, K-2, MOD)

ans = 0
for i in range(N):
    ans += i*(N-i)*M*M
    ans %= MOD

for j in range(M):
    ans += j*(M-j)*N*N
    ans %= MOD

ans *= comb
ans %= MOD

print(ans) MOD = 10**9 + 7 MOD 10**9 + 7 10**9 10 9 7 N, M, K = map(int, input().split()) N map(int, input().split()) map int input().split() ().split () input split M K def getComb(n, k, MOD):
    if n < k:
        return 0
    if n-k < k:
        k = n-k
    comb = 1
    for x in range(n-k+1, n+1):
        comb = (comb * x) % MOD
    d = 1
    for x in range(1, k+1):
        d = (d * x) % MOD
    comb *= pow(d, MOD-2, MOD)
    return comb % MOD getComb if n < k:
        return 0
     n < k n k return 0 0 if n-k < k:
        k = n-k
     n-k < k n-k n k k k = n-k k n-k n k comb = 1 comb 1 for x in range(n-k+1, n+1):
        comb = (comb * x) % MOD
     x range(n-k+1, n+1) range n-k+1 n-k n k 1 n+1 n 1 comb = (comb * x) % MOD comb (comb * x) % MOD comb * x comb x MOD d = 1 d 1 for x in range(1, k+1):
        d = (d * x) % MOD
     x range(1, k+1) range 1 k+1 k 1 d = (d * x) % MOD d (d * x) % MOD d * x d x MOD comb *= pow(d, MOD-2, MOD) comb pow(d, MOD-2, MOD) pow d MOD-2 MOD 2 MOD return comb % MOD comb % MOD comb MOD n n k k MOD MOD comb = getComb(N*M-2, K-2, MOD) comb getComb(N*M-2, K-2, MOD) getComb N*M-2 N*M N M 2 K-2 K 2 MOD ans = 0 ans 0 for i in range(N):
    ans += i*(N-i)*M*M
    ans %= MOD i range(N) range N ans += i*(N-i)*M*M ans i*(N-i)*M*M i*(N-i)*M i*(N-i) i N-i N i M M ans %= MOD ans MOD for j in range(M):
    ans += j*(M-j)*N*N
    ans %= MOD j range(M) range M ans += j*(M-j)*N*N ans j*(M-j)*N*N j*(M-j)*N j*(M-j) j M-j M j N N ans %= MOD ans MOD ans *= comb ans comb ans %= MOD ans MOD print(ans) print ans ans = 0 0 ans ans *= comb comb ans ans %= MOD MOD ans def getComb(n, k, MOD):
    if n < k:
        return 0
    if n-k < k:
        k = n-k
    comb = 1
    for x in range(n-k+1, n+1):
        comb = (comb * x) % MOD
    d = 1
    for x in range(1, k+1):
        d = (d * x) % MOD
    comb *= pow(d, MOD-2, MOD)
    return comb % MOD def getComb(n, k, MOD):
    if n < k:
        return 0
    if n-k < k:
        k = n-k
    comb = 1
    for x in range(n-k+1, n+1):
        comb = (comb * x) % MOD
    d = 1
    for x in range(1, k+1):
        d = (d * x) % MOD
    comb *= pow(d, MOD-2, MOD)
    return comb % MOD getComb K = map(int, input().split()) map(int, input().split()) K ans += i*(N-i)*M*M i*(N-i)*M*M ans M, K = map(int, input().split()) map(int, input().split()) M ans %= MOD MOD ans N, M, K = map(int, input().split()) map(int, input().split()) N ans %= MOD MOD ans ans += j*(M-j)*N*N j*(M-j)*N*N ans comb = getComb(N*M-2, K-2, MOD) getComb(N*M-2, K-2, MOD) comb MOD = 10**9 + 7 10**9 + 7 MOD
MOD = 10**9 + 7

N, M, K = map(int, input().split())

def getComb(n, k, MOD):
    if n < k:
        return 0
    if n-k < k:
        k = n-k
    comb = 1
    for x in range(n-k+1, n+1):
        comb = (comb * x) % MOD
    d = 1
    for x in range(1, k+1):
        d = (d * x) % MOD
    comb *= pow(d, MOD-2, MOD)
    return comb % MOD

comb = getComb(N*M-2, K-2, MOD)

ans = 0
for i in range(N):
    ans += i*(N-i)*M*M
    ans %= MOD

for j in range(M):
    ans += j*(M-j)*N*N
    ans %= MOD

ans *= comb
ans %= MOD

print(ans)
