N,M,K = map(int, input().split())
mod = 10**9+7
nm = N*M
n2 = (N**2)%mod
m2 = (M**2)%mod

def cmb(n, r, mod):
    if ( r<0 or r>n ):
        return 0
    r = min(r, n-r)
    return g1[n] * g2[r] * g2[n-r] % mod

g1 = [1, 1] # 元テーブル
g2 = [1, 1] #逆元テーブル
inverse = [0, 1] #逆元テーブル計算用テーブル

for i in range( 2, 10**6 +1):
    g1.append( ( g1[-1] * i ) % mod )
    inverse.append( ( -inverse[mod % i] * (mod//i) ) % mod )
    g2.append( (g2[-1] * inverse[-1]) % mod )


ans = 0
for x in range(1,M):
    ans += x*(M-x)*n2*cmb(nm-2,K-2,mod)
    ans %= mod

for y in range(1,N):
    ans += y*(N-y)*m2*cmb(nm-2,K-2,mod)
    ans %= mod

print(int(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 nm = N*M nm N*M N M n2 = (N**2)%mod n2 (N**2)%mod N**2 N 2 mod m2 = (M**2)%mod m2 (M**2)%mod M**2 M 2 mod def cmb(n, r, mod):
    if ( r<0 or r>n ):
        return 0
    r = min(r, n-r)
    return g1[n] * g2[r] * g2[n-r] % mod cmb if ( r<0 or r>n ):
        return 0
     r<0 or r>n r<0 r 0 r>n r n return 0 0 r = min(r, n-r) r min(r, n-r) min r n-r n r return g1[n] * g2[r] * g2[n-r] % mod g1[n] * g2[r] * g2[n-r] % mod g1[n] * g2[r] * g2[n-r] g1[n] * g2[r] g1[n] g1 n g2[r] g2 r g2[n-r] g2 n-r n r mod n n r r mod mod g1 = [1, 1] g1 [1, 1] 1 1 g2 = [1, 1] g2 [1, 1] 1 1 inverse = [0, 1] inverse [0, 1] 0 1 for i in range( 2, 10**6 +1):
    g1.append( ( g1[-1] * i ) % mod )
    inverse.append( ( -inverse[mod % i] * (mod//i) ) % mod )
    g2.append( (g2[-1] * inverse[-1]) % mod ) i range( 2, 10**6 +1) range 2 10**6 +1 10**6 10 6 1 g1.append( ( g1[-1] * i ) % mod ) g1.append g1 append ( g1[-1] * i ) % mod g1[-1] * i g1[-1] g1 -1 i mod inverse.append( ( -inverse[mod % i] * (mod//i) ) % mod ) inverse.append inverse append ( -inverse[mod % i] * (mod//i) ) % mod -inverse[mod % i] * (mod//i) -inverse[mod % i] inverse[mod % i] inverse mod % i mod i mod//i mod i mod g2.append( (g2[-1] * inverse[-1]) % mod ) g2.append g2 append (g2[-1] * inverse[-1]) % mod g2[-1] * inverse[-1] g2[-1] g2 -1 inverse[-1] inverse -1 mod ans = 0 ans 0 for x in range(1,M):
    ans += x*(M-x)*n2*cmb(nm-2,K-2,mod)
    ans %= mod x range(1,M) range 1 M ans += x*(M-x)*n2*cmb(nm-2,K-2,mod) ans x*(M-x)*n2*cmb(nm-2,K-2,mod) x*(M-x)*n2 x*(M-x) x M-x M x n2 cmb(nm-2,K-2,mod) cmb nm-2 nm 2 K-2 K 2 mod ans %= mod ans mod for y in range(1,N):
    ans += y*(N-y)*m2*cmb(nm-2,K-2,mod)
    ans %= mod y range(1,N) range 1 N ans += y*(N-y)*m2*cmb(nm-2,K-2,mod) ans y*(N-y)*m2*cmb(nm-2,K-2,mod) y*(N-y)*m2 y*(N-y) y N-y N y m2 cmb(nm-2,K-2,mod) cmb nm-2 nm 2 K-2 K 2 mod ans %= mod ans mod print(int(ans)) print int(ans) int ans ans %= mod mod ans g1 = [1, 1] [1, 1] g1 def cmb(n, r, mod):
    if ( r<0 or r>n ):
        return 0
    r = min(r, n-r)
    return g1[n] * g2[r] * g2[n-r] % mod def cmb(n, r, mod):
    if ( r<0 or r>n ):
        return 0
    r = min(r, n-r)
    return g1[n] * g2[r] * g2[n-r] % mod cmb ans += y*(N-y)*m2*cmb(nm-2,K-2,mod) y*(N-y)*m2*cmb(nm-2,K-2,mod) ans inverse = [0, 1] [0, 1] inverse g2 = [1, 1] [1, 1] g2 N,M,K = map(int, input().split()) map(int, input().split()) N nm = N*M N*M nm ans %= mod mod ans K = map(int, input().split()) map(int, input().split()) K M,K = map(int, input().split()) map(int, input().split()) M mod = 10**9+7 10**9+7 mod n2 = (N**2)%mod (N**2)%mod n2 m2 = (M**2)%mod (M**2)%mod m2 ans = 0 0 ans ans += x*(M-x)*n2*cmb(nm-2,K-2,mod) x*(M-x)*n2*cmb(nm-2,K-2,mod) ans
N,M,K = map(int, input().split())
mod = 10**9+7
nm = N*M
n2 = (N**2)%mod
m2 = (M**2)%mod

def cmb(n, r, mod):
    if ( r<0 or r>n ):
        return 0
    r = min(r, n-r)
    return g1[n] * g2[r] * g2[n-r] % mod

g1 = [1, 1] # 元テーブル
g2 = [1, 1] #逆元テーブル
inverse = [0, 1] #逆元テーブル計算用テーブル

for i in range( 2, 10**6 +1):
    g1.append( ( g1[-1] * i ) % mod )
    inverse.append( ( -inverse[mod % i] * (mod//i) ) % mod )
    g2.append( (g2[-1] * inverse[-1]) % mod )


ans = 0
for x in range(1,M):
    ans += x*(M-x)*n2*cmb(nm-2,K-2,mod)
    ans %= mod

for y in range(1,N):
    ans += y*(N-y)*m2*cmb(nm-2,K-2,mod)
    ans %= mod

print(int(ans))