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

mod = 10**9+7
fact = [1]*(2*N*M+1)
for n in range(1,2*N*M+1):
    fact[n] = n*fact[n-1]%mod

def inv(n):
    n %= mod
    return pow(n,mod-2,mod)

ans = M*(M**2-1)*N**2%mod+N*(N**2-1)*M**2%mod
ans *= inv(6)*fact[N*M-2]*inv(fact[K-2]*fact[N*M-K])
ans %= 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 fact = [1]*(2*N*M+1) fact [1]*(2*N*M+1) [1] 1 2*N*M+1 2*N*M 2*N 2 N M 1 for n in range(1,2*N*M+1):
    fact[n] = n*fact[n-1]%mod n range(1,2*N*M+1) range 1 2*N*M+1 2*N*M 2*N 2 N M 1 fact[n] = n*fact[n-1]%mod fact[n] fact n n*fact[n-1]%mod n*fact[n-1] n fact[n-1] fact n-1 n 1 mod def inv(n):
    n %= mod
    return pow(n,mod-2,mod) inv n %= mod n mod return pow(n,mod-2,mod) pow(n,mod-2,mod) pow n mod-2 mod 2 mod n n ans = M*(M**2-1)*N**2%mod+N*(N**2-1)*M**2%mod ans M*(M**2-1)*N**2%mod+N*(N**2-1)*M**2%mod M*(M**2-1)*N**2%mod M*(M**2-1)*N**2 M*(M**2-1) M M**2-1 M**2 M 2 1 N**2 N 2 mod N*(N**2-1)*M**2%mod N*(N**2-1)*M**2 N*(N**2-1) N N**2-1 N**2 N 2 1 M**2 M 2 mod ans *= inv(6)*fact[N*M-2]*inv(fact[K-2]*fact[N*M-K]) ans inv(6)*fact[N*M-2]*inv(fact[K-2]*fact[N*M-K]) inv(6)*fact[N*M-2] inv(6) inv 6 fact[N*M-2] fact N*M-2 N*M N M 2 inv(fact[K-2]*fact[N*M-K]) inv fact[K-2]*fact[N*M-K] fact[K-2] fact K-2 K 2 fact[N*M-K] fact N*M-K N*M N M K ans %= mod ans mod print(ans) print ans N,M,K = map(int,input().split()) map(int,input().split()) N mod = 10**9+7 10**9+7 mod ans *= inv(6)*fact[N*M-2]*inv(fact[K-2]*fact[N*M-K]) inv(6)*fact[N*M-2]*inv(fact[K-2]*fact[N*M-K]) ans K = map(int,input().split()) map(int,input().split()) K ans %= mod mod ans M,K = map(int,input().split()) map(int,input().split()) M ans = M*(M**2-1)*N**2%mod+N*(N**2-1)*M**2%mod M*(M**2-1)*N**2%mod+N*(N**2-1)*M**2%mod ans fact = [1]*(2*N*M+1) [1]*(2*N*M+1) fact def inv(n):
    n %= mod
    return pow(n,mod-2,mod) def inv(n):
    n %= mod
    return pow(n,mod-2,mod) inv
N,M,K = map(int,input().split())

mod = 10**9+7
fact = [1]*(2*N*M+1)
for n in range(1,2*N*M+1):
    fact[n] = n*fact[n-1]%mod

def inv(n):
    n %= mod
    return pow(n,mod-2,mod)

ans = M*(M**2-1)*N**2%mod+N*(N**2-1)*M**2%mod
ans *= inv(6)*fact[N*M-2]*inv(fact[K-2]*fact[N*M-K])
ans %= mod

print(ans)