import sys
input = sys.stdin.buffer.readline
sys.setrecursionlimit(10 ** 7)
mod = 10 ** 9 + 7

N, M, K, a, b = map(int, input().split())
B = [int(input()) for _ in range(N)]

p = a * pow(b, mod - 2, mod) % mod
q = (b - a) * pow(b, mod - 2, mod) % mod
warp1 = (1 + pow(q - p, K, mod)) * pow(2, mod - 2, mod) % mod
warp2 = (1 - pow(q - p, K, mod) + mod) * pow(2, mod - 2, mod) % mod

ans = 0
for i in range(M):
    ans += B[i] * warp1
    ans %= mod

for i in range(M, N):
    ans += B[i] * warp2
    ans %= mod

print(ans)