N, M, K = map(int, input().split()) ans = [0] * (K + 1) if N == 1: ans = [i * pow(i, 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)