import sys
input = sys.stdin.readline

def I(): return int(input())
def MI(): return map(int, input().split())
def LI(): return list(map(int, input().split()))

def main():
    mod=10**9+7
    N,M,K=MI()
    #2この組み合わせを全部足して，平均ぽく考える
    
    def sum2(x):
        return ((x*(x+1))//2)%mod
    
    #i,jますにいるときの相手．（i,j）よりも奥にあるやつの全組み合わせの総和
    def calc(i,j):
        temp=0
        temp=sum2(M-j-1)#該当する行
        temp=(temp+((sum2(M-j-1)+sum2(j))*(N-i-1)))%mod#下の行の階差部分，左右方向にあるから
        temp=(temp+sum2(N-i-1)*M)%mod
        return temp
        
    S=0
    for i in range(N):
        for j in range(M):
            S=(S+calc(i,j))%mod
            
    rem=(N*M)-2
    
    def cmb(n, r, mod):
        if (r < 0) or (n < r):
            return 0
        r = min(r, n - r)
        return (fact[n] * factinv[r] * factinv[n-r])%mod

    fact=[1,1]
    factinv=[1,1]
    inv=[0,1]
    
    for i in range(2, rem + 1):
        fact.append((fact[-1] * i) % mod)
        inv.append((-inv[mod % i] * (mod // i)) % mod)
        factinv.append((factinv[-1] * inv[-1]) % mod)  
        
    #print(rem,K-2)
    S=(S*cmb(rem,K-2,mod))%mod
            
    print(S)

main() import sys sys input = sys.stdin.readline input sys.stdin.readline sys.stdin sys stdin readline def I(): return int(input()) I def MI(): return map(int, input().split()) MI def LI(): return list(map(int, input().split())) LI def main():
    mod=10**9+7
    N,M,K=MI()
    #2この組み合わせを全部足して，平均ぽく考える
    
    def sum2(x):
        return ((x*(x+1))//2)%mod
    
    #i,jますにいるときの相手．（i,j）よりも奥にあるやつの全組み合わせの総和
    def calc(i,j):
        temp=0
        temp=sum2(M-j-1)#該当する行
        temp=(temp+((sum2(M-j-1)+sum2(j))*(N-i-1)))%mod#下の行の階差部分，左右方向にあるから
        temp=(temp+sum2(N-i-1)*M)%mod
        return temp
        
    S=0
    for i in range(N):
        for j in range(M):
            S=(S+calc(i,j))%mod
            
    rem=(N*M)-2
    
    def cmb(n, r, mod):
        if (r < 0) or (n < r):
            return 0
        r = min(r, n - r)
        return (fact[n] * factinv[r] * factinv[n-r])%mod

    fact=[1,1]
    factinv=[1,1]
    inv=[0,1]
    
    for i in range(2, rem + 1):
        fact.append((fact[-1] * i) % mod)
        inv.append((-inv[mod % i] * (mod // i)) % mod)
        factinv.append((factinv[-1] * inv[-1]) % mod)  
        
    #print(rem,K-2)
    S=(S*cmb(rem,K-2,mod))%mod
            
    print(S) main mod=10**9+7 mod 10**9+7 10**9 10 9 7 N,M,K=MI() N MI() MI M K def sum2(x):
        return ((x*(x+1))//2)%mod
    
    #i,jますにいるときの相手．（i,j）よりも奥にあるやつの全組み合わせの総和
     sum2 return ((x*(x+1))//2)%mod ((x*(x+1))//2)%mod (x*(x+1))//2 x*(x+1) x x+1 x 1 2 mod x x def calc(i,j):
        temp=0
        temp=sum2(M-j-1)#該当する行
        temp=(temp+((sum2(M-j-1)+sum2(j))*(N-i-1)))%mod#下の行の階差部分，左右方向にあるから
        temp=(temp+sum2(N-i-1)*M)%mod
        return temp
        
     calc temp=0 temp 0 temp=sum2(M-j-1) temp sum2(M-j-1) sum2 M-j-1 M-j M j 1 temp=(temp+((sum2(M-j-1)+sum2(j))*(N-i-1)))%mod temp (temp+((sum2(M-j-1)+sum2(j))*(N-i-1)))%mod temp+((sum2(M-j-1)+sum2(j))*(N-i-1)) temp (sum2(M-j-1)+sum2(j))*(N-i-1) sum2(M-j-1)+sum2(j) sum2(M-j-1) sum2 M-j-1 M-j M j 1 sum2(j) sum2 j N-i-1 N-i N i 1 mod temp=(temp+sum2(N-i-1)*M)%mod temp (temp+sum2(N-i-1)*M)%mod temp+sum2(N-i-1)*M temp sum2(N-i-1)*M sum2(N-i-1) sum2 N-i-1 N-i N i 1 M mod return temp temp i i j j S=0 S 0 for i in range(N):
        for j in range(M):
            S=(S+calc(i,j))%mod
            
     i range(N) range N for j in range(M):
            S=(S+calc(i,j))%mod
            
     j range(M) range M S=(S+calc(i,j))%mod S (S+calc(i,j))%mod S+calc(i,j) S calc(i,j) calc i j mod rem=(N*M)-2 rem (N*M)-2 N*M N M 2 def cmb(n, r, mod):
        if (r < 0) or (n < r):
            return 0
        r = min(r, n - r)
        return (fact[n] * factinv[r] * factinv[n-r])%mod

     cmb if (r < 0) or (n < r):
            return 0
         (r < 0) or (n < r) r < 0 r 0 n < r n r return 0 0 r = min(r, n - r) r min(r, n - r) min r n - r n r return (fact[n] * factinv[r] * factinv[n-r])%mod (fact[n] * factinv[r] * factinv[n-r])%mod fact[n] * factinv[r] * factinv[n-r] fact[n] * factinv[r] fact[n] fact n factinv[r] factinv r factinv[n-r] factinv n-r n r mod n n r r mod mod fact=[1,1] fact [1,1] 1 1 factinv=[1,1] factinv [1,1] 1 1 inv=[0,1] inv [0,1] 0 1 for i in range(2, rem + 1):
        fact.append((fact[-1] * i) % mod)
        inv.append((-inv[mod % i] * (mod // i)) % mod)
        factinv.append((factinv[-1] * inv[-1]) % mod)  
        
    #print(rem,K-2)
     i range(2, rem + 1) range 2 rem + 1 rem 1 fact.append((fact[-1] * i) % mod) fact.append fact append (fact[-1] * i) % mod fact[-1] * i fact[-1] fact -1 i mod inv.append((-inv[mod % i] * (mod // i)) % mod) inv.append inv append (-inv[mod % i] * (mod // i)) % mod -inv[mod % i] * (mod // i) -inv[mod % i] inv[mod % i] inv mod % i mod i mod // i mod i mod factinv.append((factinv[-1] * inv[-1]) % mod) factinv.append factinv append (factinv[-1] * inv[-1]) % mod factinv[-1] * inv[-1] factinv[-1] factinv -1 inv[-1] inv -1 mod S=(S*cmb(rem,K-2,mod))%mod S (S*cmb(rem,K-2,mod))%mod S*cmb(rem,K-2,mod) S cmb(rem,K-2,mod) cmb rem K-2 K 2 mod mod print(S) print S main() main def LI(): return list(map(int, input().split())) def LI(): return list(map(int, input().split())) LI def main():
    mod=10**9+7
    N,M,K=MI()
    #2この組み合わせを全部足して，平均ぽく考える
    
    def sum2(x):
        return ((x*(x+1))//2)%mod
    
    #i,jますにいるときの相手．（i,j）よりも奥にあるやつの全組み合わせの総和
    def calc(i,j):
        temp=0
        temp=sum2(M-j-1)#該当する行
        temp=(temp+((sum2(M-j-1)+sum2(j))*(N-i-1)))%mod#下の行の階差部分，左右方向にあるから
        temp=(temp+sum2(N-i-1)*M)%mod
        return temp
        
    S=0
    for i in range(N):
        for j in range(M):
            S=(S+calc(i,j))%mod
            
    rem=(N*M)-2
    
    def cmb(n, r, mod):
        if (r < 0) or (n < r):
            return 0
        r = min(r, n - r)
        return (fact[n] * factinv[r] * factinv[n-r])%mod

    fact=[1,1]
    factinv=[1,1]
    inv=[0,1]
    
    for i in range(2, rem + 1):
        fact.append((fact[-1] * i) % mod)
        inv.append((-inv[mod % i] * (mod // i)) % mod)
        factinv.append((factinv[-1] * inv[-1]) % mod)  
        
    #print(rem,K-2)
    S=(S*cmb(rem,K-2,mod))%mod
            
    print(S) def main():
    mod=10**9+7
    N,M,K=MI()
    #2この組み合わせを全部足して，平均ぽく考える
    
    def sum2(x):
        return ((x*(x+1))//2)%mod
    
    #i,jますにいるときの相手．（i,j）よりも奥にあるやつの全組み合わせの総和
    def calc(i,j):
        temp=0
        temp=sum2(M-j-1)#該当する行
        temp=(temp+((sum2(M-j-1)+sum2(j))*(N-i-1)))%mod#下の行の階差部分，左右方向にあるから
        temp=(temp+sum2(N-i-1)*M)%mod
        return temp
        
    S=0
    for i in range(N):
        for j in range(M):
            S=(S+calc(i,j))%mod
            
    rem=(N*M)-2
    
    def cmb(n, r, mod):
        if (r < 0) or (n < r):
            return 0
        r = min(r, n - r)
        return (fact[n] * factinv[r] * factinv[n-r])%mod

    fact=[1,1]
    factinv=[1,1]
    inv=[0,1]
    
    for i in range(2, rem + 1):
        fact.append((fact[-1] * i) % mod)
        inv.append((-inv[mod % i] * (mod // i)) % mod)
        factinv.append((factinv[-1] * inv[-1]) % mod)  
        
    #print(rem,K-2)
    S=(S*cmb(rem,K-2,mod))%mod
            
    print(S) main def I(): return int(input()) def I(): return int(input()) I def MI(): return map(int, input().split()) def MI(): return map(int, input().split()) MI input = sys.stdin.readline sys.stdin.readline input
import sys
input = sys.stdin.readline

def I(): return int(input())
def MI(): return map(int, input().split())
def LI(): return list(map(int, input().split()))

def main():
    mod=10**9+7
    N,M,K=MI()
    #2この組み合わせを全部足して，平均ぽく考える
    
    def sum2(x):
        return ((x*(x+1))//2)%mod
    
    #i,jますにいるときの相手．（i,j）よりも奥にあるやつの全組み合わせの総和
    def calc(i,j):
        temp=0
        temp=sum2(M-j-1)#該当する行
        temp=(temp+((sum2(M-j-1)+sum2(j))*(N-i-1)))%mod#下の行の階差部分，左右方向にあるから
        temp=(temp+sum2(N-i-1)*M)%mod
        return temp
        
    S=0
    for i in range(N):
        for j in range(M):
            S=(S+calc(i,j))%mod
            
    rem=(N*M)-2
    
    def cmb(n, r, mod):
        if (r < 0) or (n < r):
            return 0
        r = min(r, n - r)
        return (fact[n] * factinv[r] * factinv[n-r])%mod

    fact=[1,1]
    factinv=[1,1]
    inv=[0,1]
    
    for i in range(2, rem + 1):
        fact.append((fact[-1] * i) % mod)
        inv.append((-inv[mod % i] * (mod // i)) % mod)
        factinv.append((factinv[-1] * inv[-1]) % mod)  
        
    #print(rem,K-2)
    S=(S*cmb(rem,K-2,mod))%mod
            
    print(S)

main()
