Description
Solution
首先 \\(v_1,v_2,v_3...v_n,P\\) 能够构成的最小数是 \\(gcd(P,v_1,v_2,v_3...v_n)\\)
然后 \\(gcd(P,v_1,v_2,v_3...v_n)|w_i\\) 则可以构成 \\(w_i\\)
所以我们直接背包一下就可以了,设 \\(m\\) 为 \\(P\\) 的约数个数,\\(m\\) 最多是 \\(n^{\\frac{1}{3}}\\)
那么复杂度就是 \\(O(n*m*logP)\\)
容易发现如果 \\(gcd(v_i,P)=gcd(v_j,P)\\) ,那么 \\(v_i,v_j\\) 对答案的影响是一样的,把相同类型的一起处理,把答案乘上 \\(2^{cnt}-1\\) 就行了
复杂度 \\(O(m^2*logP)\\)
#include<bits/stdc++.h>
using namespace std;
template<class T>void gi(T &x){
int f;char c;
for(f=1,c=getchar();c<\'0\'||c>\'9\';c=getchar())if(c==\'-\')f=-1;
for(x=0;c<=\'9\'&&c>=\'0\';c=getchar())x=x*10+(c&15);x*=f;
}
const int N=1e6+10,mod=1e9+7;
int n,Q,P,num=0,a[N],m=0,c[N],f[2010][2010],bin[N],ans[N];
map<int,int>id;
inline void priwork(){
for(int i=1;i*i<=P;i++)
if(P%i==0){
a[++m]=i;
if(i*i!=P)a[++m]=P/i;
}
sort(a+1,a+m+1);
for(int i=1;i<=m;i++)id[a[i]]=i;
}
inline int gcd(int x,int y){return y?gcd(y,x%y):x;}
int main(){
freopen("pp.in","r",stdin);
freopen("pp.out","w",stdout);
cin>>n>>Q>>P;
priwork();
int x;bin[0]=1;
for(int i=1;i<=n;i++)gi(x),c[id[gcd(x,P)]]++,bin[i]=bin[i-1]*2%mod;
f[0][m]=1;
for(int i=1;i<=m;i++){
for(int j=1;j<=m;j++){
if(c[i]){
int p=id[gcd(a[i],a[j])];
f[i][p]=(f[i][p]+1ll*f[i-1][j]*(bin[c[i]]-1))%mod;
}
f[i][j]=(f[i][j]+f[i-1][j])%mod;
}
}
for(int i=1;i<=m;i++)
for(int j=1;j<=i;j++)
if(a[i]%a[j]==0)ans[i]=(ans[i]+f[m][j])%mod;
while(Q--){
gi(x);
printf("%d\\n",ans[id[gcd(x,P)]]);
}
return 0;
}