N,M,K=map(int,input().split())
mod=10**9+7
NM=N*M
kaijo=[0]*(NM+1)
gyaku=[0]*(NM+1)
kaijo[0]=1
for i in range(1,NM+1):
    kaijo[i]=kaijo[i-1]*i%mod
gyaku[NM]=pow(kaijo[NM],mod-2,mod)
for i in range(NM,0,-1):
    gyaku[i-1]=gyaku[i]*i%mod
C=kaijo[NM-2]*gyaku[K-2]*gyaku[NM-K]%mod
ans=0
for k in range(1,M):
    x=k*(M-k)*C*N*N%mod
    ans=(ans+x)%mod
for k in range(1,N):
    x=k*(N-k)*C*M*M%mod
    ans=(ans+x)%mod
print(ans) 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 NM=N*M NM N*M N M kaijo=[0]*(NM+1) kaijo [0]*(NM+1) [0] 0 NM+1 NM 1 gyaku=[0]*(NM+1) gyaku [0]*(NM+1) [0] 0 NM+1 NM 1 kaijo[0]=1 kaijo[0] kaijo 0 1 for i in range(1,NM+1):
    kaijo[i]=kaijo[i-1]*i%mod i range(1,NM+1) range 1 NM+1 NM 1 kaijo[i]=kaijo[i-1]*i%mod kaijo[i] kaijo i kaijo[i-1]*i%mod kaijo[i-1]*i kaijo[i-1] kaijo i-1 i 1 i mod gyaku[NM]=pow(kaijo[NM],mod-2,mod) gyaku[NM] gyaku NM pow(kaijo[NM],mod-2,mod) pow kaijo[NM] kaijo NM mod-2 mod 2 mod for i in range(NM,0,-1):
    gyaku[i-1]=gyaku[i]*i%mod i range(NM,0,-1) range NM 0 -1 gyaku[i-1]=gyaku[i]*i%mod gyaku[i-1] gyaku i-1 i 1 gyaku[i]*i%mod gyaku[i]*i gyaku[i] gyaku i i mod C=kaijo[NM-2]*gyaku[K-2]*gyaku[NM-K]%mod C kaijo[NM-2]*gyaku[K-2]*gyaku[NM-K]%mod kaijo[NM-2]*gyaku[K-2]*gyaku[NM-K] kaijo[NM-2]*gyaku[K-2] kaijo[NM-2] kaijo NM-2 NM 2 gyaku[K-2] gyaku K-2 K 2 gyaku[NM-K] gyaku NM-K NM K mod ans=0 ans 0 for k in range(1,M):
    x=k*(M-k)*C*N*N%mod
    ans=(ans+x)%mod k range(1,M) range 1 M x=k*(M-k)*C*N*N%mod x k*(M-k)*C*N*N%mod k*(M-k)*C*N*N k*(M-k)*C*N k*(M-k)*C k*(M-k) k M-k M k C N N mod ans=(ans+x)%mod ans (ans+x)%mod ans+x ans x mod for k in range(1,N):
    x=k*(N-k)*C*M*M%mod
    ans=(ans+x)%mod k range(1,N) range 1 N x=k*(N-k)*C*M*M%mod x k*(N-k)*C*M*M%mod k*(N-k)*C*M*M k*(N-k)*C*M k*(N-k)*C k*(N-k) k N-k N k C M M mod ans=(ans+x)%mod ans (ans+x)%mod ans+x ans x mod print(ans) print ans kaijo=[0]*(NM+1) [0]*(NM+1) kaijo K=map(int,input().split()) map(int,input().split()) K gyaku=[0]*(NM+1) [0]*(NM+1) gyaku N,M,K=map(int,input().split()) map(int,input().split()) N mod=10**9+7 10**9+7 mod ans=(ans+x)%mod (ans+x)%mod ans C=kaijo[NM-2]*gyaku[K-2]*gyaku[NM-K]%mod kaijo[NM-2]*gyaku[K-2]*gyaku[NM-K]%mod C x=k*(N-k)*C*M*M%mod k*(N-k)*C*M*M%mod x x=k*(M-k)*C*N*N%mod k*(M-k)*C*N*N%mod x M,K=map(int,input().split()) map(int,input().split()) M NM=N*M N*M NM ans=(ans+x)%mod (ans+x)%mod ans ans=0 0 ans
N,M,K=map(int,input().split())
mod=10**9+7
NM=N*M
kaijo=[0]*(NM+1)
gyaku=[0]*(NM+1)
kaijo[0]=1
for i in range(1,NM+1):
    kaijo[i]=kaijo[i-1]*i%mod
gyaku[NM]=pow(kaijo[NM],mod-2,mod)
for i in range(NM,0,-1):
    gyaku[i-1]=gyaku[i]*i%mod
C=kaijo[NM-2]*gyaku[K-2]*gyaku[NM-K]%mod
ans=0
for k in range(1,M):
    x=k*(M-k)*C*N*N%mod
    ans=(ans+x)%mod
for k in range(1,N):
    x=k*(N-k)*C*M*M%mod
    ans=(ans+x)%mod
print(ans)