Description
小C有一个集合S，里面的元素都是小于M的非负整数。他用程序编写了一个数列生成器，可以生成一个长度为N的数列，数列中的每个数都属于集合S。
小C用这个生成器生成了许多这样的数列。但是小C有一个问题需要你的帮助：给定整数x，求所有可以生成出的，且满足数列中所有数的乘积mod M的值等于x的不同的数列的有多少个。小C认为，两个数列{Ai}和{Bi}不同，当且仅当至少存在一个整数i，满足Ai≠Bi。另外，小C认为这个问题的答案可能很大，因此他只需要你帮助他求出答案mod 1004535809的值就可以了。
Input
一行，四个整数，N、M、x、|S|，其中|S|为集合S中元素个数。第二行，|S|个整数，表示集合S中的所有元素。
Output
一行，一个整数，表示你求出的种类数mod 1004535809的值。
Sample Input
4 3 1 2
1 2
1 2
Sample Output
8
HINT
【样例说明】
可以生成的满足要求的不同的数列有(1,1,1,1)、(1,1,2,2)、(1,2,1,2)、(1,2,2,1)、(2,1,1,2)、(2,1,2,1)、(2,2,1,1)、(2,2,2,2)。
【数据规模和约定】
对于10%的数据，1<=N<=1000；
对于30%的数据，3<=M<=100；
对于60%的数据，3<=M<=800；
对于全部的数据，1<=N<=10^9，3<=M<=8000，M为质数，1<=x<=M-1，输入数据保证集合S中元素不重复
好厉害。。。
首先有10分算法，设f[i][j]表示选i个数，乘积模M结果为x的方案数。
然后因为M为质数，我们可以求出M的原根g，这样转移就可以写成关于g的生成函数，就可以用NTT来加速了。
接下来发现每次乘的都是相同的多项式，那么多项式快速幂即可。
#include<cstdio>
#include<cctype>
#include<queue>
#include<cstring>
#include<algorithm>
#define rep(i,s,t) for(int i=s;i<=t;i++)
#define dwn(i,s,t) for(int i=s;i>=t;i--)
#define ren for(int i=first[x];i;i=next[i])
using namespace std;
const int BufferSize=1<<16;
char buffer[BufferSize],*head,*tail;
inline char Getchar() {
if(head==tail) {
int l=fread(buffer,1,BufferSize,stdin);
tail=(head=buffer)+l;
}
return *head++;
}
inline int read() {
int x=0,f=1;char c=Getchar();
for(;!isdigit(c);c=Getchar()) if(c=='-') f=-1;
for(;isdigit(c);c=Getchar()) x=x*10+c-'0';
return x*f;
}
const int maxn=18010;
const int p=1004535809;
const int G=3;
const int NUM=15;
typedef long long ll;
ll wn[maxn],inv;
ll pow(ll n,ll m,ll mod=p) {
ll ans=1;
for(;m;m>>=1,(n*=n)%=mod) if(m&1) (ans*=n)%=mod;
return ans;
}
void NTT(ll* A,int len,int tp) {
int j=len>>1,c=0;
rep(i,1,len-2) {
if(i<j) swap(A[i],A[j]);int k=len>>1;
while(j>=k) j-=k,k>>=1;j+=k;
}
for(int i=2;i<=len;i<<=1) {
c++;
for(int j=0;j<len;j+=i) {
ll w=1;
for(int k=j;k<j+(i>>1);k++) {
ll u=A[k],t=w*A[k+(i>>1)]%p;
A[k]=(u+t)%p;A[k+(i>>1)]=(u-t+p)%p;
(w*=wn[c])%=p;
}
}
}
if(tp<0) {
rep(i,1,len/2-1) swap(A[i],A[len-i]);
ll inv=pow(len,p-2);
rep(i,0,len-1) (A[i]*=inv)%=p;
}
}
int check(int g,int m) {
for(int i=2;i*i<m;i++) if((m-1)%i==0&&(pow(g,i,m)==1||pow(g,(m-1)/i,m)==1)) return 0;
return 1;
}
int n,m,X,S,len,gs=2,c[maxn];
ll D[maxn];
void mul(ll* A,ll* B) {
rep(i,0,len-1) D[i]=B[i];
NTT(A,len,1);NTT(D,len,1);
rep(i,0,len-1) (A[i]*=D[i])%=p;
NTT(A,len,-1);
dwn(i,len-1,m-1) (A[i-m+1]+=A[i])%=p,A[i]=0;
}
ll A[maxn],B[maxn];
int main() {
rep(i,0,NUM-1) wn[i]=pow(G,(p-1)/(1<<i));
n=read();m=read();X=read();S=read();
len=1;while(len<=(m<<1)) len<<=1;
while(!check(gs,m)) gs++;
int m0=1;c[1]=0;
rep(i,1,m-2) (m0*=gs)%=m,c[m0]=i;
rep(i,1,S) {
int x=read()%m;
if(x) A[c[x]]=1;
}
B[0]=1;
for(;n;mul(A,A),n>>=1) if(n&1) mul(B,A);
printf("%lld
",B[c[X]]);
return 0;
}