mod=10**9+7
MAX=2*10**5+100

g1=[1,1]
g2=[1,1]
for i in range(2,MAX+1):
    num_1=g1[-1]*i%mod
    g1.append(num_1)
    g2.append(pow(num_1,mod-2,mod))
    
def cmb(n,r):
    return g1[n]*g2[r]*g2[n-r]%mod


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

ans=0
for t in range(1,M):
    ans=(ans+t*(N**2)*(M-t))%mod
for t in range(1,N):
    ans=(ans+t*(M**2)*(N-t))%mod

print(ans*cmb(N*M-2,K-2)%mod) mod=10**9+7 mod 10**9+7 10**9 10 9 7 MAX=2*10**5+100 MAX 2*10**5+100 2*10**5 2 10**5 10 5 100 g1=[1,1] g1 [1,1] 1 1 g2=[1,1] g2 [1,1] 1 1 for i in range(2,MAX+1):
    num_1=g1[-1]*i%mod
    g1.append(num_1)
    g2.append(pow(num_1,mod-2,mod))
     i range(2,MAX+1) range 2 MAX+1 MAX 1 num_1=g1[-1]*i%mod num_1 g1[-1]*i%mod g1[-1]*i g1[-1] g1 -1 i mod g1.append(num_1) g1.append g1 append num_1 g2.append(pow(num_1,mod-2,mod)) g2.append g2 append pow(num_1,mod-2,mod) pow num_1 mod-2 mod 2 mod def cmb(n,r):
    return g1[n]*g2[r]*g2[n-r]%mod cmb return g1[n]*g2[r]*g2[n-r]%mod g1[n]*g2[r]*g2[n-r]%mod g1[n]*g2[r]*g2[n-r] g1[n]*g2[r] g1[n] g1 n g2[r] g2 r g2[n-r] g2 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 for t in range(1,M):
    ans=(ans+t*(N**2)*(M-t))%mod t range(1,M) range 1 M ans=(ans+t*(N**2)*(M-t))%mod ans (ans+t*(N**2)*(M-t))%mod ans+t*(N**2)*(M-t) ans t*(N**2)*(M-t) t*(N**2) t N**2 N 2 M-t M t mod for t in range(1,N):
    ans=(ans+t*(M**2)*(N-t))%mod t range(1,N) range 1 N ans=(ans+t*(M**2)*(N-t))%mod ans (ans+t*(M**2)*(N-t))%mod ans+t*(M**2)*(N-t) ans t*(M**2)*(N-t) t*(M**2) t M**2 M 2 N-t N t mod print(ans*cmb(N*M-2,K-2)%mod) print ans*cmb(N*M-2,K-2)%mod ans*cmb(N*M-2,K-2) ans cmb(N*M-2,K-2) cmb N*M-2 N*M N M 2 K-2 K 2 mod ans=(ans+t*(M**2)*(N-t))%mod (ans+t*(M**2)*(N-t))%mod ans def cmb(n,r):
    return g1[n]*g2[r]*g2[n-r]%mod def cmb(n,r):
    return g1[n]*g2[r]*g2[n-r]%mod cmb N,M,K=map(int,input().split()) map(int,input().split()) N M,K=map(int,input().split()) map(int,input().split()) M MAX=2*10**5+100 2*10**5+100 MAX K=map(int,input().split()) map(int,input().split()) K g1=[1,1] [1,1] g1 num_1=g1[-1]*i%mod g1[-1]*i%mod num_1 ans=(ans+t*(N**2)*(M-t))%mod (ans+t*(N**2)*(M-t))%mod ans ans=0 0 ans mod=10**9+7 10**9+7 mod g2=[1,1] [1,1] g2
mod=10**9+7
MAX=2*10**5+100

g1=[1,1]
g2=[1,1]
for i in range(2,MAX+1):
    num_1=g1[-1]*i%mod
    g1.append(num_1)
    g2.append(pow(num_1,mod-2,mod))
    
def cmb(n,r):
    return g1[n]*g2[r]*g2[n-r]%mod


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

ans=0
for t in range(1,M):
    ans=(ans+t*(N**2)*(M-t))%mod
for t in range(1,N):
    ans=(ans+t*(M**2)*(N-t))%mod

print(ans*cmb(N*M-2,K-2)%mod)