N, M, K = map(int, input().split())

ans = [0] * (K + 1)

if N == 1:
    ans = [i * pow(i - 1, M - 1, 1000000007) % 1000000007 for i in range(K + 1)]
elif N == 2:
    ans = [i * (i - 1) * pow((i - 3) * i + 3, M - 1, 1000000007) % 1000000007 for i in range(K + 1)]
else:
    import numpy as np
    ans = np.array([(i - 2) * ((i - 3) * i + 5) % 1000000007 for i in range(K + 1)])
    b = np.array([(((i - 6) * i + 13) * i - 11) % 1000000007 * (i - 1) % 1000000007 for i in range(K + 1)])
    c = np.array([i * (i - 1) * (i - 1) % 1000000007 for i in range(K + 1)])
    d = np.array([((i - 3) * i + 1) % 1000000007 * i * (i - 1) % 1000000007 for i in range(K + 1)])
    n = M - 1
    while n:
        if n & 1:
            c = ans * c - d
            c = c % 1000000007
            d = b * d
        else:
            d = ans * d - b * c
        d = d % 1000000007
        ans = ans * ans - 2 * b
        ans = ans % 1000000007
        b = b * b
        b = b % 1000000007
        n >>= 1
    ans = c.tolist()

coef = [1] * (K + 1)

for i in range(K):
    coef[i + 1] = coef[i] * (K - i) % 1000000007

c = 1
for i in range(K + 1, 0, -1):
    c = -c * i % 1000000007
    coef[i - 1] = coef[i - 1] * c % 1000000007

ret = 0
for i in range(K + 1):
    ret += ans[i] * coef[i]

ret *= pow(coef[-1], 1000000005, 1000000007)
print(ret % 1000000007)