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

def FastConvolution(A,B): # A,B: 2^c 成分を想定
    return ifft(fft(A)*fft(B))

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 = int(("1" * 2 + "0" * 30) * p, 2)
pa = (1 << 30) - ((1 << 30) % MOD)

for ai in a:
    y = 0
    for i,c in enumerate(lst):
        if c!=0:
            y += c<<(i*ai%p*96)
    #print(bin(y)[2:])
    x = x*y
    x = (x>>(NUM)) + (x&ALL)
    x -= ((x & m) >> 96) * pa
    

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