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

mod = 10**9 + 7

def chs(a,b, md):
    res = 1
    for i in range(a-b+1,a+1):
        res *=  i
        res %= md
    rev = 1
    for i in range(1,b+1):
        rev *= i
        rev %= md
    rev2 = pow(rev,md-2,md)
    
    res *= rev2
    res %= md
    
    return res

ans = 0
c = chs(n*m-2,k-2,mod)

for i in range(1,n):
    ans += i*(n-i)*m**2
    ans %= mod

        

for i in range(1,m):
    ans +=i*(m-i)*n**2
    ans %=mod
ans *= 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 def chs(a,b, md):
    res = 1
    for i in range(a-b+1,a+1):
        res *=  i
        res %= md
    rev = 1
    for i in range(1,b+1):
        rev *= i
        rev %= md
    rev2 = pow(rev,md-2,md)
    
    res *= rev2
    res %= md
    
    return res chs res = 1 res 1 for i in range(a-b+1,a+1):
        res *=  i
        res %= md
     i range(a-b+1,a+1) range a-b+1 a-b a b 1 a+1 a 1 res *=  i res i res %= md res md rev = 1 rev 1 for i in range(1,b+1):
        rev *= i
        rev %= md
     i range(1,b+1) range 1 b+1 b 1 rev *= i rev i rev %= md rev md rev2 = pow(rev,md-2,md) rev2 pow(rev,md-2,md) pow rev md-2 md 2 md res *= rev2 res rev2 res %= md res md return res res a a b b md md ans = 0 ans 0 c = chs(n*m-2,k-2,mod) c chs(n*m-2,k-2,mod) chs n*m-2 n*m n m 2 k-2 k 2 mod for i in range(1,n):
    ans += i*(n-i)*m**2
    ans %= mod

         i range(1,n) range 1 n ans += i*(n-i)*m**2 ans i*(n-i)*m**2 i*(n-i) i n-i n i m**2 m 2 ans %= mod ans mod for i in range(1,m):
    ans +=i*(m-i)*n**2
    ans %=mod i range(1,m) range 1 m ans +=i*(m-i)*n**2 ans i*(m-i)*n**2 i*(m-i) i m-i m i n**2 n 2 ans %=mod ans mod ans *= c ans c ans %= mod ans mod print(ans) print ans ans = 0 0 ans ans += i*(n-i)*m**2 i*(n-i)*m**2 ans m, k = map(int, input().split( )) map(int, input().split( )) m def chs(a,b, md):
    res = 1
    for i in range(a-b+1,a+1):
        res *=  i
        res %= md
    rev = 1
    for i in range(1,b+1):
        rev *= i
        rev %= md
    rev2 = pow(rev,md-2,md)
    
    res *= rev2
    res %= md
    
    return res def chs(a,b, md):
    res = 1
    for i in range(a-b+1,a+1):
        res *=  i
        res %= md
    rev = 1
    for i in range(1,b+1):
        rev *= i
        rev %= md
    rev2 = pow(rev,md-2,md)
    
    res *= rev2
    res %= md
    
    return res chs k = map(int, input().split( )) map(int, input().split( )) k ans +=i*(m-i)*n**2 i*(m-i)*n**2 ans ans %= mod mod ans mod = 10**9 + 7 10**9 + 7 mod c = chs(n*m-2,k-2,mod) chs(n*m-2,k-2,mod) c ans *= c c ans ans %=mod mod ans n, m, k = map(int, input().split( )) map(int, input().split( )) n ans %= mod mod ans
n, m, k = map(int, input().split( ))

mod = 10**9 + 7

def chs(a,b, md):
    res = 1
    for i in range(a-b+1,a+1):
        res *=  i
        res %= md
    rev = 1
    for i in range(1,b+1):
        rev *= i
        rev %= md
    rev2 = pow(rev,md-2,md)
    
    res *= rev2
    res %= md
    
    return res

ans = 0
c = chs(n*m-2,k-2,mod)

for i in range(1,n):
    ans += i*(n-i)*m**2
    ans %= mod

        

for i in range(1,m):
    ans +=i*(m-i)*n**2
    ans %=mod
ans *= c
ans %= mod
print(ans)
