# coding: utf-8
# Your code here!
import sys
sys.setrecursionlimit(10**6)
readline = sys.stdin.readline
read = sys.stdin.read


p,n,k,b,*a = [int(i) for i in read().split()]

lst = [0]*p
for i in range(p): lst[pow(i,k,p)] += 1

MOD = 10**9+7

x = 1
M = 1<<96
NUM = 96*p
ALL = (1<<NUM)-1
m = mm = ((1<<32)-1)<<64
for _ in range(n):
    m <<= 96
    m ^= mm
pa = (1 << 64) - ((1 << 64) % MOD)

for ai in a:
    y = 0
    for i,c in enumerate(lst):
        y += c<<(i*ai%p*96)
    x = x*y
    x = (x>>NUM) + (x&ALL)
    x -= ((x & m) >> 64) * pa

ans = (x>>(b*96))&((1<<96)-1)
print(ans%MOD)