MOD=10**9+7

def com(n,k,mod):
    res=1
    tmp=1
    for i in range(1,k+1):
        res=res*(n-i+1)%mod
        tmp=tmp*i%mod
    a=pow(tmp,mod-2,mod)
    return res*a%mod

N,M,K=map(int,input().split())
C=com(N*M-2,K-2,MOD)

X,Y=0,0
for d in range(1,N):
    X=(X+((N-d)*M**2*d))%MOD
for d in range(1,M):
    Y=(Y+((M-d)*N**2*d))%MOD
print((X+Y)*C%MOD) MOD=10**9+7 MOD 10**9+7 10**9 10 9 7 def com(n,k,mod):
    res=1
    tmp=1
    for i in range(1,k+1):
        res=res*(n-i+1)%mod
        tmp=tmp*i%mod
    a=pow(tmp,mod-2,mod)
    return res*a%mod com res=1 res 1 tmp=1 tmp 1 for i in range(1,k+1):
        res=res*(n-i+1)%mod
        tmp=tmp*i%mod
     i range(1,k+1) range 1 k+1 k 1 res=res*(n-i+1)%mod res res*(n-i+1)%mod res*(n-i+1) res n-i+1 n-i n i 1 mod tmp=tmp*i%mod tmp tmp*i%mod tmp*i tmp i mod a=pow(tmp,mod-2,mod) a pow(tmp,mod-2,mod) pow tmp mod-2 mod 2 mod return res*a%mod res*a%mod res*a res a mod n n k k mod mod N,M,K=map(int,input().split()) N map(int,input().split()) map int input().split() ().split () input split M K C=com(N*M-2,K-2,MOD) C com(N*M-2,K-2,MOD) com N*M-2 N*M N M 2 K-2 K 2 MOD X,Y=0,0 X 0 Y 0 for d in range(1,N):
    X=(X+((N-d)*M**2*d))%MOD d range(1,N) range 1 N X=(X+((N-d)*M**2*d))%MOD X (X+((N-d)*M**2*d))%MOD X+((N-d)*M**2*d) X (N-d)*M**2*d (N-d)*M**2 N-d N d M**2 M 2 d MOD for d in range(1,M):
    Y=(Y+((M-d)*N**2*d))%MOD d range(1,M) range 1 M Y=(Y+((M-d)*N**2*d))%MOD Y (Y+((M-d)*N**2*d))%MOD Y+((M-d)*N**2*d) Y (M-d)*N**2*d (M-d)*N**2 M-d M d N**2 N 2 d MOD print((X+Y)*C%MOD) print (X+Y)*C%MOD (X+Y)*C X+Y X Y C MOD X=(X+((N-d)*M**2*d))%MOD (X+((N-d)*M**2*d))%MOD X N,M,K=map(int,input().split()) map(int,input().split()) N Y=0,0 0 Y M,K=map(int,input().split()) map(int,input().split()) M K=map(int,input().split()) map(int,input().split()) K Y=(Y+((M-d)*N**2*d))%MOD (Y+((M-d)*N**2*d))%MOD Y X,Y=0,0 0 X MOD=10**9+7 10**9+7 MOD C=com(N*M-2,K-2,MOD) com(N*M-2,K-2,MOD) C def com(n,k,mod):
    res=1
    tmp=1
    for i in range(1,k+1):
        res=res*(n-i+1)%mod
        tmp=tmp*i%mod
    a=pow(tmp,mod-2,mod)
    return res*a%mod def com(n,k,mod):
    res=1
    tmp=1
    for i in range(1,k+1):
        res=res*(n-i+1)%mod
        tmp=tmp*i%mod
    a=pow(tmp,mod-2,mod)
    return res*a%mod com
MOD=10**9+7

def com(n,k,mod):
    res=1
    tmp=1
    for i in range(1,k+1):
        res=res*(n-i+1)%mod
        tmp=tmp*i%mod
    a=pow(tmp,mod-2,mod)
    return res*a%mod

N,M,K=map(int,input().split())
C=com(N*M-2,K-2,MOD)

X,Y=0,0
for d in range(1,N):
    X=(X+((N-d)*M**2*d))%MOD
for d in range(1,M):
    Y=(Y+((M-d)*N**2*d))%MOD
print((X+Y)*C%MOD)