MOD=10**9+7
UPPERLIMIT=3*10**5
MODMUL=[1, 1]+[0]*(UPPERLIMIT-1)
for i in range(2, UPPERLIMIT+1):
  MODMUL[i]=MODMUL[i-1]*i%MOD
MODDIV=[1]*UPPERLIMIT+[pow(MODMUL[-1], MOD-2, MOD)]
for i in range(UPPERLIMIT, 0, -1):
  MODDIV[i-1]=MODDIV[i]*i%MOD
def MODCOMB(n, r):
  if n<r or n<1:
    return 0
  else:
    return (((MODMUL[n]*MODDIV[r])%MOD)*MODDIV[n-r])%MOD

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

ans=0
k=MODCOMB(N*M-2, K-2)
for sa in range(1, N):
  ans=(ans+sa*(N-sa)*M*M*k)%MOD
    
for sa in range(1, M):
  ans=(ans+sa*(M-sa)*N*N*k)%MOD
    
print(ans) MOD=10**9+7 MOD 10**9+7 10**9 10 9 7 UPPERLIMIT=3*10**5 UPPERLIMIT 3*10**5 3 10**5 10 5 MODMUL=[1, 1]+[0]*(UPPERLIMIT-1) MODMUL [1, 1]+[0]*(UPPERLIMIT-1) [1, 1] 1 1 [0]*(UPPERLIMIT-1) [0] 0 UPPERLIMIT-1 UPPERLIMIT 1 for i in range(2, UPPERLIMIT+1):
  MODMUL[i]=MODMUL[i-1]*i%MOD i range(2, UPPERLIMIT+1) range 2 UPPERLIMIT+1 UPPERLIMIT 1 MODMUL[i]=MODMUL[i-1]*i%MOD MODMUL[i] MODMUL i MODMUL[i-1]*i%MOD MODMUL[i-1]*i MODMUL[i-1] MODMUL i-1 i 1 i MOD MODDIV=[1]*UPPERLIMIT+[pow(MODMUL[-1], MOD-2, MOD)] MODDIV [1]*UPPERLIMIT+[pow(MODMUL[-1], MOD-2, MOD)] [1]*UPPERLIMIT [1] 1 UPPERLIMIT [pow(MODMUL[-1], MOD-2, MOD)] pow(MODMUL[-1], MOD-2, MOD) pow MODMUL[-1] MODMUL -1 MOD-2 MOD 2 MOD for i in range(UPPERLIMIT, 0, -1):
  MODDIV[i-1]=MODDIV[i]*i%MOD i range(UPPERLIMIT, 0, -1) range UPPERLIMIT 0 -1 MODDIV[i-1]=MODDIV[i]*i%MOD MODDIV[i-1] MODDIV i-1 i 1 MODDIV[i]*i%MOD MODDIV[i]*i MODDIV[i] MODDIV i i MOD def MODCOMB(n, r):
  if n<r or n<1:
    return 0
  else:
    return (((MODMUL[n]*MODDIV[r])%MOD)*MODDIV[n-r])%MOD MODCOMB if n<r or n<1:
    return 0
  else:
    return (((MODMUL[n]*MODDIV[r])%MOD)*MODDIV[n-r])%MOD n<r or n<1 n<r n r n<1 n 1 return 0 0 return (((MODMUL[n]*MODDIV[r])%MOD)*MODDIV[n-r])%MOD (((MODMUL[n]*MODDIV[r])%MOD)*MODDIV[n-r])%MOD ((MODMUL[n]*MODDIV[r])%MOD)*MODDIV[n-r] (MODMUL[n]*MODDIV[r])%MOD MODMUL[n]*MODDIV[r] MODMUL[n] MODMUL n MODDIV[r] MODDIV r MOD MODDIV[n-r] MODDIV n-r n r MOD n n r r N, M, K=map(int, input().split()) N map(int, input().split()) map int input().split() ().split () input split M K ans=0 ans 0 k=MODCOMB(N*M-2, K-2) k MODCOMB(N*M-2, K-2) MODCOMB N*M-2 N*M N M 2 K-2 K 2 for sa in range(1, N):
  ans=(ans+sa*(N-sa)*M*M*k)%MOD
     sa range(1, N) range 1 N ans=(ans+sa*(N-sa)*M*M*k)%MOD ans (ans+sa*(N-sa)*M*M*k)%MOD ans+sa*(N-sa)*M*M*k ans sa*(N-sa)*M*M*k sa*(N-sa)*M*M sa*(N-sa)*M sa*(N-sa) sa N-sa N sa M M k MOD for sa in range(1, M):
  ans=(ans+sa*(M-sa)*N*N*k)%MOD
     sa range(1, M) range 1 M ans=(ans+sa*(M-sa)*N*N*k)%MOD ans (ans+sa*(M-sa)*N*N*k)%MOD ans+sa*(M-sa)*N*N*k ans sa*(M-sa)*N*N*k sa*(M-sa)*N*N sa*(M-sa)*N sa*(M-sa) sa M-sa M sa N N k MOD print(ans) print ans K=map(int, input().split()) map(int, input().split()) K ans=(ans+sa*(N-sa)*M*M*k)%MOD (ans+sa*(N-sa)*M*M*k)%MOD ans ans=(ans+sa*(M-sa)*N*N*k)%MOD (ans+sa*(M-sa)*N*N*k)%MOD ans MOD=10**9+7 10**9+7 MOD def MODCOMB(n, r):
  if n<r or n<1:
    return 0
  else:
    return (((MODMUL[n]*MODDIV[r])%MOD)*MODDIV[n-r])%MOD def MODCOMB(n, r):
  if n<r or n<1:
    return 0
  else:
    return (((MODMUL[n]*MODDIV[r])%MOD)*MODDIV[n-r])%MOD MODCOMB N, M, K=map(int, input().split()) map(int, input().split()) N ans=0 0 ans M, K=map(int, input().split()) map(int, input().split()) M MODMUL=[1, 1]+[0]*(UPPERLIMIT-1) [1, 1]+[0]*(UPPERLIMIT-1) MODMUL k=MODCOMB(N*M-2, K-2) MODCOMB(N*M-2, K-2) k UPPERLIMIT=3*10**5 3*10**5 UPPERLIMIT MODDIV=[1]*UPPERLIMIT+[pow(MODMUL[-1], MOD-2, MOD)] [1]*UPPERLIMIT+[pow(MODMUL[-1], MOD-2, MOD)] MODDIV
MOD=10**9+7
UPPERLIMIT=3*10**5
MODMUL=[1, 1]+[0]*(UPPERLIMIT-1)
for i in range(2, UPPERLIMIT+1):
  MODMUL[i]=MODMUL[i-1]*i%MOD
MODDIV=[1]*UPPERLIMIT+[pow(MODMUL[-1], MOD-2, MOD)]
for i in range(UPPERLIMIT, 0, -1):
  MODDIV[i-1]=MODDIV[i]*i%MOD
def MODCOMB(n, r):
  if n<r or n<1:
    return 0
  else:
    return (((MODMUL[n]*MODDIV[r])%MOD)*MODDIV[n-r])%MOD

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

ans=0
k=MODCOMB(N*M-2, K-2)
for sa in range(1, N):
  ans=(ans+sa*(N-sa)*M*M*k)%MOD
    
for sa in range(1, M):
  ans=(ans+sa*(M-sa)*N*N*k)%MOD
    
print(ans)