bzoj3992

题解：

求模素数 p 原根的方法：对 p-1 进行素因子分解，记pi为p-1的第i个因子，若恒有a^((p-1)/pi) mod p ≠ 1  成立，则 a 就是 p 的原根。（对于合数求原根，只需把 p-1 换成 phi(p) 即可）

首先比较暴力是f[i][j]表示前i个，乘积为j

然后是n*m^2的

我们可以利用原根将每个数变成g^x

这样就变成了加和

可以利用生成函数的思想

然后分治+fft就可以了

代码：

#include <bits/stdc++.h>
using namespace std;
#define rint register int
#define rep(i,h,t) for(rint i=h;i<=t;i++)
#define dep(i,t,h) for(rint i=t;i>=h;i--)
#define IL inline
const double pi=acos(-1.0);
const int N=3e4;
int n,m,l,r[N],a[N],b[N],w[N],c[N],g,kk,G=;
int p=;
IL int fsp(int x,int y,int p)
{
  int ans=;
  while (y)
  {
    if(y&) ans=1ll*ans*x%p;
    x=1ll*x*x%p;
    y>>=;
  }
  return ans;
}
IL void jf(int &x,int y)
{
  x+=y;
  if (x>p) x-=p;
}
IL void fft(int *a,int o)
{
  rep(i,,n-) if (i>r[i]) swap(a[i],a[r[i]]);
  for (int i=;i<n;i*=)
  {
    int wn=fsp(G,(p-)/(i*),p); w[]=;
    rep(j,,i-) w[j]=(1ll*w[j-]*wn)%p;
    for (int j=;j<n;j+=(i*))
    {
      int *x=a+j,*y=a+i+j;
      for (rint k=;k<i;k++)
      {
        const int t=(1ll*w[k]*y[k])%p;
        y[k]=x[k]-t; if (y[k]<) y[k]+=p;
        x[k]+=t; if (x[k]>p) x[k]-=p;
      }
    }
  }
  if (o==-)
  {
    reverse(&a[],&a[n]);
    for (int i=,inv=fsp(n,p-,p);i<n;i++)
      a[i]=1ll*a[i]*inv%p;
    rep(i,kk,m)
      jf(a[i%kk],a[i]),a[i]=;
  }
}
IL void query()
{
  l=;
  for (n=;n<=m;n<<=) l++;
  rep(i,,n-) r[i]=(r[i/]/)|((i&)<<(l-));
  fft(a,); fft(b,);
  rep(i,,n-) a[i]=1ll*a[i]*b[i]%p;
  fft(a,-);
}
IL void get_g(int x)
{
  rep(i,,x-)
  {
    int j;
    for (j=;j<m;j++)
      if (fsp(i,j,m)==) break;
    if (j==m-)
    {
      g=i; break;
    }
  }
}
struct re{
  int a[N];
}now;
re fsp2(int x)
{
  if (x==) return(now);
  re tmp=fsp2(x/);
  memcpy(b,tmp.a,sizeof(tmp.a));
  memcpy(a,tmp.a,sizeof(tmp.a));
  query();
  if (x%)
  {
    memcpy(b,c,sizeof(c));
    query();
  }
  memcpy(tmp.a,a,sizeof(a));
  return(tmp);
}
int f[N];
int main()
{
  freopen("1.in","r",stdin);
  freopen("1.out","w",stdout);
  ios::sync_with_stdio(false);
  int x,s;
  cin>>n>>m>>x>>s; kk=m-;
  rep(i,,s) cin>>a[i];
  get_g(m);
  rep(i,,m-) f[fsp(g,i,m)]=i;
  m*=;
  rep(i,,s) if (a[i]) now.a[f[a[i]]]++,c[f[a[i]]]++;
  fsp2(n);
  cout<<(a[f[x]]+p)%p;
  return ;
}