def egcd(a, b):
    (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, lasty, a)

# ax ≡ 1 (mod m)なるxを返す
def modinv(a):
    m=MOD
    (inv, q, gcd_val) = egcd(a, m)
    return inv % m

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

b=1
for i in range(N*M-2,N*M-K,-1):
    b=(b*i)%MOD

for i in range(1,K-1):
    b=(b*modinv(i))%MOD

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=b*(x+y)%MOD
print(ans) def egcd(a, b):
    (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, lasty, a)

# ax ≡ 1 (mod m)なるxを返す egcd (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, lasty, a) (lastx, lasty, a) lastx lasty a a a b b def modinv(a):
    m=MOD
    (inv, q, gcd_val) = egcd(a, m)
    return inv % m modinv m=MOD m MOD (inv, q, gcd_val) = egcd(a, m) (inv, q, gcd_val) inv q gcd_val egcd(a, m) egcd a m return inv % m inv % m inv m 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 b=1 b 1 for i in range(N*M-2,N*M-K,-1):
    b=(b*i)%MOD i range(N*M-2,N*M-K,-1) range N*M-2 N*M N M 2 N*M-K N*M N M K -1 b=(b*i)%MOD b (b*i)%MOD b*i b i MOD for i in range(1,K-1):
    b=(b*modinv(i))%MOD i range(1,K-1) range 1 K-1 K 1 b=(b*modinv(i))%MOD b (b*modinv(i))%MOD b*modinv(i) b modinv(i) modinv i MOD 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=b*(x+y)%MOD ans b*(x+y)%MOD b*(x+y) b x+y x y MOD print(ans) print ans K=map(int,input().split()) map(int,input().split()) K N,M,K=map(int,input().split()) map(int,input().split()) N y=y+j*(M-j) y+j*(M-j) y y=0 0 y x=0 0 x b=(b*i)%MOD (b*i)%MOD b y=(y*N*N)%MOD (y*N*N)%MOD y x=(x*M*M)%MOD (x*M*M)%MOD x M,K=map(int,input().split()) map(int,input().split()) M def modinv(a):
    m=MOD
    (inv, q, gcd_val) = egcd(a, m)
    return inv % m def modinv(a):
    m=MOD
    (inv, q, gcd_val) = egcd(a, m)
    return inv % m modinv MOD=10**9+7 10**9+7 MOD b=1 1 b def egcd(a, b):
    (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, lasty, a)

# ax ≡ 1 (mod m)なるxを返す def egcd(a, b):
    (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, lasty, a)

# ax ≡ 1 (mod m)なるxを返す egcd x=x+i*(N-i) x+i*(N-i) x ans=b*(x+y)%MOD b*(x+y)%MOD ans b=(b*modinv(i))%MOD (b*modinv(i))%MOD b
def egcd(a, b):
    (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, lasty, a)

# ax ≡ 1 (mod m)なるxを返す
def modinv(a):
    m=MOD
    (inv, q, gcd_val) = egcd(a, m)
    return inv % m

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

b=1
for i in range(N*M-2,N*M-K,-1):
    b=(b*i)%MOD

for i in range(1,K-1):
    b=(b*modinv(i))%MOD

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=b*(x+y)%MOD
print(ans)