def com(n,r):
    re=fct[n]*pow(fct[r],md-2,md)*pow(fct[n-r],md-2,md)
    re%=md
    return re

n,m,k=map(int,input().split())
md=10**9+7
ans=0
fct=[1]*(n*m-1)
f=1
for i in range(1,n*m-1):
    f=f*i%md
    fct[i]=f
for d in range(1,m):
    ans=(ans+d*(m-d)*n*n)%md
for d in range(1,n):
    ans=(ans+d*(n-d)*m*m)%md
ans=(ans*com(n*m-2,k-2))%md
print(ans) def com(n,r):
    re=fct[n]*pow(fct[r],md-2,md)*pow(fct[n-r],md-2,md)
    re%=md
    return re com re=fct[n]*pow(fct[r],md-2,md)*pow(fct[n-r],md-2,md) re fct[n]*pow(fct[r],md-2,md)*pow(fct[n-r],md-2,md) fct[n]*pow(fct[r],md-2,md) fct[n] fct n pow(fct[r],md-2,md) pow fct[r] fct r md-2 md 2 md pow(fct[n-r],md-2,md) pow fct[n-r] fct n-r n r md-2 md 2 md re%=md re md return re re n n r r n,m,k=map(int,input().split()) n map(int,input().split()) map int input().split() ().split () input split m k md=10**9+7 md 10**9+7 10**9 10 9 7 ans=0 ans 0 fct=[1]*(n*m-1) fct [1]*(n*m-1) [1] 1 n*m-1 n*m n m 1 f=1 f 1 for i in range(1,n*m-1):
    f=f*i%md
    fct[i]=f i range(1,n*m-1) range 1 n*m-1 n*m n m 1 f=f*i%md f f*i%md f*i f i md fct[i]=f fct[i] fct i f for d in range(1,m):
    ans=(ans+d*(m-d)*n*n)%md d range(1,m) range 1 m ans=(ans+d*(m-d)*n*n)%md ans (ans+d*(m-d)*n*n)%md ans+d*(m-d)*n*n ans d*(m-d)*n*n d*(m-d)*n d*(m-d) d m-d m d n n md for d in range(1,n):
    ans=(ans+d*(n-d)*m*m)%md d range(1,n) range 1 n ans=(ans+d*(n-d)*m*m)%md ans (ans+d*(n-d)*m*m)%md ans+d*(n-d)*m*m ans d*(n-d)*m*m d*(n-d)*m d*(n-d) d n-d n d m m md ans=(ans*com(n*m-2,k-2))%md ans (ans*com(n*m-2,k-2))%md ans*com(n*m-2,k-2) ans com(n*m-2,k-2) com n*m-2 n*m n m 2 k-2 k 2 md print(ans) print ans ans=0 0 ans ans=(ans*com(n*m-2,k-2))%md (ans*com(n*m-2,k-2))%md ans ans=(ans+d*(n-d)*m*m)%md (ans+d*(n-d)*m*m)%md ans fct=[1]*(n*m-1) [1]*(n*m-1) fct ans=(ans+d*(m-d)*n*n)%md (ans+d*(m-d)*n*n)%md ans def com(n,r):
    re=fct[n]*pow(fct[r],md-2,md)*pow(fct[n-r],md-2,md)
    re%=md
    return re def com(n,r):
    re=fct[n]*pow(fct[r],md-2,md)*pow(fct[n-r],md-2,md)
    re%=md
    return re com n,m,k=map(int,input().split()) map(int,input().split()) n k=map(int,input().split()) map(int,input().split()) k md=10**9+7 10**9+7 md f=f*i%md f*i%md f m,k=map(int,input().split()) map(int,input().split()) m f=1 1 f
def com(n,r):
    re=fct[n]*pow(fct[r],md-2,md)*pow(fct[n-r],md-2,md)
    re%=md
    return re

n,m,k=map(int,input().split())
md=10**9+7
ans=0
fct=[1]*(n*m-1)
f=1
for i in range(1,n*m-1):
    f=f*i%md
    fct[i]=f
for d in range(1,m):
    ans=(ans+d*(m-d)*n*n)%md
for d in range(1,n):
    ans=(ans+d*(n-d)*m*m)%md
ans=(ans*com(n*m-2,k-2))%md
print(ans)
