MOD = 10**9 + 7

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

invs = [1] * (K-1)
for x in range(2, K-1):
    invs[x] = (-(MOD//x) * invs[MOD%x]) % MOD

C = 1
for x in range(N*M-K+1, N*M-1):
    C *= x
    C %= MOD
for x in range(1, K-1):
    C *= invs[x]
    C %= MOD

ans1 = 0
for d in range(1, N):
    ans1 += (N-d) * d
    ans1 %= MOD
ans1 *= M**2
ans1 %= MOD

ans2 = 0
for d in range(1, M):
    ans2 += (M-d) * d
    ans2 %= MOD
ans2 *= N**2
ans2 %= MOD

print((ans1+ans2) * C % MOD) 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 invs = [1] * (K-1) invs [1] * (K-1) [1] 1 K-1 K 1 for x in range(2, K-1):
    invs[x] = (-(MOD//x) * invs[MOD%x]) % MOD x range(2, K-1) range 2 K-1 K 1 invs[x] = (-(MOD//x) * invs[MOD%x]) % MOD invs[x] invs x (-(MOD//x) * invs[MOD%x]) % MOD -(MOD//x) * invs[MOD%x] -(MOD//x) MOD//x MOD x invs[MOD%x] invs MOD%x MOD x MOD C = 1 C 1 for x in range(N*M-K+1, N*M-1):
    C *= x
    C %= MOD x range(N*M-K+1, N*M-1) range N*M-K+1 N*M-K N*M N M K 1 N*M-1 N*M N M 1 C *= x C x C %= MOD C MOD for x in range(1, K-1):
    C *= invs[x]
    C %= MOD x range(1, K-1) range 1 K-1 K 1 C *= invs[x] C invs[x] invs x C %= MOD C MOD ans1 = 0 ans1 0 for d in range(1, N):
    ans1 += (N-d) * d
    ans1 %= MOD d range(1, N) range 1 N ans1 += (N-d) * d ans1 (N-d) * d N-d N d d ans1 %= MOD ans1 MOD ans1 *= M**2 ans1 M**2 M 2 ans1 %= MOD ans1 MOD ans2 = 0 ans2 0 for d in range(1, M):
    ans2 += (M-d) * d
    ans2 %= MOD d range(1, M) range 1 M ans2 += (M-d) * d ans2 (M-d) * d M-d M d d ans2 %= MOD ans2 MOD ans2 *= N**2 ans2 N**2 N 2 ans2 %= MOD ans2 MOD print((ans1+ans2) * C % MOD) print (ans1+ans2) * C % MOD (ans1+ans2) * C ans1+ans2 ans1 ans2 C MOD MOD = 10**9 + 7 10**9 + 7 MOD C %= MOD MOD C C *= invs[x] invs[x] C ans1 = 0 0 ans1 ans1 %= MOD MOD ans1 C = 1 1 C invs = [1] * (K-1) [1] * (K-1) invs ans1 %= MOD MOD ans1 ans2 = 0 0 ans2 ans2 *= N**2 N**2 ans2 C %= MOD MOD C C *= x x C ans2 %= MOD MOD ans2 ans1 *= M**2 M**2 ans1 ans2 %= MOD MOD ans2 N, M, K = map(int, input().split()) map(int, input().split()) N ans2 += (M-d) * d (M-d) * d ans2 M, K = map(int, input().split()) map(int, input().split()) M K = map(int, input().split()) map(int, input().split()) K ans1 += (N-d) * d (N-d) * d ans1
MOD = 10**9 + 7

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

invs = [1] * (K-1)
for x in range(2, K-1):
    invs[x] = (-(MOD//x) * invs[MOD%x]) % MOD

C = 1
for x in range(N*M-K+1, N*M-1):
    C *= x
    C %= MOD
for x in range(1, K-1):
    C *= invs[x]
    C %= MOD

ans1 = 0
for d in range(1, N):
    ans1 += (N-d) * d
    ans1 %= MOD
ans1 *= M**2
ans1 %= MOD

ans2 = 0
for d in range(1, M):
    ans2 += (M-d) * d
    ans2 %= MOD
ans2 *= N**2
ans2 %= MOD

print((ans1+ans2) * C % MOD)
