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

mod = 10**9+7
a = 1
b = 1
for i in range(min(K-2, N*M-K)):
    a *= (N*M-2-i)
    a %= mod
    b *= i+1
    b %= mod
comb = a * pow(b, mod - 2, mod) % mod

ans = (((N*N*(N+1)//2-N*(N+1)*(2*N+1)//6)*M**2 + (M*M*(M+1)//2-M*(M+1)*(2*M+1)//6)*N**2)*comb)%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 a = 1 a 1 b = 1 b 1 for i in range(min(K-2, N*M-K)):
    a *= (N*M-2-i)
    a %= mod
    b *= i+1
    b %= mod i range(min(K-2, N*M-K)) range min(K-2, N*M-K) min K-2 K 2 N*M-K N*M N M K a *= (N*M-2-i) a N*M-2-i N*M-2 N*M N M 2 i a %= mod a mod b *= i+1 b i+1 i 1 b %= mod b mod comb = a * pow(b, mod - 2, mod) % mod comb a * pow(b, mod - 2, mod) % mod a * pow(b, mod - 2, mod) a pow(b, mod - 2, mod) pow b mod - 2 mod 2 mod mod ans = (((N*N*(N+1)//2-N*(N+1)*(2*N+1)//6)*M**2 + (M*M*(M+1)//2-M*(M+1)*(2*M+1)//6)*N**2)*comb)%mod ans (((N*N*(N+1)//2-N*(N+1)*(2*N+1)//6)*M**2 + (M*M*(M+1)//2-M*(M+1)*(2*M+1)//6)*N**2)*comb)%mod ((N*N*(N+1)//2-N*(N+1)*(2*N+1)//6)*M**2 + (M*M*(M+1)//2-M*(M+1)*(2*M+1)//6)*N**2)*comb (N*N*(N+1)//2-N*(N+1)*(2*N+1)//6)*M**2 + (M*M*(M+1)//2-M*(M+1)*(2*M+1)//6)*N**2 (N*N*(N+1)//2-N*(N+1)*(2*N+1)//6)*M**2 N*N*(N+1)//2-N*(N+1)*(2*N+1)//6 N*N*(N+1)//2 N*N*(N+1) N*N N N 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 M**2 M 2 (M*M*(M+1)//2-M*(M+1)*(2*M+1)//6)*N**2 M*M*(M+1)//2-M*(M+1)*(2*M+1)//6 M*M*(M+1)//2 M*M*(M+1) M*M M M 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 N**2 N 2 comb mod print(ans) print ans b = 1 1 b K = map(int, input().split()) map(int, input().split()) K a %= mod mod a b *= i+1 i+1 b a *= (N*M-2-i) N*M-2-i a M,K = map(int, input().split()) map(int, input().split()) M N,M,K = map(int, input().split()) map(int, input().split()) N a = 1 1 a ans = (((N*N*(N+1)//2-N*(N+1)*(2*N+1)//6)*M**2 + (M*M*(M+1)//2-M*(M+1)*(2*M+1)//6)*N**2)*comb)%mod (((N*N*(N+1)//2-N*(N+1)*(2*N+1)//6)*M**2 + (M*M*(M+1)//2-M*(M+1)*(2*M+1)//6)*N**2)*comb)%mod ans comb = a * pow(b, mod - 2, mod) % mod a * pow(b, mod - 2, mod) % mod comb mod = 10**9+7 10**9+7 mod b %= mod mod b
N,M,K = map(int, input().split())

mod = 10**9+7
a = 1
b = 1
for i in range(min(K-2, N*M-K)):
    a *= (N*M-2-i)
    a %= mod
    b *= i+1
    b %= mod
comb = a * pow(b, mod - 2, mod) % mod

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