N,M,K=map(int,input().split())
c=[0 for i in range(N*M)]
c[1]=1
c[0]=1
mod=10**9+7
for i in range(1,N*M-1):
  c[i+1]=(i+1)*c[i]%mod
def comb(a,b):
  return c[a]*pow(c[b],mod-2,mod)*pow(c[a-b],mod-2,mod)%mod
ans=0
for i in range(N):
  ans+=i*(N-i)*pow(M,2)
  ans%=mod
for i in range(M):
  ans+=i*(M-i)*pow(N,2)%mod
  ans%=mod
ans*=comb(N*M-2,K-2)

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

ans%=mod
print(ans)