N,M,K=map(int,input().split())
mod=10**9+7
X=N*M%mod
C=[1]
for i in range(K+1):
  C.append(C[i]*(X-2-i)*pow(i+1,mod-2,mod)%mod)
P=0
for i in range(N):
  P=(P+(i*(i+1)//2)+((N-i)*(N-i-1)//2))%mod
P=P*M*M%mod
Q=0
for i in range(M):
  Q=(Q+(i*(i+1)//2)+((M-i)*(M-i-1)//2))%mod
Q=Q*N*N%mod
print((P+Q)*C[K-2]*((mod+1)//2)%mod) 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 X=N*M%mod X N*M%mod N*M N M mod C=[1] C [1] 1 for i in range(K+1):
  C.append(C[i]*(X-2-i)*pow(i+1,mod-2,mod)%mod) i range(K+1) range K+1 K 1 C.append(C[i]*(X-2-i)*pow(i+1,mod-2,mod)%mod) C.append C append C[i]*(X-2-i)*pow(i+1,mod-2,mod)%mod C[i]*(X-2-i)*pow(i+1,mod-2,mod) C[i]*(X-2-i) C[i] C i X-2-i X-2 X 2 i pow(i+1,mod-2,mod) pow i+1 i 1 mod-2 mod 2 mod mod P=0 P 0 for i in range(N):
  P=(P+(i*(i+1)//2)+((N-i)*(N-i-1)//2))%mod i range(N) range N P=(P+(i*(i+1)//2)+((N-i)*(N-i-1)//2))%mod P (P+(i*(i+1)//2)+((N-i)*(N-i-1)//2))%mod P+(i*(i+1)//2)+((N-i)*(N-i-1)//2) P+(i*(i+1)//2) P i*(i+1)//2 i*(i+1) i i+1 i 1 2 (N-i)*(N-i-1)//2 (N-i)*(N-i-1) N-i N i N-i-1 N-i N i 1 2 mod P=P*M*M%mod P P*M*M%mod P*M*M P*M P M M mod Q=0 Q 0 for i in range(M):
  Q=(Q+(i*(i+1)//2)+((M-i)*(M-i-1)//2))%mod i range(M) range M Q=(Q+(i*(i+1)//2)+((M-i)*(M-i-1)//2))%mod Q (Q+(i*(i+1)//2)+((M-i)*(M-i-1)//2))%mod Q+(i*(i+1)//2)+((M-i)*(M-i-1)//2) Q+(i*(i+1)//2) Q i*(i+1)//2 i*(i+1) i i+1 i 1 2 (M-i)*(M-i-1)//2 (M-i)*(M-i-1) M-i M i M-i-1 M-i M i 1 2 mod Q=Q*N*N%mod Q Q*N*N%mod Q*N*N Q*N Q N N mod print((P+Q)*C[K-2]*((mod+1)//2)%mod) print (P+Q)*C[K-2]*((mod+1)//2)%mod (P+Q)*C[K-2]*((mod+1)//2) (P+Q)*C[K-2] P+Q P Q C[K-2] C K-2 K 2 (mod+1)//2 mod+1 mod 1 2 mod P=P*M*M%mod P*M*M%mod P C=[1] [1] C P=(P+(i*(i+1)//2)+((N-i)*(N-i-1)//2))%mod (P+(i*(i+1)//2)+((N-i)*(N-i-1)//2))%mod P Q=(Q+(i*(i+1)//2)+((M-i)*(M-i-1)//2))%mod (Q+(i*(i+1)//2)+((M-i)*(M-i-1)//2))%mod Q K=map(int,input().split()) map(int,input().split()) K M,K=map(int,input().split()) map(int,input().split()) M P=0 0 P X=N*M%mod N*M%mod X Q=0 0 Q mod=10**9+7 10**9+7 mod N,M,K=map(int,input().split()) map(int,input().split()) N Q=Q*N*N%mod Q*N*N%mod Q
N,M,K=map(int,input().split())
mod=10**9+7
X=N*M%mod
C=[1]
for i in range(K+1):
  C.append(C[i]*(X-2-i)*pow(i+1,mod-2,mod)%mod)
P=0
for i in range(N):
  P=(P+(i*(i+1)//2)+((N-i)*(N-i-1)//2))%mod
P=P*M*M%mod
Q=0
for i in range(M):
  Q=(Q+(i*(i+1)//2)+((M-i)*(M-i-1)//2))%mod
Q=Q*N*N%mod
print((P+Q)*C[K-2]*((mod+1)//2)%mod)