mod = 10**9 + 7
 
 
def comb(n, k):
    x, y = 1, 1
    k = min(k, n - k)
    for i in range(k):
        x *= n - i
        x %= mod
        y *= i + 1
        y %= mod
    return x * pow(y, mod - 2, mod) % mod

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

c = comb(N*M-2,K-2)

ans = c * (M**2 * (N**2 * (N+1)//2 - N*(N+1)*(2*N+1)//6) + N**2 * (M**2 * (M+1)//2 - M*(M+1)*(2*M+1)//6)) % mod

print(ans) mod = 10**9 + 7 mod 10**9 + 7 10**9 10 9 7 def comb(n, k):
    x, y = 1, 1
    k = min(k, n - k)
    for i in range(k):
        x *= n - i
        x %= mod
        y *= i + 1
        y %= mod
    return x * pow(y, mod - 2, mod) % mod comb x, y = 1, 1 x 1 y 1 k = min(k, n - k) k min(k, n - k) min k n - k n k for i in range(k):
        x *= n - i
        x %= mod
        y *= i + 1
        y %= mod
     i range(k) range k x *= n - i x n - i n i x %= mod x mod y *= i + 1 y i + 1 i 1 y %= mod y mod return x * pow(y, mod - 2, mod) % mod x * pow(y, mod - 2, mod) % mod x * pow(y, mod - 2, mod) x pow(y, mod - 2, mod) pow y mod - 2 mod 2 mod mod n n k k N, M, K = map(int, input().split()) N map(int, input().split()) map int input().split() ().split () input split M K c = comb(N*M-2,K-2) c comb(N*M-2,K-2) comb N*M-2 N*M N M 2 K-2 K 2 ans = c * (M**2 * (N**2 * (N+1)//2 - N*(N+1)*(2*N+1)//6) + N**2 * (M**2 * (M+1)//2 - M*(M+1)*(2*M+1)//6)) % mod ans c * (M**2 * (N**2 * (N+1)//2 - N*(N+1)*(2*N+1)//6) + N**2 * (M**2 * (M+1)//2 - M*(M+1)*(2*M+1)//6)) % mod c * (M**2 * (N**2 * (N+1)//2 - N*(N+1)*(2*N+1)//6) + N**2 * (M**2 * (M+1)//2 - M*(M+1)*(2*M+1)//6)) c M**2 * (N**2 * (N+1)//2 - N*(N+1)*(2*N+1)//6) + N**2 * (M**2 * (M+1)//2 - M*(M+1)*(2*M+1)//6) M**2 * (N**2 * (N+1)//2 - N*(N+1)*(2*N+1)//6) M**2 M 2 N**2 * (N+1)//2 - N*(N+1)*(2*N+1)//6 N**2 * (N+1)//2 N**2 * (N+1) N**2 N 2 N+1 N 1 2 N*(N+1)*(2*N+1)//6 N*(N+1)*(2*N+1) N*(N+1) N N+1 N 1 2*N+1 2*N 2 N 1 6 N**2 * (M**2 * (M+1)//2 - M*(M+1)*(2*M+1)//6) N**2 N 2 M**2 * (M+1)//2 - M*(M+1)*(2*M+1)//6 M**2 * (M+1)//2 M**2 * (M+1) M**2 M 2 M+1 M 1 2 M*(M+1)*(2*M+1)//6 M*(M+1)*(2*M+1) M*(M+1) M M+1 M 1 2*M+1 2*M 2 M 1 6 mod print(ans) print ans ans = c * (M**2 * (N**2 * (N+1)//2 - N*(N+1)*(2*N+1)//6) + N**2 * (M**2 * (M+1)//2 - M*(M+1)*(2*M+1)//6)) % mod c * (M**2 * (N**2 * (N+1)//2 - N*(N+1)*(2*N+1)//6) + N**2 * (M**2 * (M+1)//2 - M*(M+1)*(2*M+1)//6)) % mod ans mod = 10**9 + 7 10**9 + 7 mod def comb(n, k):
    x, y = 1, 1
    k = min(k, n - k)
    for i in range(k):
        x *= n - i
        x %= mod
        y *= i + 1
        y %= mod
    return x * pow(y, mod - 2, mod) % mod def comb(n, k):
    x, y = 1, 1
    k = min(k, n - k)
    for i in range(k):
        x *= n - i
        x %= mod
        y *= i + 1
        y %= mod
    return x * pow(y, mod - 2, mod) % mod comb c = comb(N*M-2,K-2) comb(N*M-2,K-2) c M, K = map(int, input().split()) map(int, input().split()) M K = map(int, input().split()) map(int, input().split()) K N, M, K = map(int, input().split()) map(int, input().split()) N
mod = 10**9 + 7
 
 
def comb(n, k):
    x, y = 1, 1
    k = min(k, n - k)
    for i in range(k):
        x *= n - i
        x %= mod
        y *= i + 1
        y %= mod
    return x * pow(y, mod - 2, mod) % mod

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

c = comb(N*M-2,K-2)

ans = c * (M**2 * (N**2 * (N+1)//2 - N*(N+1)*(2*N+1)//6) + N**2 * (M**2 * (M+1)//2 - M*(M+1)*(2*M+1)//6)) % mod

print(ans)