n,m,k=map(int,input().split())
mod=10**9+7
x=1
y=1
for i in range(1,k-1):
    y*=i
    y%=mod
for i in range(k-2):
    x*=n*m-2-i
    x%=mod
x*=pow(y,mod-2,mod)
x%=mod
def f(n,m):
    s=0
    for i in range(1,n):
        s+=i*(n-i)
        s%=mod
    return s*m*m%mod
print((f(n,m)+f(m,n))*x%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 x=1 x 1 y=1 y 1 for i in range(1,k-1):
    y*=i
    y%=mod i range(1,k-1) range 1 k-1 k 1 y*=i y i y%=mod y mod for i in range(k-2):
    x*=n*m-2-i
    x%=mod i range(k-2) range k-2 k 2 x*=n*m-2-i x n*m-2-i n*m-2 n*m n m 2 i x%=mod x mod x*=pow(y,mod-2,mod) x pow(y,mod-2,mod) pow y mod-2 mod 2 mod x%=mod x mod def f(n,m):
    s=0
    for i in range(1,n):
        s+=i*(n-i)
        s%=mod
    return s*m*m%mod f s=0 s 0 for i in range(1,n):
        s+=i*(n-i)
        s%=mod
     i range(1,n) range 1 n s+=i*(n-i) s i*(n-i) i n-i n i s%=mod s mod return s*m*m%mod s*m*m%mod s*m*m s*m s m m mod n n m m print((f(n,m)+f(m,n))*x%mod) print (f(n,m)+f(m,n))*x%mod (f(n,m)+f(m,n))*x f(n,m)+f(m,n) f(n,m) f n m f(m,n) f m n x mod mod=10**9+7 10**9+7 mod x*=n*m-2-i n*m-2-i x k=map(int,input().split()) map(int,input().split()) k m,k=map(int,input().split()) map(int,input().split()) m y%=mod mod y x%=mod mod x x%=mod mod x def f(n,m):
    s=0
    for i in range(1,n):
        s+=i*(n-i)
        s%=mod
    return s*m*m%mod def f(n,m):
    s=0
    for i in range(1,n):
        s+=i*(n-i)
        s%=mod
    return s*m*m%mod f x*=pow(y,mod-2,mod) pow(y,mod-2,mod) x n,m,k=map(int,input().split()) map(int,input().split()) n y*=i i y y=1 1 y x=1 1 x
n,m,k=map(int,input().split())
mod=10**9+7
x=1
y=1
for i in range(1,k-1):
    y*=i
    y%=mod
for i in range(k-2):
    x*=n*m-2-i
    x%=mod
x*=pow(y,mod-2,mod)
x%=mod
def f(n,m):
    s=0
    for i in range(1,n):
        s+=i*(n-i)
        s%=mod
    return s*m*m%mod
print((f(n,m)+f(m,n))*x%mod)