N,M,K=map(int,input().split())
def cmb(x,y):
  X=1
  Y=1
  for i in range(1,y+1):
    X= X*(x+1-i)%MOD
    Y= Y*i%MOD
  return (X*pow(Y,MOD-2,MOD))%MOD
MOD=10**9+7
A=cmb(M*N-2, K-2)
ans=0
for d in range(1,N):
  ans+=d*(N-d)*M*M*A%MOD
for d in range(1,M):
  ans+=d*(M-d)*N*N*A%MOD
print(ans%MOD) N,M,K=map(int,input().split()) N map(int,input().split()) map int input().split() ().split () input split M K def cmb(x,y):
  X=1
  Y=1
  for i in range(1,y+1):
    X= X*(x+1-i)%MOD
    Y= Y*i%MOD
  return (X*pow(Y,MOD-2,MOD))%MOD cmb X=1 X 1 Y=1 Y 1 for i in range(1,y+1):
    X= X*(x+1-i)%MOD
    Y= Y*i%MOD
   i range(1,y+1) range 1 y+1 y 1 X= X*(x+1-i)%MOD X X*(x+1-i)%MOD X*(x+1-i) X x+1-i x+1 x 1 i MOD Y= Y*i%MOD Y Y*i%MOD Y*i Y i 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 x x y y MOD=10**9+7 MOD 10**9+7 10**9 10 9 7 A=cmb(M*N-2, K-2) A cmb(M*N-2, K-2) cmb M*N-2 M*N M N 2 K-2 K 2 ans=0 ans 0 for d in range(1,N):
  ans+=d*(N-d)*M*M*A%MOD d range(1,N) range 1 N ans+=d*(N-d)*M*M*A%MOD ans d*(N-d)*M*M*A%MOD d*(N-d)*M*M*A d*(N-d)*M*M d*(N-d)*M d*(N-d) d N-d N d M M A MOD for d in range(1,M):
  ans+=d*(M-d)*N*N*A%MOD d range(1,M) range 1 M ans+=d*(M-d)*N*N*A%MOD ans d*(M-d)*N*N*A%MOD d*(M-d)*N*N*A d*(M-d)*N*N d*(M-d)*N d*(M-d) d M-d M d N N A MOD print(ans%MOD) print ans%MOD ans MOD ans+=d*(N-d)*M*M*A%MOD d*(N-d)*M*M*A%MOD ans A=cmb(M*N-2, K-2) cmb(M*N-2, K-2) A ans+=d*(M-d)*N*N*A%MOD d*(M-d)*N*N*A%MOD ans ans=0 0 ans N,M,K=map(int,input().split()) map(int,input().split()) N K=map(int,input().split()) map(int,input().split()) K def cmb(x,y):
  X=1
  Y=1
  for i in range(1,y+1):
    X= X*(x+1-i)%MOD
    Y= Y*i%MOD
  return (X*pow(Y,MOD-2,MOD))%MOD def cmb(x,y):
  X=1
  Y=1
  for i in range(1,y+1):
    X= X*(x+1-i)%MOD
    Y= Y*i%MOD
  return (X*pow(Y,MOD-2,MOD))%MOD cmb 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())
def cmb(x,y):
  X=1
  Y=1
  for i in range(1,y+1):
    X= X*(x+1-i)%MOD
    Y= Y*i%MOD
  return (X*pow(Y,MOD-2,MOD))%MOD
MOD=10**9+7
A=cmb(M*N-2, K-2)
ans=0
for d in range(1,N):
  ans+=d*(N-d)*M*M*A%MOD
for d in range(1,M):
  ans+=d*(M-d)*N*N*A%MOD
print(ans%MOD)