import sys
import heapq
import math
import fractions
import bisect
import itertools
from collections import Counter
from collections import deque
from operator import itemgetter
def input(): return sys.stdin.readline().strip()
def mp(): return map(int,input().split())
def lmp(): return list(map(int,input().split()))

mod = 10**9 + 7

# x ** a をmodで割った余りを、O(log(a))時間で求める。
def power(x, a):
    if a == 0:
        return 1
    elif a == 1:
        return x
    elif a % 2 == 0:
        return power(x, a//2) **2 % mod
    else:
        return power(x, a//2) **2 * x % mod

# xの逆元を求める。フェルマーの小定理より、 x の逆元は x ^ (mod - 2) に等しい。計算時間はO(log(mod))程度。
# https://qiita.com/Yaruki00/items/fd1fc269ff7fe40d09a6
def modinv(x):
    return power(x, mod-2)

def binomial_coefficients(n, k):
    numera = 1  # 分子
    denomi = 1  # 分母

    for i in range(k):
        numera *= n-i
        numera %= mod
        denomi *= i+1
        denomi %= mod
    return numera * modinv(denomi) % mod

ans=0
n,m,k=mp()
for i in range(n):
    for j in range(m):
        if i!=0 or j!=0:
            if i==0 or j==0:
                ans+=(n-i)*(m-j)*(i+j)%mod
            else:
                ans+=2*(n-i)*(m-j)*(i+j)%mod
            ans%=mod
print(ans*binomial_coefficients(n*m-2, k-2)%mod)     import sys sys import heapq heapq import math math import fractions fractions import bisect bisect import itertools itertools from collections import Counter from collections import deque from operator import itemgetter def input(): return sys.stdin.readline().strip() input def mp(): return map(int,input().split()) mp def lmp(): return list(map(int,input().split())) lmp mod = 10**9 + 7 mod 10**9 + 7 10**9 10 9 7 def power(x, a):
    if a == 0:
        return 1
    elif a == 1:
        return x
    elif a % 2 == 0:
        return power(x, a//2) **2 % mod
    else:
        return power(x, a//2) **2 * x % mod

# xの逆元を求める。フェルマーの小定理より、 x の逆元は x ^ (mod - 2) に等しい。計算時間はO(log(mod))程度。
# https://qiita.com/Yaruki00/items/fd1fc269ff7fe40d09a6 power if a == 0:
        return 1
    elif a == 1:
        return x
    elif a % 2 == 0:
        return power(x, a//2) **2 % mod
    else:
        return power(x, a//2) **2 * x % mod

# xの逆元を求める。フェルマーの小定理より、 x の逆元は x ^ (mod - 2) に等しい。計算時間はO(log(mod))程度。
# https://qiita.com/Yaruki00/items/fd1fc269ff7fe40d09a6 a == 0 a 0 return 1 1 elif a == 1:
        return x
     a == 1 a 1 return x x elif a % 2 == 0:
        return power(x, a//2) **2 % mod
     a % 2 == 0 a % 2 a 2 0 return power(x, a//2) **2 % mod power(x, a//2) **2 % mod power(x, a//2) **2 power(x, a//2) power x a//2 a 2 2 mod return power(x, a//2) **2 * x % mod power(x, a//2) **2 * x % mod power(x, a//2) **2 * x power(x, a//2) **2 power(x, a//2) power x a//2 a 2 2 x mod x x a a def modinv(x):
    return power(x, mod-2) modinv return power(x, mod-2) power(x, mod-2) power x mod-2 mod 2 x x def binomial_coefficients(n, k):
    numera = 1  # 分子
    denomi = 1  # 分母

    for i in range(k):
        numera *= n-i
        numera %= mod
        denomi *= i+1
        denomi %= mod
    return numera * modinv(denomi) % mod binomial_coefficients numera = 1 numera 1 denomi = 1 denomi 1 for i in range(k):
        numera *= n-i
        numera %= mod
        denomi *= i+1
        denomi %= mod
     i range(k) range k numera *= n-i numera n-i n i numera %= mod numera mod denomi *= i+1 denomi i+1 i 1 denomi %= mod denomi mod return numera * modinv(denomi) % mod numera * modinv(denomi) % mod numera * modinv(denomi) numera modinv(denomi) modinv denomi mod n n k k ans=0 ans 0 n,m,k=mp() n mp() mp m k for i in range(n):
    for j in range(m):
        if i!=0 or j!=0:
            if i==0 or j==0:
                ans+=(n-i)*(m-j)*(i+j)%mod
            else:
                ans+=2*(n-i)*(m-j)*(i+j)%mod
            ans%=mod i range(n) range n for j in range(m):
        if i!=0 or j!=0:
            if i==0 or j==0:
                ans+=(n-i)*(m-j)*(i+j)%mod
            else:
                ans+=2*(n-i)*(m-j)*(i+j)%mod
            ans%=mod j range(m) range m if i!=0 or j!=0:
            if i==0 or j==0:
                ans+=(n-i)*(m-j)*(i+j)%mod
            else:
                ans+=2*(n-i)*(m-j)*(i+j)%mod
            ans%=mod i!=0 or j!=0 i!=0 i 0 j!=0 j 0 if i==0 or j==0:
                ans+=(n-i)*(m-j)*(i+j)%mod
            else:
                ans+=2*(n-i)*(m-j)*(i+j)%mod
             i==0 or j==0 i==0 i 0 j==0 j 0 ans+=(n-i)*(m-j)*(i+j)%mod ans (n-i)*(m-j)*(i+j)%mod (n-i)*(m-j)*(i+j) (n-i)*(m-j) n-i n i m-j m j i+j i j mod ans+=2*(n-i)*(m-j)*(i+j)%mod ans 2*(n-i)*(m-j)*(i+j)%mod 2*(n-i)*(m-j)*(i+j) 2*(n-i)*(m-j) 2*(n-i) 2 n-i n i m-j m j i+j i j mod ans%=mod ans mod print(ans*binomial_coefficients(n*m-2, k-2)%mod) print ans*binomial_coefficients(n*m-2, k-2)%mod ans*binomial_coefficients(n*m-2, k-2) ans binomial_coefficients(n*m-2, k-2) binomial_coefficients n*m-2 n*m n m 2 k-2 k 2 mod def lmp(): return list(map(int,input().split())) def lmp(): return list(map(int,input().split())) lmp k=mp() mp() k def mp(): return map(int,input().split()) def mp(): return map(int,input().split()) mp n,m,k=mp() mp() n m,k=mp() mp() m def input(): return sys.stdin.readline().strip() def input(): return sys.stdin.readline().strip() input ans=0 0 ans ans+=2*(n-i)*(m-j)*(i+j)%mod 2*(n-i)*(m-j)*(i+j)%mod ans def modinv(x):
    return power(x, mod-2) def modinv(x):
    return power(x, mod-2) modinv mod = 10**9 + 7 10**9 + 7 mod def power(x, a):
    if a == 0:
        return 1
    elif a == 1:
        return x
    elif a % 2 == 0:
        return power(x, a//2) **2 % mod
    else:
        return power(x, a//2) **2 * x % mod

# xの逆元を求める。フェルマーの小定理より、 x の逆元は x ^ (mod - 2) に等しい。計算時間はO(log(mod))程度。
# https://qiita.com/Yaruki00/items/fd1fc269ff7fe40d09a6 def power(x, a):
    if a == 0:
        return 1
    elif a == 1:
        return x
    elif a % 2 == 0:
        return power(x, a//2) **2 % mod
    else:
        return power(x, a//2) **2 * x % mod

# xの逆元を求める。フェルマーの小定理より、 x の逆元は x ^ (mod - 2) に等しい。計算時間はO(log(mod))程度。
# https://qiita.com/Yaruki00/items/fd1fc269ff7fe40d09a6 power ans+=(n-i)*(m-j)*(i+j)%mod (n-i)*(m-j)*(i+j)%mod ans ans%=mod mod ans def binomial_coefficients(n, k):
    numera = 1  # 分子
    denomi = 1  # 分母

    for i in range(k):
        numera *= n-i
        numera %= mod
        denomi *= i+1
        denomi %= mod
    return numera * modinv(denomi) % mod def binomial_coefficients(n, k):
    numera = 1  # 分子
    denomi = 1  # 分母

    for i in range(k):
        numera *= n-i
        numera %= mod
        denomi *= i+1
        denomi %= mod
    return numera * modinv(denomi) % mod binomial_coefficients
import sys
import heapq
import math
import fractions
import bisect
import itertools
from collections import Counter
from collections import deque
from operator import itemgetter
def input(): return sys.stdin.readline().strip()
def mp(): return map(int,input().split())
def lmp(): return list(map(int,input().split()))

mod = 10**9 + 7

# x ** a をmodで割った余りを、O(log(a))時間で求める。
def power(x, a):
    if a == 0:
        return 1
    elif a == 1:
        return x
    elif a % 2 == 0:
        return power(x, a//2) **2 % mod
    else:
        return power(x, a//2) **2 * x % mod

# xの逆元を求める。フェルマーの小定理より、 x の逆元は x ^ (mod - 2) に等しい。計算時間はO(log(mod))程度。
# https://qiita.com/Yaruki00/items/fd1fc269ff7fe40d09a6
def modinv(x):
    return power(x, mod-2)

def binomial_coefficients(n, k):
    numera = 1  # 分子
    denomi = 1  # 分母

    for i in range(k):
        numera *= n-i
        numera %= mod
        denomi *= i+1
        denomi %= mod
    return numera * modinv(denomi) % mod

ans=0
n,m,k=mp()
for i in range(n):
    for j in range(m):
        if i!=0 or j!=0:
            if i==0 or j==0:
                ans+=(n-i)*(m-j)*(i+j)%mod
            else:
                ans+=2*(n-i)*(m-j)*(i+j)%mod
            ans%=mod
print(ans*binomial_coefficients(n*m-2, k-2)%mod)    