N,M,K=map(int,input().split())
L=N*M+10
mod=10**9+7
table=[1]*(L+3)
t=1
for i in range(1,L+3):
    t*=i
    t%=mod
    table[i]=t
rtable=[1]*(L+3)
t=1
for i in range(1,L+3):
    t*=pow(i,mod-2,mod)
    t%=mod
    rtable[i]=t
ans=0
t=table[N*M-2]*rtable[K-2]*rtable[N*M-K]
for i in range(1,M):
    ans+=N*N*(M-i)*i
    ans%=mod
for i in range(1,N):
    ans+=M*M*(N-i)*i
    ans%=mod
print((ans*t)%mod) N,M,K=map(int,input().split()) N map(int,input().split()) map int input().split() ().split () input split M K L=N*M+10 L N*M+10 N*M N M 10 mod=10**9+7 mod 10**9+7 10**9 10 9 7 table=[1]*(L+3) table [1]*(L+3) [1] 1 L+3 L 3 t=1 t 1 for i in range(1,L+3):
    t*=i
    t%=mod
    table[i]=t i range(1,L+3) range 1 L+3 L 3 t*=i t i t%=mod t mod table[i]=t table[i] table i t rtable=[1]*(L+3) rtable [1]*(L+3) [1] 1 L+3 L 3 t=1 t 1 for i in range(1,L+3):
    t*=pow(i,mod-2,mod)
    t%=mod
    rtable[i]=t i range(1,L+3) range 1 L+3 L 3 t*=pow(i,mod-2,mod) t pow(i,mod-2,mod) pow i mod-2 mod 2 mod t%=mod t mod rtable[i]=t rtable[i] rtable i t ans=0 ans 0 t=table[N*M-2]*rtable[K-2]*rtable[N*M-K] t table[N*M-2]*rtable[K-2]*rtable[N*M-K] table[N*M-2]*rtable[K-2] table[N*M-2] table N*M-2 N*M N M 2 rtable[K-2] rtable K-2 K 2 rtable[N*M-K] rtable N*M-K N*M N M K for i in range(1,M):
    ans+=N*N*(M-i)*i
    ans%=mod i range(1,M) range 1 M ans+=N*N*(M-i)*i ans N*N*(M-i)*i N*N*(M-i) N*N N N M-i M i i ans%=mod ans mod for i in range(1,N):
    ans+=M*M*(N-i)*i
    ans%=mod i range(1,N) range 1 N ans+=M*M*(N-i)*i ans M*M*(N-i)*i M*M*(N-i) M*M M M N-i N i i ans%=mod ans mod print((ans*t)%mod) print (ans*t)%mod ans*t ans t mod table=[1]*(L+3) [1]*(L+3) table t=1 1 t N,M,K=map(int,input().split()) map(int,input().split()) N ans=0 0 ans t=1 1 t L=N*M+10 N*M+10 L ans%=mod mod ans t%=mod mod t t*=pow(i,mod-2,mod) pow(i,mod-2,mod) t K=map(int,input().split()) map(int,input().split()) K ans+=N*N*(M-i)*i N*N*(M-i)*i ans rtable=[1]*(L+3) [1]*(L+3) rtable M,K=map(int,input().split()) map(int,input().split()) M ans+=M*M*(N-i)*i M*M*(N-i)*i ans mod=10**9+7 10**9+7 mod t*=i i t ans%=mod mod ans t=table[N*M-2]*rtable[K-2]*rtable[N*M-K] table[N*M-2]*rtable[K-2]*rtable[N*M-K] t t%=mod mod t
N,M,K=map(int,input().split())
L=N*M+10
mod=10**9+7
table=[1]*(L+3)
t=1
for i in range(1,L+3):
    t*=i
    t%=mod
    table[i]=t
rtable=[1]*(L+3)
t=1
for i in range(1,L+3):
    t*=pow(i,mod-2,mod)
    t%=mod
    rtable[i]=t
ans=0
t=table[N*M-2]*rtable[K-2]*rtable[N*M-K]
for i in range(1,M):
    ans+=N*N*(M-i)*i
    ans%=mod
for i in range(1,N):
    ans+=M*M*(N-i)*i
    ans%=mod
print((ans*t)%mod)