n,m,k=map(int,input().split())

mod = 10**9+7
f = [1]
for i in range(n*m):
  f += [f[-1]*(i+1) % mod]
def comb(a, b):
  return f[a]*pow(f[b], mod-2, mod)*pow(f[a-b], mod-2, mod) % mod

ans=0
for dist in range(1,n):
  c=n-dist
  p=(c*m*m)%mod
  ans+=p*dist
for dist in range(1,m):
  c=m-dist
  p=(c*n*n)%mod
  ans+=p*dist
ans%=mod
print(ans*comb(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 mod = 10**9+7 mod 10**9+7 10**9 10 9 7 f = [1] f [1] 1 for i in range(n*m):
  f += [f[-1]*(i+1) % mod] i range(n*m) range n*m n m f += [f[-1]*(i+1) % mod] f [f[-1]*(i+1) % mod] f[-1]*(i+1) % mod f[-1]*(i+1) f[-1] f -1 i+1 i 1 mod def comb(a, b):
  return f[a]*pow(f[b], mod-2, mod)*pow(f[a-b], mod-2, mod) % mod comb return f[a]*pow(f[b], mod-2, mod)*pow(f[a-b], mod-2, mod) % mod f[a]*pow(f[b], mod-2, mod)*pow(f[a-b], mod-2, mod) % mod f[a]*pow(f[b], mod-2, mod)*pow(f[a-b], mod-2, mod) f[a]*pow(f[b], mod-2, mod) f[a] f a pow(f[b], mod-2, mod) pow f[b] f b mod-2 mod 2 mod pow(f[a-b], mod-2, mod) pow f[a-b] f a-b a b mod-2 mod 2 mod mod a a b b ans=0 ans 0 for dist in range(1,n):
  c=n-dist
  p=(c*m*m)%mod
  ans+=p*dist dist range(1,n) range 1 n c=n-dist c n-dist n dist p=(c*m*m)%mod p (c*m*m)%mod c*m*m c*m c m m mod ans+=p*dist ans p*dist p dist for dist in range(1,m):
  c=m-dist
  p=(c*n*n)%mod
  ans+=p*dist dist range(1,m) range 1 m c=m-dist c m-dist m dist p=(c*n*n)%mod p (c*n*n)%mod c*n*n c*n c n n mod ans+=p*dist ans p*dist p dist ans%=mod ans mod print(ans*comb(n*m-2,k-2)%mod) print ans*comb(n*m-2,k-2)%mod ans*comb(n*m-2,k-2) ans comb(n*m-2,k-2) comb n*m-2 n*m n m 2 k-2 k 2 mod c=n-dist n-dist c p=(c*n*n)%mod (c*n*n)%mod p m,k=map(int,input().split()) map(int,input().split()) m mod = 10**9+7 10**9+7 mod ans%=mod mod ans c=m-dist m-dist c f += [f[-1]*(i+1) % mod] [f[-1]*(i+1) % mod] f n,m,k=map(int,input().split()) map(int,input().split()) n f = [1] [1] f ans=0 0 ans def comb(a, b):
  return f[a]*pow(f[b], mod-2, mod)*pow(f[a-b], mod-2, mod) % mod def comb(a, b):
  return f[a]*pow(f[b], mod-2, mod)*pow(f[a-b], mod-2, mod) % mod comb ans+=p*dist p*dist ans ans+=p*dist p*dist ans k=map(int,input().split()) map(int,input().split()) k p=(c*m*m)%mod (c*m*m)%mod p
n,m,k=map(int,input().split())

mod = 10**9+7
f = [1]
for i in range(n*m):
  f += [f[-1]*(i+1) % mod]
def comb(a, b):
  return f[a]*pow(f[b], mod-2, mod)*pow(f[a-b], mod-2, mod) % mod

ans=0
for dist in range(1,n):
  c=n-dist
  p=(c*m*m)%mod
  ans+=p*dist
for dist in range(1,m):
  c=m-dist
  p=(c*n*n)%mod
  ans+=p*dist
ans%=mod
print(ans*comb(n*m-2,k-2)%mod)