【SNOI2017】【洛谷P5366】【BZOJ5019】【LOJ2257】遗失的答案（数论基础）（状压DP）（FMT）

LOJ传送门

---

题解：

（超大常数在线丢人）

其实可以不用脑子先把非法判掉。

实际上不同的质因子非常少。

那么实际上我们就是要求选择若干个数，其中存在指数达到$gcd$的下界，且存在指数达到$lcm$的上界。

发现是个带位或转移的状压DP。

FMT算出DP的前后缀就能支持强制令某个数出现，随便水了。

---

代码：

#include<bits/stdc++.h>
#define ll long long
#define re register
#define gc get_char
#define cs const

namespace IO{
	inline char get_char(){
		static cs int Rlen=1<<22|1;
		static char buf[Rlen],*p1,*p2;
		return (p1==p2)&&(p2=(p1=buf)+fread(buf,1,Rlen,stdin),p1==p2)?EOF:*p1++;
	}
	
	template<typename T>
	inline T get(){
		char c;
		while(!isdigit(c=gc()));T num=c^48;
		while(isdigit(c=gc()))num=(num+(num<<2)<<1)+(c^48);
		return num;
	}
	inline int getint(){return get<int>();}
}
using namespace IO;

using std::cerr;
using std::cout;

cs int mod=1e9+7;
inline int add(int a,int b){return (a+=b)>=mod?a-mod:a;}
inline int dec(int a,int b){return (a-=b)<0?a+mod:a;}
inline int mul(int a,int b){static ll r;r=(ll)a*b;return r>=mod?r%mod:r;}
inline int power(int a,int b,int res=1){
	for(;b;b>>=1,a=mul(a,a))(b&1)&&(res=mul(res,a));
	return res;
}
inline void Inc(int &a,int b){(a+=b)>=mod&&(a-=mod);}
inline void Dec(int &a,int b){(a-=b)<0&&(a+=mod);}
inline void Mul(int &a,int b){a=mul(a,b);}

cs int N=650,SIZE=1<<16|1; 

int n,G,L,Q;

int S; 

inline void FMT(int *A){
	for(int re i=1;i<S;i<<=1)
	for(int re j=0;j<S;j+=i<<1)
	for(int re k=0;k<i;++k)Inc(A[i|j|k],A[j|k]);
}
inline void IFMT(int *A){
	for(int re i=1;i<S;i<<=1)
	for(int re j=0;j<S;j+=i<<1)
	for(int re k=0;k<i;++k)Dec(A[i|j|k],A[j|k]);
}

int p[N],pcnt,tim[N];

inline void get_factor(int x){
	for(int re i=2;(ll)i*i<=x;++i)if(x%i==0){
		p[++pcnt]=i;
		while(x%i==0)x/=i,++tim[pcnt];
	}
	if(x>1){
		p[++pcnt]=x;
		tim[pcnt]=1;
	}
}

int s[SIZE],t[SIZE],f[N][SIZE],g[N][SIZE],h[N][SIZE];
std::map<int,int> id;
void dfs(int now,int cur,int state){
	if(now>pcnt){
		++s[state];id[cur]=state;
		return ;
	}
	for(int re i=0;i<=tim[now];++i){
		dfs(now+1,cur,state|((i==0)<<(now-1))|((i==tim[now])<<(now+pcnt-1)));
		if((cur*=p[now])>n)return ;
	}
}

int num[SIZE],ans[SIZE],cnt;

signed main(){
#ifdef zxyoi
	freopen("lost.in","r",stdin);
#endif
	n=getint(),G=getint(),L=getint(),Q=getint();
	if(L%G){while(Q--)puts("0");return 0;}
	n/=G,L/=G;get_factor(L);
	dfs(1,1,0);S=1<<(pcnt<<1);
	for(int re i=0;i<S;++i)if(s[i]){
		num[++cnt]=i;
		t[cnt]=power(2,s[i])-1;
	}
	f[0][0]=1,g[cnt+1][0]=1;
	for(int re i=0;i<cnt;++i)
	for(int re j=0;j<S;++j){
		if(!f[i][j])continue;
		Inc(f[i+1][j],f[i][j]);
		Inc(f[i+1][j|num[i+1]],mul(f[i][j],t[i+1]));
	}
	for(int re i=cnt+1;i>1;--i)
	for(int re j=0;j<S;++j){
		if(!g[i][j])continue;
		Inc(g[i-1][j],g[i][j]);
		Inc(g[i-1][j|num[i-1]],mul(g[i][j],t[i-1]));
	}
	for(int re i=0;i<=cnt+1;++i)FMT(f[i]),FMT(g[i]);
	for(int re i=1;i<=cnt;++i){
		for(int re j=0;j<S;++j)Inc(h[i][j],mul(f[i-1][j],g[i+1][j]));
		IFMT(h[i]);
		for(int re j=0;j<S;++j)if((j|num[i])==S-1)Inc(ans[i],h[i][j]);
		Mul(ans[i],power(2,s[num[i]]-1));
	}
	while(Q--){
		int x=getint();
		if(x%G){cout<<"0\n";continue;}x/=G;
		if(L%x||x>n){cout<<"0\n";continue;}
		int y=std::lower_bound(num+1,num+cnt+1,id[x])-num;
		cout<<ans[y]<<"\n";
	}
	return 0;
}