标签：P8292 省选 ll sqrt 联考 leq sum include 质数

正题

题目链接:https://www.luogu.com.cn/problem/P8292

---

题目大意

有\(n\)张卡牌，第\(i\)张上的数字是\(s_i\)。\(m\)次询问给出\(c_i\)个质数，要求选择一些卡使得这些卡的乘积是这些质数的倍数，求方案数。

\(1\leq n\leq 10^6,1\leq s_i\leq 2000,1\leq m\leq 1500,\sum_{i=1}^m c_i\leq 18000\)

---

解题思路

对于一个数字\(x\)来说，大于\(\sqrt x\)的质因子只会有一个，也就是对于大于所有的\(\sqrt s_i\)的质数\(p\)，它们的倍数之间不会产生重复，可以分开计数。

而小于\(\sqrt s_i\)的质数个数只有\(14\)个，我们可以考虑状压。

处理出\(g_s\)表示在集合\(s\)中的质数都不选择（即没有这些数的倍数）时的数字个数，然后再预处理出一个\(f_{i,s}\)表示质数\(i\)的倍数中集合\(s\)中的质数都不选择时的数字个数。

那么对于一个询问，我们先处理出小于等于\(\sqrt s\)的质数被加进去的状态\(S\)，然后对于大于\(\sqrt s\)的每个质数\(p\)，那么答案就是

\[\sum_{T\sube S}\left(2^{g_T}\times \sum_{p}\left(1-\frac{1}{2^{f_{p,T}}}\right)\right)(-1)^{|T|} \]

然后求和就可以了

---

code

#include<cstdio>
#include<cstring>
#include<algorithm>
#include<vector>
#include<cctype>
#define ll long long
using namespace std;
const ll N=2100,M=1e6+10,L=14,P=998244353;
ll n,m,cnt,pri[320],pw[M],iw[M],ct[1<<L];
ll w[N],f[1<<L][320],g[1<<L],q[N*10],v[N];
vector<ll> r;bool usd[N];
ll read(){
	ll x=0,f=1;char c=getchar();
	while(!isdigit(c)){if(c=='-')f=-f;c=getchar();}
	while(isdigit(c)){x=(x<<1)+(x<<3)+c-48;c=getchar();}
	return x*f;
}
void Prime(){
	for(ll i=2;i<N;i++){
		if(!v[i])pri[cnt]=i,cnt++,v[i]=cnt-1;
		for(ll j=0;j<cnt&&i*pri[j]<N;j++){
			v[i*pri[j]]=1;
			if(i%pri[j]==0)break;
		}
	}
	return;
}
signed main()
{
	Prime();
	n=read();pw[0]=iw[0]=1;
	for(ll i=1;i<=n;i++)pw[i]=pw[i-1]*2ll%P;
	for(ll i=1;i<=n;i++)iw[i]=iw[i-1]*(P+1)/2ll%P;
	for(ll i=1,x;i<=n;i++)x=read(),w[x]++;
	ll MS=(1<<L);
	for(ll s=1;s<MS;s++)ct[s]=ct[s-(s&-s)]+1;
	for(ll s=0;s<MS;s++){
		r.clear();
		memset(usd,0,sizeof(usd));
		for(ll i=0;i<L;i++)
			if((s>>i)&1){
				for(ll j=pri[i];j<=2000;j+=pri[i])
					usd[j]=1;
			}
		for(ll i=1;i<=2000;i++)
			if(!usd[i])g[s]+=w[i],r.push_back(i);
		g[s]=pw[g[s]];
		for(ll i=L;i<cnt;i++){
			f[s][i]=0;
			for(ll j=0;r[j]*pri[i]<=2000;j++)
				f[s][i]+=w[r[j]*pri[i]];
			ll k=f[s][i];
			f[s][i]=iw[k]*(pw[k]-1)%P;
		}
	}
	m=read();
	while(m--){
		ll k=read();
		ll t=0,tot=0,ans=0;
		for(ll i=1,x;i<=k;i++){
			x=read();
			if(v[x]<L)t|=1<<v[x];
			else q[++tot]=v[x];
		}
		sort(q+1,q+1+tot);
		tot=unique(q+1,q+1+tot)-q-1;
		for(ll s=t;s>=0;s=(s-1)&t){
			ll k=(ct[s]&1)?(P-g[s]):g[s];
			for(ll i=1;i<=tot;i++)
				k=k*f[s][q[i]]%P;
			(ans+=k)%=P;
			if(!s)break;
		}
		printf("%lld\n",(ans+P)%P);
	}
	return 0;
}

标签：P8292,省选,ll,sqrt,联考,leq,sum,include,质数
来源： https://www.cnblogs.com/QuantAsk/p/16191815.html

本站声明： 1. iCode9 技术分享网（下文简称本站）提供的所有内容，仅供技术学习、探讨和分享；
2. 关于本站的所有留言、评论、转载及引用，纯属内容发起人的个人观点，与本站观点和立场无关；
3. 关于本站的所有言论和文字，纯属内容发起人的个人观点，与本站观点和立场无关；
4. 本站文章均是网友提供，不完全保证技术分享内容的完整性、准确性、时效性、风险性和版权归属；如您发现该文章侵犯了您的权益，可联系我们第一时间进行删除；
5. 本站为非盈利性的个人网站，所有内容不会用来进行牟利，也不会利用任何形式的广告来间接获益，纯粹是为了广大技术爱好者提供技术内容和技术思想的分享性交流网站。