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

def comb(n, k):
    ret = 1
    for i in range(1, k + 1):
        ret = ret * (n - i + 1) % MOD        
        ret = ret * pow(i, MOD - 2, MOD) % MOD
        
    return ret

ans = 0
for d in range(N):
    ans += d * (N - d) * M ** 2
    ans %= MOD
    
for d in range(M):
    ans += d * (M - d) * N ** 2
    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 def comb(n, k):
    ret = 1
    for i in range(1, k + 1):
        ret = ret * (n - i + 1) % MOD        
        ret = ret * pow(i, MOD - 2, MOD) % MOD
        
    return ret comb ret = 1 ret 1 for i in range(1, k + 1):
        ret = ret * (n - i + 1) % MOD        
        ret = ret * pow(i, MOD - 2, MOD) % MOD
        
     i range(1, k + 1) range 1 k + 1 k 1 ret = ret * (n - i + 1) % MOD ret ret * (n - i + 1) % MOD ret * (n - i + 1) ret n - i + 1 n - i n i 1 MOD ret = ret * pow(i, MOD - 2, MOD) % MOD ret ret * pow(i, MOD - 2, MOD) % MOD ret * pow(i, MOD - 2, MOD) ret pow(i, MOD - 2, MOD) pow i MOD - 2 MOD 2 MOD MOD return ret ret n n k k ans = 0 ans 0 for d in range(N):
    ans += d * (N - d) * M ** 2
    ans %= MOD
     d range(N) range N ans += d * (N - d) * M ** 2 ans d * (N - d) * M ** 2 d * (N - d) d N - d N d M ** 2 M 2 ans %= MOD ans MOD for d in range(M):
    ans += d * (M - d) * N ** 2
    ans %= MOD
     d range(M) range M ans += d * (M - d) * N ** 2 ans d * (M - d) * N ** 2 d * (M - d) d M - d M d N ** 2 N 2 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 K = map(int, input().split()) map(int, input().split()) K ans += d * (M - d) * N ** 2 d * (M - d) * N ** 2 ans ans %= MOD MOD ans ans = 0 0 ans N, M, K = map(int, input().split()) map(int, input().split()) N ans += d * (N - d) * M ** 2 d * (N - d) * M ** 2 ans ans %= MOD MOD ans M, K = map(int, input().split()) map(int, input().split()) M def comb(n, k):
    ret = 1
    for i in range(1, k + 1):
        ret = ret * (n - i + 1) % MOD        
        ret = ret * pow(i, MOD - 2, MOD) % MOD
        
    return ret def comb(n, k):
    ret = 1
    for i in range(1, k + 1):
        ret = ret * (n - i + 1) % MOD        
        ret = ret * pow(i, MOD - 2, MOD) % MOD
        
    return ret comb MOD = 10 ** 9 + 7 10 ** 9 + 7 MOD
N, M, K = map(int, input().split())
MOD = 10 ** 9 + 7

def comb(n, k):
    ret = 1
    for i in range(1, k + 1):
        ret = ret * (n - i + 1) % MOD        
        ret = ret * pow(i, MOD - 2, MOD) % MOD
        
    return ret

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

