N, M, K = map(int, input().split())

MOD = 10**9+7
n = N*M-2
fac = [1]*(n+1)
rev = [1]*(n+1)
 
for i in range(1,n+1):
  fac[i] = i*fac[i-1]%MOD
  rev[i] = pow(fac[i], MOD-2, MOD)
 
comb = lambda a,b:(fac[a]*rev[a-b]*rev[b])%MOD

C = comb(N*M-2, K-2)
ans = 0
for d in range(N):
    ans += d*(N-d)*(M**2)*C
    ans %= MOD

for d in range(M):
    ans += d*(M-d)*(N**2)*C
    ans %= MOD

print(ans) 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 n = N*M-2 n N*M-2 N*M N M 2 fac = [1]*(n+1) fac [1]*(n+1) [1] 1 n+1 n 1 rev = [1]*(n+1) rev [1]*(n+1) [1] 1 n+1 n 1 for i in range(1,n+1):
  fac[i] = i*fac[i-1]%MOD
  rev[i] = pow(fac[i], MOD-2, MOD)
  i range(1,n+1) range 1 n+1 n 1 fac[i] = i*fac[i-1]%MOD fac[i] fac i i*fac[i-1]%MOD i*fac[i-1] i fac[i-1] fac i-1 i 1 MOD rev[i] = pow(fac[i], MOD-2, MOD) rev[i] rev i pow(fac[i], MOD-2, MOD) pow fac[i] fac i MOD-2 MOD 2 MOD comb = lambda a,b:(fac[a]*rev[a-b]*rev[b])%MOD comb lambda a,b:(fac[a]*rev[a-b]*rev[b])%MOD (fac[a]*rev[a-b]*rev[b])%MOD fac[a]*rev[a-b]*rev[b] fac[a]*rev[a-b] fac[a] fac a rev[a-b] rev a-b a b rev[b] rev b MOD a b C = comb(N*M-2, K-2) C comb(N*M-2, K-2) comb N*M-2 N*M N M 2 K-2 K 2 ans = 0 ans 0 for d in range(N):
    ans += d*(N-d)*(M**2)*C
    ans %= MOD d range(N) range N ans += d*(N-d)*(M**2)*C ans d*(N-d)*(M**2)*C d*(N-d)*(M**2) d*(N-d) d N-d N d M**2 M 2 C ans %= MOD ans MOD for d in range(M):
    ans += d*(M-d)*(N**2)*C
    ans %= MOD d range(M) range M ans += d*(M-d)*(N**2)*C ans d*(M-d)*(N**2)*C d*(M-d)*(N**2) d*(M-d) d M-d M d N**2 N 2 C ans %= MOD ans MOD print(ans) print ans ans += d*(N-d)*(M**2)*C d*(N-d)*(M**2)*C ans ans += d*(M-d)*(N**2)*C d*(M-d)*(N**2)*C ans fac = [1]*(n+1) [1]*(n+1) fac n = N*M-2 N*M-2 n C = comb(N*M-2, K-2) comb(N*M-2, K-2) C M, K = map(int, input().split()) map(int, input().split()) M comb = lambda a,b:(fac[a]*rev[a-b]*rev[b])%MOD lambda a,b:(fac[a]*rev[a-b]*rev[b])%MOD comb K = map(int, input().split()) map(int, input().split()) K rev = [1]*(n+1) [1]*(n+1) rev ans = 0 0 ans ans %= MOD MOD ans MOD = 10**9+7 10**9+7 MOD ans %= MOD MOD ans N, M, K = map(int, input().split()) map(int, input().split()) N
N, M, K = map(int, input().split())

MOD = 10**9+7
n = N*M-2
fac = [1]*(n+1)
rev = [1]*(n+1)
 
for i in range(1,n+1):
  fac[i] = i*fac[i-1]%MOD
  rev[i] = pow(fac[i], MOD-2, MOD)
 
comb = lambda a,b:(fac[a]*rev[a-b]*rev[b])%MOD

C = comb(N*M-2, K-2)
ans = 0
for d in range(N):
    ans += d*(N-d)*(M**2)*C
    ans %= MOD

for d in range(M):
    ans += d*(M-d)*(N**2)*C
    ans %= MOD

print(ans)