import numpy as np N, W, K = map(int, input().split()) A = tuple(map(int, input().split())) MOD = 10 ** 9 + 7 def pow_matrix_mod(x, n, mod=MOD): x = x.astype(object) if not n: return np.eye(len(x), dtype=object) if n % 2 == 0: return pow_matrix_mod(x @ x % mod, n // 2) % mod else: return x @ pow_matrix_mod(x @ x % mod, (n - 1) // 2) % mod dp1 = [0] * (2 * W + 1) dp1[0] = 1 for i in range(W + 1): dp1[i] %= MOD x = dp1[i] for a in A: if i + a <= W: dp1[i + a] += x dp2 = [0] * (2 * W + 1) dp2[0] = 1 for i in range(2 * W + 1): dp2[i] %= MOD if i == W: dp2[i] = 0 continue x = dp2[i] for a in A: if i + a <= 2 * W: dp2[i + a] += x matr = np.array([[0, dp2[2 * W]], [1, dp1[W]]], dtype=int) print(pow_matrix_mod(matr, K)[1, 1])