Link
https://www.lydsy.com/JudgeOnline/problem.php?id=3992
題意
給定一個集合 , 中元素均爲小於 的非負整數。( 爲質數)
問:用集合中的數能夠構造出多少個不同的數列(允許多次使用同一個數),
滿足數列中所有數的乘積 ?
答案模 。
設計:填了前 個元素,乘積 的方案數
轉移
注意到 是質數就算了, ……幹什麼用的?奇素數??
那先不管具體怎麼樣,我們思考用用原根
表示填前 個元素,乘積
設
(實際上 )
但是這樣的話
所以我們應該設
即
那麼
注意在卷 卷卷的時候要不停把下標饃饃
#include<cstdio>
#include<iostream>
#include<algorithm>
#include<cstdlib>
#include<cmath>
#include<cstring>
#include<cctype>
#include<ctime>
using namespace std;
#define getchar() (frS==frT&&(frT=(frS=frBB)+fread(frBB,1,1<<12,stdin),frS==frT)?EOF:*frS++)
char frBB[1<<12], *frS=frBB, *frT=frBB;
inline void read(int& x)
{
x = 0; char ch = getchar(); bool w = 0;
while (!isdigit(ch)) w |= (ch == '-'), ch = getchar();
while (isdigit(ch)) x = x * 10 + (ch ^ 48), ch = getchar();
w ? (x = -x) : 0;
}
const int MAXN = 16384+16;
const long long p = 1004535809;
inline long long qpow(long long a, long long b, const long long& p)
{
long long ret = 1;
a %= p;
while(b)
{
if (b & 1)
{
ret *= a;
ret %= p;
}
a *= a;
a %= p;
b >>= 1;
}
return ret;
}
const long long g = 3;
int n, m, x, s, Lim = 1, Log;
long long q;
long long Wn[MAXN][2];
int Rev[MAXN];
inline long long Adjust(long long x)
{
x = (x % p) + p;
return (x >= p) ? (x - p) : x;
}
void NTT(long long *a, const bool& Type = 0)
{
for (register int i = 0; i < Lim; ++i) if (Rev[i] > i) swap(a[Rev[i]], a[i]);
register long long x;
for (register int Len, Mid = 1, qwq; Mid < Lim; Mid <<= 1)
{
Len = Mid << 1;
qwq = Lim / Len;
for (register int Pos = 0; Pos < Lim; Pos += Len)
{
for (register int Sub = 0; Sub < Mid; ++Sub)
{
x = Wn[qwq * Sub][Type] * a[Pos + Mid + Sub] % p;
a[Pos + Mid + Sub] = Adjust(a[Pos + Sub] - x);
a[Pos + Sub] = (a[Pos + Sub] + x) % p;
}
}
}
}
int gM;
int PMStack[8005], Ind[8005];
inline void PMSplit(int x)
{
for (register int i = 2; i * i <= x; ++i)
{
if (x % i == 0)
{
PMStack[++PMStack[0]] = i;
while (x % i == 0) x /= i;
}
}
}
inline int GetPrimeRoot(int x, int phi)
{
PMSplit(phi);
register bool flag;
for (register int i = 2; i < x; ++i)
{
flag = 0;
for (register int j = 1; j <= PMStack[0]; ++j)
{
if (qpow(i, phi/PMStack[j], m) == 1)
{
flag = 1;
break;
}
}
if (!flag) return i;
}
}
long long a[MAXN];
long long ta[MAXN], tb[MAXN], inv;
inline void PolyMul(long long *x, long long *y)
{
for (register int i = 0; i < Lim; ++i) ta[i] = x[i] % p;
for (register int i = 0; i < Lim; ++i) tb[i] = y[i] % p;
NTT(ta); NTT(tb);
for (register int i = 0; i < Lim; ++i) ta[i] = ta[i] * tb[i] % p;
NTT(ta, 1);
for (register int i = 0; i < Lim; ++i) ta[i] = ta[i] * inv % p;
for (register int i = 0; i < m-1; ++i) x[i] = (ta[i] + ta[i+m-1]) % p;
//同樣要小心自乘的情況
}
long long Ans[MAXN];
void PolyPow(long long *a, int Index)
{
Ans[0] = 1;
while (Index)
{
if (Index & 1)
{
PolyMul(Ans, a);
}
PolyMul(a, a);
Index >>= 1;
}
}
int main()
{
read(n), read(m), read(x), read(s);
gM = GetPrimeRoot(m, m-1);
for (register int i = 0, cur = 1; i < m - 1; ++i) Ind[cur] = i, cur *= gM, cur %= m;
for (register int t, i = 1; i <= s; ++i) read(t), t ? (a[Ind[t]] = 1) : 0;
int mm = m << 1;
while (Lim <= mm) Lim <<= 1, ++Log;
for (register int i = 0; i < Lim; ++i) Rev[i] = (Rev[i>>1]>>1)|((i&1)<<(Log-1));
inv = qpow(Lim, p-2, p);
q = (p - 1) / Lim;
register long long fafa = qpow(g, q, p);
Wn[0][0] = Wn[Lim][0] = Wn[0][1] = Wn[Lim][1] = 1;
for (register int i = 1; i < Lim; ++i) Wn[i][0] = Wn[i-1][0] * fafa % p;
for (register int i = 1; i < Lim; ++i) Wn[i][1] = Wn[Lim-i][0];
PolyPow(a, n);
printf("%lld", Ans[Ind[x]]);
return 0;
}