#include #include #include #define llint long long #define mod 1000000007 using namespace std; llint modpow(llint a, llint n) { if(n == 0) return 1; if(n % 2){ return ((a%mod) * (modpow(a, n-1)%mod)) % mod; } else{ return modpow((a*a)%mod, n/2) % mod; } } struct Matrix{ int h, w; vector mat; Matrix(){ h = w = 0; } Matrix(int h, int w){ this->h = h, this->w = w; mat.resize(h*w); } llint& at(int i, int j){ return mat[w*(i-1)+(j-1)]; } static Matrix ident(int size){ Matrix ret(size, size); for(int i = 1; i <= size; i++) ret.at(i, i) = 1; return ret; } Matrix operator*(Matrix& ope){ Matrix ret(h, ope.w); for(int i = 1; i <= h; i++){ for(int j = 1; j <= ope.w; j++){ for(int k = 1; k <= w; k++){ ret.at(i, j) += at(i, k) * ope.at(k, j) % mod; ret.at(i, j) %= mod; } } } return ret; } Matrix pow(llint n){ if(n == 0) return ident(h); if(n % 2){ return pow(n-1) * (*this); } else{ Matrix tmp = pow(n/2); return tmp * tmp; } } void print(){ for(int i = 1; i <= h; i++){ for(int j = 1; j <= w; j++){ cout << at(i, j) << " "; } cout << endl; } } }; llint n, w, k; llint a[105]; llint dp[105], dp2[105]; llint x1, x2; Matrix mat; int main(void) { ios::sync_with_stdio(0); cin.tie(0); cin >> n >> w >> k; for(int i = 1; i <= n; i++) cin >> a[i]; dp[0] = 1; for(int i = 0; i < w; i++){ for(int j = 1; j <= n; j++){ if(i+a[j] <= w) dp[i+a[j]] += dp[i], dp[i+a[j]] %= mod; } } x1 = dp[w]; dp2[0] = 1; for(int i = 0; i < 2*w; i++){ if(i == w) continue; for(int j = 1; j <= n; j++){ if(i+a[j] <= 2*w) dp2[i+a[j]] += dp2[i], dp2[i+a[j]] %= mod; } } x2 = dp2[2*w]; mat = Matrix(2, 2); mat.at(1, 1) = x1, mat.at(1, 2) = x2; mat.at(2, 1) = 1, mat.at(2, 2) = 0; Matrix pmat = mat.pow(k); cout << pmat.at(1, 1) << endl; return 0; }