#pragma GCC optimize ("O3") #pragma GCC target ("tune=native") #pragma GCC target ("avx") #include // 汎用マクロ #define ALL_OF(x) (x).begin(), (x).end() #define REP(i,n) for (long long i=0, i##_len=(n); i=i##_end; i--) #define UNIQUE(v) { sort((v).begin(), (v).end()); (v).erase(unique((v).begin(), (v).end()), (v).end()); } template bool chmax(T &a, const T &b) {if (a < b) {a = b; return true;} return false; } template bool chmin(T &a, const T &b) {if (a > b) {a = b; return true;} return false; } #define INF 0x7FFFFFFF #define LINF 0x7FFFFFFFFFFFFFFFLL #define Yes(q) (q ? "Yes" : "No") #define YES(q) (q ? "YES" : "NO") #define DUMP(q) cerr << "[DEBUG] " #q ": " << (q) << " at " __FILE__ ":" << __LINE__ << endl #define DUMPALL(q) { cerr << "[DEBUG] " #q ": ["; REP(dumpall_i, (q).size()) { cerr << q[dumpall_i] << (dumpall_i == (q).size() - 1 ? "" : ", "); } cerr << "] at " __FILE__ ":" << __LINE__ << endl; } template T gcd(const T &a, const T &b) { return a < b ? gcd(b, a) : b ? gcd(b, a % b) : a; } template T lcm(const T &a, const T &b) { return a / gcd(a, b) * b; } // gcc拡張マクロ #define popcount __builtin_popcount #define popcountll __builtin_popcountll // エイリアス using ll = long long; using ull = unsigned long long; using ld = long double; using namespace std; // モジュール // 剰余演算 constexpr ll p = 1000000007LL; // 10^9+7 ll powll(ll a, ll n) { if (n == 0) { return 1LL; } else if (n == 1) { return a % p; } else { ll temp = powll(a, n/2); temp = temp * temp % p; return n % 2 ? temp * a % p : temp; } } inline ll invll(ll a) { return powll(a, p-2); } vector fac_cache = {1}, invfac_cache = {1}; void make_fac_cache(ll a) { ll old_max = fac_cache.size() - 1; if (a > old_max) { fac_cache .resize(a+1); invfac_cache.resize(a+1); for (ll i = old_max + 1; i <= a; i++) { fac_cache[i] = fac_cache[i-1] * i % p; } invfac_cache[a] = invll(fac_cache[a]); for (ll i = a-1; i > old_max; i--) { invfac_cache[i] = invfac_cache[i+1] * (i + 1) % p; } } } inline ll facll(ll a) { make_fac_cache(a); return fac_cache[a]; } inline ll invfacll(ll a) { make_fac_cache(a); return invfac_cache[a]; } inline bool isoutll(ll n, ll r) { return n < 0 || r < 0 || n < r; } inline ll nPr(ll n, ll r) { return isoutll(n, r) ? 0 : facll(n) * invfacll(n-r) % p; } inline ll nCr(ll n, ll r) { return isoutll(n, r) ? 0 : facll(n) * invfacll(n-r) % p * invfacll(r) % p; } // 行列累乗; 1e9+7で割ったあまりを求めたい場合は適宜コメントアウトを解除 // 単位行列; idtm<型>(大きさ) template vector> idtm(const size_t n) { vector> r(n, vector(n, 0)); REP(i, n) r[i][i] = 1; return r; } // 行列の積 template vector> mltm(const vector> &a, const vector> &b) { ull m = a.size(), q = b.size(), r = b[0].size(); vector> c(m, vector(r, (T)0)); REP(i, m) REP(j, r) { REP(k, m) { // c[i][j] += a[i][k] * b[k][j]; c[i][j] += a[i][k] * b[k][j] % p; c[i][j] %= p; } } return c; } // 行列の累乗 template vector> powm(const vector> &a, const ll n) { if (n == 0) { return idtm(a.size()); } else if (n == 1) { return a; } else if (n % 2 == 0) { vector> t = powm(a, n/2); return mltm(t, t); } else { vector> t = powm(a, n/2); return mltm(mltm(t, t), a); } } // 処理内容 int main() { clock_t sclk = clock(); ll n, m, k, pp, q; cin >> n >> m >> k >> pp >> q; vector b(n); REP(i, n) cin >> b[i]; ll bmlk = 0, bktc = 0; REP(i, n) { if (i < m){ bmlk += b[i]; bmlk %= p; } else { bktc += b[i]; bktc %= p; } } ll pq2 = pp * invll(q) % p * 2 % p; ll half = invll(2); ll term = (1 + p - pq2) % p; ll pmlk = (1 + powll(term, k) ) % p * half % p; ll pktc = (1 - powll(term, k) + p) % p * half % p; // DUMP(n); // DUMP(m); // DUMP(k); // DUMP(pp); // DUMP(q); // DUMP(bmlk); // DUMP(bktc); // DUMP(pq2); // DUMP(pmlk); // DUMP(pktc); cout << (bmlk * pmlk % p + bktc * pktc % p) % p << endl; cerr << "[DEBUG] Done in " << (double)(clock() - sclk) / CLOCKS_PER_SEC * 1e3 << " ms" << endl; }