n,m,k=map(int,input().split())
num=n*m
MOD = 10**9+7

def inv_mod(a, p=MOD):
    def inv_mod_sub(a, p):
        if a == 1:
            return 1, 0
        else:
            d, r = p//a, p%a
            x, y = inv_mod_sub(r, a)
            return y-d*x, x
    if p < 0: p = -p
    a %= p
    return inv_mod_sub(a,p)[0] % p

def comb_mod(n, k):
    if k < 0 or k > n:
        return 0
    else:
        return f_mod[n]*f_mod_inverse[k]*f_mod_inverse[n-k] % MOD
f_mod=[1]*(num+1)
f_mod_inverse=[1]*(num+1)
for i in range(1,num+1):
  f_mod[i]=f_mod[i-1]*i%MOD
  f_mod_inverse[i]=f_mod_inverse[i-1]*inv_mod(i)%MOD
N=0;M=0
for i in range(1,n):
  N+=i*(n-i)
  N%=MOD
for i in range(1,m):
  M+=i*(m-i)
  M%=MOD
print((N*(m**2)%MOD+M*(n**2)%MOD)*comb_mod(n*m-2,k-2)%MOD) n,m,k=map(int,input().split()) n map(int,input().split()) map int input().split() ().split () input split m k num=n*m num n*m n m MOD = 10**9+7 MOD 10**9+7 10**9 10 9 7 def inv_mod(a, p=MOD):
    def inv_mod_sub(a, p):
        if a == 1:
            return 1, 0
        else:
            d, r = p//a, p%a
            x, y = inv_mod_sub(r, a)
            return y-d*x, x
    if p < 0: p = -p
    a %= p
    return inv_mod_sub(a,p)[0] % p inv_mod def inv_mod_sub(a, p):
        if a == 1:
            return 1, 0
        else:
            d, r = p//a, p%a
            x, y = inv_mod_sub(r, a)
            return y-d*x, x
     inv_mod_sub if a == 1:
            return 1, 0
        else:
            d, r = p//a, p%a
            x, y = inv_mod_sub(r, a)
            return y-d*x, x
     a == 1 a 1 return 1, 0 return 1, 0 1 0 d, r = p//a, p%a d p//a p a r p%a p a x, y = inv_mod_sub(r, a) x inv_mod_sub(r, a) inv_mod_sub r a y return y-d*x, x return y-d*x, x y-d*x y d*x d x x a a p p if p < 0: p = -p
     p < 0 p 0 a %= p a p return inv_mod_sub(a,p)[0] % p inv_mod_sub(a,p)[0] % p inv_mod_sub(a,p)[0] (a,p) inv_mod_sub a p 0 p a a p=MOD p MOD def comb_mod(n, k):
    if k < 0 or k > n:
        return 0
    else:
        return f_mod[n]*f_mod_inverse[k]*f_mod_inverse[n-k] % MOD comb_mod if k < 0 or k > n:
        return 0
    else:
        return f_mod[n]*f_mod_inverse[k]*f_mod_inverse[n-k] % MOD k < 0 or k > n k < 0 k 0 k > n k n return 0 0 return f_mod[n]*f_mod_inverse[k]*f_mod_inverse[n-k] % MOD f_mod[n]*f_mod_inverse[k]*f_mod_inverse[n-k] % MOD f_mod[n]*f_mod_inverse[k]*f_mod_inverse[n-k] f_mod[n]*f_mod_inverse[k] f_mod[n] f_mod n f_mod_inverse[k] f_mod_inverse k f_mod_inverse[n-k] f_mod_inverse n-k n k MOD n n k k f_mod=[1]*(num+1) f_mod [1]*(num+1) [1] 1 num+1 num 1 f_mod_inverse=[1]*(num+1) f_mod_inverse [1]*(num+1) [1] 1 num+1 num 1 for i in range(1,num+1):
  f_mod[i]=f_mod[i-1]*i%MOD
  f_mod_inverse[i]=f_mod_inverse[i-1]*inv_mod(i)%MOD i range(1,num+1) range 1 num+1 num 1 f_mod[i]=f_mod[i-1]*i%MOD f_mod[i] f_mod i f_mod[i-1]*i%MOD f_mod[i-1]*i f_mod[i-1] f_mod i-1 i 1 i MOD f_mod_inverse[i]=f_mod_inverse[i-1]*inv_mod(i)%MOD f_mod_inverse[i] f_mod_inverse i f_mod_inverse[i-1]*inv_mod(i)%MOD f_mod_inverse[i-1]*inv_mod(i) f_mod_inverse[i-1] f_mod_inverse i-1 i 1 inv_mod(i) inv_mod i MOD N=0 N 0 M=0 M 0 for i in range(1,n):
  N+=i*(n-i)
  N%=MOD i range(1,n) range 1 n N+=i*(n-i) N i*(n-i) i n-i n i N%=MOD N MOD for i in range(1,m):
  M+=i*(m-i)
  M%=MOD i range(1,m) range 1 m M+=i*(m-i) M i*(m-i) i m-i m i M%=MOD M MOD print((N*(m**2)%MOD+M*(n**2)%MOD)*comb_mod(n*m-2,k-2)%MOD) print (N*(m**2)%MOD+M*(n**2)%MOD)*comb_mod(n*m-2,k-2)%MOD (N*(m**2)%MOD+M*(n**2)%MOD)*comb_mod(n*m-2,k-2) N*(m**2)%MOD+M*(n**2)%MOD N*(m**2)%MOD N*(m**2) N m**2 m 2 MOD M*(n**2)%MOD M*(n**2) M n**2 n 2 MOD comb_mod(n*m-2,k-2) comb_mod n*m-2 n*m n m 2 k-2 k 2 MOD M+=i*(m-i) i*(m-i) M def inv_mod(a, p=MOD):
    def inv_mod_sub(a, p):
        if a == 1:
            return 1, 0
        else:
            d, r = p//a, p%a
            x, y = inv_mod_sub(r, a)
            return y-d*x, x
    if p < 0: p = -p
    a %= p
    return inv_mod_sub(a,p)[0] % p def inv_mod(a, p=MOD):
    def inv_mod_sub(a, p):
        if a == 1:
            return 1, 0
        else:
            d, r = p//a, p%a
            x, y = inv_mod_sub(r, a)
            return y-d*x, x
    if p < 0: p = -p
    a %= p
    return inv_mod_sub(a,p)[0] % p inv_mod MOD = 10**9+7 10**9+7 MOD N+=i*(n-i) i*(n-i) N k=map(int,input().split()) map(int,input().split()) k f_mod_inverse=[1]*(num+1) [1]*(num+1) f_mod_inverse M%=MOD MOD M n,m,k=map(int,input().split()) map(int,input().split()) n M=0 0 M m,k=map(int,input().split()) map(int,input().split()) m num=n*m n*m num def comb_mod(n, k):
    if k < 0 or k > n:
        return 0
    else:
        return f_mod[n]*f_mod_inverse[k]*f_mod_inverse[n-k] % MOD def comb_mod(n, k):
    if k < 0 or k > n:
        return 0
    else:
        return f_mod[n]*f_mod_inverse[k]*f_mod_inverse[n-k] % MOD comb_mod N=0 0 N N%=MOD MOD N f_mod=[1]*(num+1) [1]*(num+1) f_mod
n,m,k=map(int,input().split())
num=n*m
MOD = 10**9+7

def inv_mod(a, p=MOD):
    def inv_mod_sub(a, p):
        if a == 1:
            return 1, 0
        else:
            d, r = p//a, p%a
            x, y = inv_mod_sub(r, a)
            return y-d*x, x
    if p < 0: p = -p
    a %= p
    return inv_mod_sub(a,p)[0] % p

def comb_mod(n, k):
    if k < 0 or k > n:
        return 0
    else:
        return f_mod[n]*f_mod_inverse[k]*f_mod_inverse[n-k] % MOD
f_mod=[1]*(num+1)
f_mod_inverse=[1]*(num+1)
for i in range(1,num+1):
  f_mod[i]=f_mod[i-1]*i%MOD
  f_mod_inverse[i]=f_mod_inverse[i-1]*inv_mod(i)%MOD
N=0;M=0
for i in range(1,n):
  N+=i*(n-i)
  N%=MOD
for i in range(1,m):
  M+=i*(m-i)
  M%=MOD
print((N*(m**2)%MOD+M*(n**2)%MOD)*comb_mod(n*m-2,k-2)%MOD)