def comb(n, k):
    x, y = 1, 1
    k = min(k, n-k)
    for i in range(k):
        x *= n - i
        x %= MOD
        y *= i + 1
        y %= MOD
    return x*modinv(y) % MOD

def modinv(a):
    b=MOD
    (x, lastx) = (0, 1)
    (y, lasty) = (1, 0)
    while b != 0:
        q = a // b
        (a, b) = (b, a % b)
        (x, lastx) = (lastx - q * x, x)
        (y, lasty) = (lasty - q * y, y)
    return lastx%MOD

MOD=10**9+7
N,M,K=map(int,input().split())


x=0
for i in range(1,N):
   x=x+i*(N-i)
x=(x*M*M)%MOD

y=0
for j in range(1,M):
    y=y+j*(M-j)
y=(y*N*N)%MOD
    
ans=comb(N*M-2,K-2)*(x+y)%MOD
print(ans) def comb(n, k):
    x, y = 1, 1
    k = min(k, n-k)
    for i in range(k):
        x *= n - i
        x %= MOD
        y *= i + 1
        y %= MOD
    return x*modinv(y) % MOD comb x, y = 1, 1 x 1 y 1 k = min(k, n-k) k min(k, n-k) min k n-k n k for i in range(k):
        x *= n - i
        x %= MOD
        y *= i + 1
        y %= MOD
     i range(k) range k x *= n - i x n - i n i x %= MOD x MOD y *= i + 1 y i + 1 i 1 y %= MOD y MOD return x*modinv(y) % MOD x*modinv(y) % MOD x*modinv(y) x modinv(y) modinv y MOD n n k k def modinv(a):
    b=MOD
    (x, lastx) = (0, 1)
    (y, lasty) = (1, 0)
    while b != 0:
        q = a // b
        (a, b) = (b, a % b)
        (x, lastx) = (lastx - q * x, x)
        (y, lasty) = (lasty - q * y, y)
    return lastx%MOD modinv b=MOD b MOD (x, lastx) = (0, 1) (x, lastx) x lastx (0, 1) 0 1 (y, lasty) = (1, 0) (y, lasty) y lasty (1, 0) 1 0 while b != 0:
        q = a // b
        (a, b) = (b, a % b)
        (x, lastx) = (lastx - q * x, x)
        (y, lasty) = (lasty - q * y, y)
     b != 0 b 0 q = a // b q a // b a b (a, b) = (b, a % b) (a, b) a b (b, a % b) b a % b a b (x, lastx) = (lastx - q * x, x) (x, lastx) x lastx (lastx - q * x, x) lastx - q * x lastx q * x q x x (y, lasty) = (lasty - q * y, y) (y, lasty) y lasty (lasty - q * y, y) lasty - q * y lasty q * y q y y return lastx%MOD lastx%MOD lastx MOD a a MOD=10**9+7 MOD 10**9+7 10**9 10 9 7 N,M,K=map(int,input().split()) N map(int,input().split()) map int input().split() ().split () input split M K x=0 x 0 for i in range(1,N):
   x=x+i*(N-i) i range(1,N) range 1 N x=x+i*(N-i) x x+i*(N-i) x i*(N-i) i N-i N i x=(x*M*M)%MOD x (x*M*M)%MOD x*M*M x*M x M M MOD y=0 y 0 for j in range(1,M):
    y=y+j*(M-j) j range(1,M) range 1 M y=y+j*(M-j) y y+j*(M-j) y j*(M-j) j M-j M j y=(y*N*N)%MOD y (y*N*N)%MOD y*N*N y*N y N N MOD ans=comb(N*M-2,K-2)*(x+y)%MOD ans comb(N*M-2,K-2)*(x+y)%MOD comb(N*M-2,K-2)*(x+y) comb(N*M-2,K-2) comb N*M-2 N*M N M 2 K-2 K 2 x+y x y MOD print(ans) print ans def comb(n, k):
    x, y = 1, 1
    k = min(k, n-k)
    for i in range(k):
        x *= n - i
        x %= MOD
        y *= i + 1
        y %= MOD
    return x*modinv(y) % MOD def comb(n, k):
    x, y = 1, 1
    k = min(k, n-k)
    for i in range(k):
        x *= n - i
        x %= MOD
        y *= i + 1
        y %= MOD
    return x*modinv(y) % MOD comb def modinv(a):
    b=MOD
    (x, lastx) = (0, 1)
    (y, lasty) = (1, 0)
    while b != 0:
        q = a // b
        (a, b) = (b, a % b)
        (x, lastx) = (lastx - q * x, x)
        (y, lasty) = (lasty - q * y, y)
    return lastx%MOD def modinv(a):
    b=MOD
    (x, lastx) = (0, 1)
    (y, lasty) = (1, 0)
    while b != 0:
        q = a // b
        (a, b) = (b, a % b)
        (x, lastx) = (lastx - q * x, x)
        (y, lasty) = (lasty - q * y, y)
    return lastx%MOD modinv N,M,K=map(int,input().split()) map(int,input().split()) N y=0 0 y x=x+i*(N-i) x+i*(N-i) x M,K=map(int,input().split()) map(int,input().split()) M ans=comb(N*M-2,K-2)*(x+y)%MOD comb(N*M-2,K-2)*(x+y)%MOD ans y=(y*N*N)%MOD (y*N*N)%MOD y x=(x*M*M)%MOD (x*M*M)%MOD x K=map(int,input().split()) map(int,input().split()) K MOD=10**9+7 10**9+7 MOD x=0 0 x y=y+j*(M-j) y+j*(M-j) y
def comb(n, k):
    x, y = 1, 1
    k = min(k, n-k)
    for i in range(k):
        x *= n - i
        x %= MOD
        y *= i + 1
        y %= MOD
    return x*modinv(y) % MOD

def modinv(a):
    b=MOD
    (x, lastx) = (0, 1)
    (y, lasty) = (1, 0)
    while b != 0:
        q = a // b
        (a, b) = (b, a % b)
        (x, lastx) = (lastx - q * x, x)
        (y, lasty) = (lasty - q * y, y)
    return lastx%MOD

MOD=10**9+7
N,M,K=map(int,input().split())


x=0
for i in range(1,N):
   x=x+i*(N-i)
x=(x*M*M)%MOD

y=0
for j in range(1,M):
    y=y+j*(M-j)
y=(y*N*N)%MOD
    
ans=comb(N*M-2,K-2)*(x+y)%MOD
print(ans)